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Da bi se rezultati eksperimentalnoga rada u poljoprivredi mogli 
pravilno interpretirati i ocijeniti prije nego budu predani poljoprivrednoj 
praksi potrebno je, da se brojčani podaci ispravno matematski obrade. 
U tu nara svrhu služi biomelrika. 

God. 1929 napisao sam naputak za izračunavanje raznih vari- 
jaciono statističkih vrijednosti pod naslovom „ Variaciona - statistika u 
eksperimentalnoj poljoprivredi 1 . Knjižica je izašla iz štampe bez mog 
imprimatura na korekture sa mnogim štamparskim pogreškama, ali je 
i ovakova ubrzo bila rasprodana. Budući da nam je strana literatura 
o biometrici slabo pristupačna i da su i tamo pojedine partije ove 
discipline većinom podijeljene po različitim specijalnim priručnicima, a 
za primjere uzeti ponajviše organizmi i svojstva, koja nemaju poljo¬ 
privredne važnosti, odlučio sam da za naše potrebe pretežno na bili- 
nogojskom materijalu iz vlastitih pokusa prikažem primienu različnih 
biometričkih metoda. 

Knjiga je namijenjena u prvom redu poljoprivrednim stručnjacima 
i studentima poljoprivrede ali držim, da će se njome moći koristiti 
i stručnjaci i studenti drugih bioloških disciplina. 

Ona naročito ima da bude poticaj za što širu i što češću primjenu 
biometrike u eksperimentalnom radu na cijelom području prirodoslovnih 
nauka kod nas. 

Toplo zahvaljujem uredniku Znanstvene poljoprivredne knjižnice 
kolegi univ. prof. Dr. Mihovilu Gračaninu za mnoge dragocjene sugestije 
u pogledu pedagoške i stilističke obrade sastavaka. 

Za savjesnu kontrolu računskih operacija i štamparske korekture 
zahvaljujem svom asistentu Ing. Mariji Kump. 

Zagreb, svibnja 1946. 



UVOD 


Intenzivnim studijem fizikalnih i kemijskih svojstava mrtve 
tvari — naročito tokom zadnjeg rata — došlo se do mnogih 
novih otkrića, koja ljudski um kao nadmoćan gospodar može 
brzo iskoristiti i primijeniti u svim granama tehnike. Mrtvu 
tvar on modelira i transformira po svojoj volji u različne agre¬ 
gate podložene zakonima fizike i kemije. Tamo pak, gdje po¬ 
kušava promijeniti kvalitet i kvantitet proizvodnje živih bića, 
kao na pr. u poljoprivredi, nailazi na vrlo velike zapreke i tek 
polagano, uz nebrojene napore i materijalne žrtve, krči put i 
nastoji, da iskoristivši poznate faktore proizvodnje i različne 
nove biološke procese, stvori organizme, koji će u određenom 
proizvodnom području dati što veću gospodarsku korist. 

Vanjski izgled biljke i životinje rezultanta je genetskih oso¬ 
bina individua i njegove sposobnosti, kako on reagira na vanjske 
i unutarnje faktore proizvodnje. 

Točno djelovanje ovih faktora, ni pojedinačno ni u skupi¬ 
nama, nije još istraženo. Intenzivnim eksperimentalnim radom 
na tom području, uz eksaktnu obradu brojčanih podataka, do¬ 
lazimo do pravilnih dedukcija i bržih praktičnih uspjeha. 

Opažajući i mjereći pojedina svojstva na skupini genetski 
jednakih organizama ustanovit ćemo da ona nisu posve jednaka 
nego različita. Ova je razlika uzrokovana time, što se intenzitet 
pojedinih vegetacijskih faktora ne može udesiti tako, da bi bio 
svim organizmima, pa ni samo jednoj maloj skupini, u jednakoj 
mjeri pristupačan. Prema tome je nejednaki intenzitet pro¬ 
izvodnih faktora uzrok, da svako svojstvo organizma u izvje¬ 
snim granicama varira. 

G a u s s je prvi ustanovio, da kod mjerenja, vaganja ili 
brojenja pravimo redovno izvjesne pogreške, i da se brojčana 
vrijednost tog svojstva dade najbolje ustanoviti pomoću aritme¬ 
tičke sredine velikog broja pojedinačnih mjerenja. Taj je 
Gaussov zakon o pogrešci bio najprije primjenjivan u astrono¬ 
miji, pa zatim u fizici i tehnici. 
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Belgijski antropolog Q u e t e 1 e t je bio prvi, koji je u 
svom djelu »Anthropometric«, iz god. 1871. utvrdio, da Gaussov 
zakon o pogrešci vrijedi i kod mjerenja različnih svojstava živih 
organizama. Queteletova su iskustva u brojčanom obra¬ 
đivanju svojstava kod čovjeka potakla mnoge druge prirodo¬ 
slovce, da matematički obrade brojčane podatke dobivene svo¬ 
jim pokusima s biljkama i životinjama. Iz tih prvih nastojanja, 
da se matematika primijeni na istraživanje varijabiliteta poje¬ 
dinih svojstava organizama, razvila se kasnije posebna nauka, 
biometrika, naročito zaslugom engleskog matematičara 
Pearsona, koji je god. 1900. izdao znamenito djelo »T h e 
grammar of science«. 

Pod bio metrikom dakle razumijevamo primjenu ma¬ 
tematičkih metoda kod istraživanja raznolikosti organizma i 
međusobne ovisnosti pojedinih svojstava, kako bi se ustanovio 
stupanj signifikantnosti eksperimentom dobivenih brojčanih 
podataka. 

Račun vjerojatnosti i kombinatorika su matematički te¬ 
melji, na kojima je izgrađena biometrika. 

Ovdje ćemo se pozabaviti naročito biometrikom varijacije 
ili varijacionom statistikom i biometrikom korelacije kod poljo¬ 
privrednih biljaka i domaćih životinja, da upoznamo metodiku 
u poljoprivrednim istraživanjima. 


A. VARIJABILITET 


a) Općenito 

Pod varijabilitetom nekog organizma razumije¬ 
vamo promjenu njegovih svojstava izazvanu vanjskim i unutar¬ 
njim faktorima života. 

Ako promotrimo potomke istih roditelja, ne ćemo nikada 
naći ni dva individua potpuno jednaka. Uzroci nejednakosti, 
odnosno varijabiliteta svojstava uglavnom su dvojaki i to: 

a) roditelji u svojim nasljednim osnovama nisu bili jednaki, 
ili 

b) roditelji su imali podjednake nasljedne osnove, a poje¬ 
dini su potomci rasli i razvijali se pod različnim životnim 
uvjetima. 

Svako je svojstvo organizma uvjetovano djelovanjem i vanj¬ 
skih i unutarnjih (nasljednih i dr.) faktora života. 

Skup svih nasljednih osnova sačinjava njegov idiotip. Ove 
se osnove nalaze u: a) staničnoj jezgri, b) c i to¬ 
pi a z m i i c) p 1 a s t i d i m a. Prema tome se idiotip sa¬ 
stoji od: 

a) genotipa: osnove lokalizirane u staničnoj jezgri, 

b) plazmotipa: „ „ „ citoplazmi, i 

c) plastidotipa: „ „ „ plastidima. 

Najvažnija su nasljedna svojstva organizma uvjetovana 

genotipom. 

G e n o t i p sačinjava velik broj temeljnih nasljedstvenih 
jedinica, koje zovemo nasljedstvenim ili genetskim fakto¬ 
ri m a ili g e n i m a. 

Geni određuju osobine svakog organizma, a vanjski 
faktori odlučuju, do koje će se mjere pojedina svojstva oči¬ 
tovati. 

Ako imamo dva individua jednakoga genotipa, od kojih se 
jedan razvija pod jednim, a drugi pod drugim životnim uvje- 
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tima, onda njihove osobine ne će i ne mogu biti istovjetne. 
Konačni izgled izvjesnog svojstva je rezultanta međusobne dje¬ 
latnosti gena i vanjskih faktora. 

^JSit nasljeđivanja leži u tome, da se roditeljski geni preko 
rasplodnih stanica prenose na potomstvo. 

Da se neko svojstvo očituje, nije dovoljno prisustvo samih 
gena, nego je potreban i stanoviti intenzitet vanjskih faktora 
života. Već prema intenzitetu vanjskih faktora jedan će se isti 
organizam različito razvijati. 

Skup svih gena jednog organizma sačinjavaju, kako je već 
spomenuto, njegov g e n o t i p, a skup svih značajka, svojstava 
i sposobnosti t. j. nasljednih i nenasljednih, njegov prividni tip 
ili f e n o t i p. 

Zbog modifikacije istome fenotipu mogu pripadati i različni 
genotipovi, a i individui istoga genotipa ne moraju pripadati 
istome fenotipu. 

Od genotipa zavise granice, do kojih se neko svojstvo može 
promijeniti — varirati. Ove su granice uvjetovane reakcionom 
sposobnošću ili reakcionom normom genotipa. Ova je reakciona 
norma nasljedna. 

Po samome svojstvu, dakle po fenotipu, obično još ne 
možemo znati, kakav je njegov genotip. Naročito to vrijedi za 
kvantitativna svojstva, koja jače variraju od kvalitativnih. 

i/ Da bismo mogli brojkama karakterizirati izvjesno svojstvo, 
treba ga proučiti varijaciono-statistički i to tako, da ga istra¬ 
žimo na prosječnom uzorku i na dovoljno velikom broju orga¬ 
nizama. I onda, ako za takva istraživanja upotrebimo skup 
individua istog genotipa, naravno je, da će obzirom na istra¬ 
ženo svojstvo među njima biti i individua sa ekstremno veli¬ 
kim vrijednostima — plus varijante — i ekstremno ma¬ 
lim — minus varijante. Osim toga naći ćemo pojedine 
varijante, koje imaju različne gradacije između jednog i drugog 
ekstrema. Raz mak između jednog i drugog ekstrem a zovemo 
var i ja ci ono m . ili m o d i f i k a c i onom š i r inom. 
Svaki individuum ima u toj širini, već prema brojčanoj vrijed¬ 
nosti svog svojstva, posve određeno mjesto. U toj varijacionoj 
širini svaki individuum se očituje kao jedna varijanta. 

Ako brojčanu vrijednost pojedinih varijanata svrstamo u 
skupine s jednakim razmakom, a te skupine opet svrstamo po 
veličini, od najmanje do najveće, dobit ćemo vrlo pregleđnu i 
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karakterističnu podjelu varijanata. Vidjeti ćemo, da se broj 
varijanata u pojedinim skupinama diže od minus ekstrema 
otprilike do sredine, a onda opet, više ili manje simetrično pada 
do plus ekstrema.jjMajviše varijanata nalazimo dakle u srednjim 
skupinama.^ 

Uzrok, što kod pojedinih varijanata, koje pripadaju jednom 
istom genotipu, isto svojstvo nije jednako nego različito — to 
jest da varira — leži u tome, što intenzitet vanjskih 
faktora (na pr. kod bilina: svjetlost, toplina, vlaga, hranjive 
tvari i t. d.) nije i ne može biti za svaku jedinicu potpuno jednak. 
Osim toga isti faktori djeluju nekad povoljno, nekad nepo¬ 
voljno već prema svom intenzitetu. 


Tab. 1. 


Kombinacije, 

njihova vrijednost i frekvencija 




a b c d 


0 


1 




1 
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Na pr.: neka visina organizma zavisi o 4 vegetacijska 
faktora A, B, C i D, od kojih svaki vrijedi 1. Odsustvo svakoga 
od njih označit ćemo s a, odnosno b, c i d, a odsustvo svih četi¬ 
riju s 0. U ovom slučaju ima 16 različitih kombinacija, a nji¬ 
hova brojčana vrijednost varira od 0 do 4, kako to vidimo na 
skrižaljci Tab. I. 

Ova podjela odgovara raspodjeli dobivenoj po binominal- 
noj formuli: (a + b) 4 , gdje je a = b = 1. 

Ako na abscisi označimo 5 jednako udaljenih točaka i na 
njih nanesemo frekvenciju, i onda te točke frekvencije među¬ 
sobno povežemo krivuljom, dobivamo modifikacionu 
ili varijacionu krivulju. (SI. 1). 



SI. 1. 

l/ Što je veći broj individua ili varijanata, ta se krivulja sve 
to više približuje binominalnoj krivulji konstruiranoj iz for¬ 
mule (a + b) n gdje je a = b = 1 dakle: 


(a + b) 1 = a+b . 1 + 1 

(a+b) 2 = a 2 + 2ab + b 2 . 1+2+1 

(a+b) 3 = a 2 + 3a 2 b+3ab 2 +b 3 . 1 + 3+3+1 

(a+b) 4 =a 4 + 4a 3 b + 6a 2 b 2 + 4ab 3 +b 4 .1+4 + 6 + 4+1 


Kad je n = oo, onda varijaciona krivulja prelazi u n o r- 
malnu, Gaussovu, Queteletovu, idealnu ili b i- 


-ir 



nominalnu krivulju (SI. 2.). Ovu krivulju možemo kon¬ 
struirati po formuli: y ~~ 


X 1 

k C 2 aa 



y je u našem slučaju visina krivulje kod udaljenosti x od 
sredine krivulje; a je konstanta za širinu krivulje. 

Budući da se utjecajem vanjskih i unutarnjih faktora života 
pojedina svojstva organizma mijenjaju u izvjesnim granicama, 
potrebno je, da bismo mogli procijeniti ta svojstva, naći njihovu 
prosječnu vrijednost ili aritmetičku sredinu 
i granice njihova varijabiliteta. 

U našem radu i nastojanjima oko unapređenja poljopri¬ 
vrede naći ćemo se pred čitavim nizom takvih problema, koje 
treba varijaciono-statistički obraditi i proučiti i to na području 
biljne proizvodnje, živinogojstva, gospodarske uprave i t. d. 

Treba na pr. pronaći, koja domaća ili oplemenjena sorta 
neke vrste biljke najbolje odgovara izvjesnim klimatskim i 
edafskim prilikama, ili koje su pasmine pojedinih vrsta doma¬ 
ćih životinja po kvaliteti i kvantiteti svojih produkata najbolje 
za izvjesno proizvodno područje ili koja je obrada tla najbolja, 
kako treba gnojiti izvjesno polje, koja je najbolja dubina i gu¬ 
stoća sjetve, kakav je plodored najbolji, odnosno najbolja obrada 
tla za vrijeme rasta pojedinih vrsti i sorata bilja, kakvo je drža- 
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nje i kakva je hranidba pojedinih pasmina domaćih životinja 
najbolja i t. d. 

Želimo li pokusom dobiti valjane zaključke, treba primije¬ 
niti ispravnu metodiku i tehniku pokusa. Ne možemo zamisliti 
pravu ocjenu brojčane vrijednosti nekog svojstva bez varija- 
ciono-statističke obrade. Pokus izveden samo na jednoj biljci 
ili životinji, na jednoj pokusnoj parceli, na jednom gospodar¬ 
stvu i t. d. ne može nam dati pouzdanu sliku o nekom svojstvu. 
Stoga se utječemo zakonu o velikom broju, koji veli: što je veći 
broj varijanata (t. j. pokusnih biljaka, životinja, parcela, go¬ 
spodarstva i t. d.), to će prosjek, odnosno aritmetička sredina 
dati pravilnije podatke za ocjenu istraživanog svojstva. 


b) Svojstva varijanata 

Varijaciono-statistički možemo istražiti: a) kvantitativna i 
b) kvalitativna svojstva stvari i živih bića. Njihova je varijacija 
prema tome a) kvantitativna ili b) kvalitativna. 

Pojedine pokusne životinje, biljke, parcele, gospodarstva 
i t. d., na kojima proučavamo neko svojstvo, zovemo varijan¬ 
tama. Prema tome, da li proučavamo kvantitativna ili kvalita¬ 
tivna svojstva, dijelimo i varijante na kvantitativne i 
kvalitativne. 

Kvantitativne varijante dijelimo na dvije glavne skupine i 
to na: 

a) cijele (diskretne) i b) razredne varijante. 

'/Cijele (diskretne) su varijante one, kod kojih možemo 
istraživano svojstvo izraziti cijelim brojem, kao na pr. broj vlati 
na biljci, broj klasića ili zrna na klasu, broj mahuna na biljci, 
broj sjemena u mahuni i t. d. 

^Razredne varijante su pak one, kojih svojstvo ne mo¬ 
žemo izraziti cijelim brojem. To redovno biva kad istražujemo 
dimenzije, težinu, kemijski sastav i si. pojedinih individua od¬ 
nosno njihovih organa. U razredne varijante ubrajamo na pr. 
duljinu klasa, težinu zrna na klasu, količinu priroda po hektaru, 
količinu mlijeka, koju pojedine krave daju u jednom danu ili 
godini, visinu i težinu životinje, debljinu što je ima nit vune 
i t. d. 





c) Valjanost biometričkih istraživanja 

Iz rezultata, dobivenih biometričkom obradbom, možemo 
samo onda izvesti valjane zaključke, ako je: 

’iS a) uzorak tipičan za materijal, koji treba obraditi, 

b) broj varijanata, koje nam valja istražiti, dovoljno velik 
(zakon o velikom broju), 

c) materijal genetski jednoličan; naprotiv smjesa različnih 
sorata (genotipova), zatim smjesa sjemena istog genotipa ali 
proizvedenog na različitim tlima ili u različitim godinama; kod 
životinja pripadnici različitih pasmina ili iste pasmine ali razli¬ 
čite starosti, težine i dr. manje su podesan materijal za varija- 
ciono-statistička istraživanja. 

Koliko ćemo pokusnih biljaka, životinja, parcelica ili go¬ 
spodarstava uvrstiti u pokus, da dobijemo valjane rezultate, 
zavisi o tome, kako jako varira svojstvo, koje proučavamo: ako 
varira jače, treba uzeti više varijanata, varira li slabije, dostajat 
će i manji broj. 

Broj varijanata potrebnih za pokus zavisi i o tome, u kojoj 
mjeri moraju biti rezultati precizni. O tome će još kasnije biti 
govora. 




B. BIOMETRIČKA ISTRAŽIVANJA 
KVANTITATIVNIH SVOJSTAVA 

U poljoprivredi su od naročite važnosti kvantitativna svoi- 
stva bilja i životinja, stoga treba njihovu varijaciono-statistič- 
kom proučavanju posvetiti naročitu pažnju. Na različnim ćemo 
primjerima vidjeti, koje vrijednosti treba izračunati, da bismo 
neko svojstvo pravilno i točno ocijenili. 


B-l. BIOMETRIKA VARIJACIJE 
a) Varijacioni red 


, --naoc uumaue sorte KuKuruza 

Hrvatice, izmjerit ćemo duljinu od kojih 50 primjeraka. Naj- 
manji je klip dug 16.2 cm, a najdulji 23.7 cm. 

Kad duljine klipova poredamo u skupine ili razrede s razma¬ 
kom ili arealom od, recimo, 1 cm, dobit ćemo slijedeću po¬ 
djelu varijanata nazvanu varijacioni red: 

Tab. II. 
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Srednji razredi imaju najviše varijanata, a na ekstremne 
razrede lijevo i desno otpao je najmanji broj. 

/?) Drugi primjer: na jednoj čistoj liniji (genotipu) Šomot- 
ske pšenice-golice istraživala se prosječna duljina 400 klasova 
s glavnih vlati. (Čista linija je potomstvo jedne homocigotne 
biljke nastale autogamnom konju gacijo m). Najmanji je 
klas dug 6.3 cm, najdulji 14.7 cm. Kad dobivene duljine svih 
400 klasova poredamo u skupine ili razrede s izvjesnim razma¬ 
kom (arealom) na pr. od 0.5 cm, dobijemo slijedeći varijacioni 
red: 


Tab. 111. 


Duljina g_g_5—7— 7 .5— 8 - 8.5— 9— 9.5—10—10.5 — 11 — 13.5-12—12.5—13—13.5—14—14.5—15 


klasa u cmj 


t i 






i 

1 



1 ' 




1 


Broj varijanata u 
pojedinim razre¬ 
dima (frekven¬ 
cija) p 

3 

8 

12 

23 

27 

39 

48 

56 

45 

37 

26 

21 

18 

14 

10 

7 

4 

2 



Suma 0 3 11 23 46 73 112 160 216 261 298 324 345 363 377 387 394 398 400 

frekvencija 
do granica 
razreda 

Kao apsolutnu vrijednost svakog pojedinog razreda uzi¬ 
mamo kod razrednih varijanata sredinu razreda. Na pr. apso¬ 
lutna vrijednost razreda 9.5—10 cm jest v = 9.75 cm. 


b) Grafički prikaz podjele varijanata 


r 





đonom redu sumiramo varijante od razreda do razreda tako, 
da broju varijanata u jednom razredu pribrojimo broj vari¬ 
janata svih prijašnjih razreda, kako je to prikazano na tab. III. 
na str. 18. Ako sumarne frekvencije za svaki razred posebno na- 
nesemo na ordinate razrednih granica, i potom ove točke spo¬ 
jimo crtom, dobit ćemo sumarni poligon. (SI. 4. puna 
linija). 


Predočimo li rezultate gornjih istraživanja tako, da na 
apscisi označimo granice pojedinih razreda, a na ordinati fre¬ 
kvenciju varijanata, pa onda te točke spojimo, dobit ćemo 
stepe nasti poligon (si. 3 puna linija). Ako na apscisu 
nanesemo točke za apsolutnu vrijednost razreda, a na njihove 
ordinate frekvencije, i potom te točke spojimo crtom, dobit 
ćemo varijacioni ili frekvencioni poligon (si. 3, 
iscrtkana linija). Spojimo li spomenute točke krivuljom, nastaje 
varijaciona krivulja. 

Podjelu varijanata u varijacionom redu možemo prikazati 
i pomoću sumarnog poligona. Zato treba da u varija- 


BIOMETRIČKE VRIJEDNOSTI 

Prve i najvažnije statističke vrijednosti, koje moramo izra¬ 
čunati za svako svojstvo jesu apsolutna vrijednost 
svojstva i varijabilitet svojstva. 

I. APSOLUTNA VRIJEDNOST SVOJSTVA 

Apsolutnu vrijednost svojstva može donekle karakterizirati 
vrijednost zvana medijana, a precizno mjerilo je aritme¬ 
tička sredina ili srednja vrijednost. 
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MED. Q 2 

SI. 4. 

1. Medijana — Med 

Za karakteristiku prosječne vrijednosti nekog svojstva upo¬ 
trebljavala se prije vrijednost medijana — Med. Danas se 
ona za tu svrhu rijetko upotrebljava, ali jer ćemo je trebati kod 
izračunavanja nekih drugih važnih varijaciono-statističkih vri¬ 
jednosti, treba da se s njom već sada upoznamo. 

Medijana je u varijacionom redu apsolutna vrijednost 
srednje varijante. Ako na pr. u varijacionom redu za duljinu 
klasa u našem primjeru s pšenicom (str. 18) ima n = 400 vari¬ 
janata, onda je medijana apsolutna vrijednost 200-te varijante 
toga reda. Ovu vrijednost izračunamo ovako: 
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Do razredne granice 9.5 cm ima ukupno 160 varijanata, a 
do razredne granice 10 cm ima ih 216 (vidi str. 18. Tabl, III.). 
Prema tome se mora medijana ili 200-ta varijanta nalaziti u 
razredu između 9.5 i 10 cm. Kako je razredni razmak ili areal 
(a) između 10 — 9.5 cm = 0.5 cm, i u njemu ima 56 varijanata, 

0 5 V 40 

to na 200 — 160 = 40 varijanata otpada : -4g — = 0.357.cm. 

Dakle je vrijednost srednje varijante našeg varijacionog reda 
t. j. Med = q 2 = 9,5 + 0,357 = 9,857 cm. (si. 4.). 

2. Srednja vrijednost : — M 

Mnogo točnije karakterizira apsolutnu vrijednost nekog 
svojstva njegova srednja vrijednost, koju označujemo slovom M. 

Prosječna ili srednja vrijednost je aritmetička sredina, 
koju dobijemo, ako zbroj vrijednosti pojedinih mjerenja (V) 
podijelimo s brojem varijanata (n). 



n 


gdje 2 znači sumu, V apsolutnu duljinu pojedine varijante, a 
n ukupni broj svih varijanata. U varijacionom redu za duljinu 
kukuruznog klipa n = 50. 

Ima i drugi način, kako se izračunava srednja vrijednost: 

U varijacionom redu u tab. III. smjestili smo klasove duljine 
6 do 6.5 cm u prvu skupinu ili razred, klasove s duljinom 6.5 
do 7 cm u drugu skupinu i t. d. Za razredni razmak ili areal 
uzeli smo 0.5 cm. Veličinu tog razmaka određujemo već prema 
varijacionoj širini i karakteru istraživanog svojstva. Varija- 
ciono-statističke vrijednosti su dakako to točnije, što je razredni 
razmak manji. U praksi se radi obično sa 8 do 12 razreda. 

Broj varijanata u jednom varijacionom razredu zovemo 
razrednom frekvencijom. Iz varijacionog razreda izraču¬ 
nat ćemo srednju vrijednost (M) tako, da svaku poje¬ 
dinu apsolutnu vrijednost razreda (v) pomnožimo s odgovara¬ 
jućom razrednom frekvencijom (p), a zbroj (2) ovih produkata 
podijelimo sa zbrojem svih varijanata (n). U tom slučaju- glasi 

formula za srednju vrijednost M = ^ V 

n 

Za izračunavanje srednje vrijednosti upotrebljavamo u 
praksi najčešće slijedeći način: 
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Kao polaznu točku A (aproksimativnu srednju vrijednost) 
uzimamo apsolutnu vrijednost ma kojeg razreda, ali obično vri¬ 
jednost onoga razreda, u kome se otprilike može očekivati 
srednja vrijednost. Za apsolutnu vrijednost razreda uzima se u 
praksi redovno sredina razreda. U našem primjeru za duljinu 
klasa (Tab. III.) uzet ćemo za A apsolutnu vrijednost varijacio- 
nog razreda 10 — 10.5 cm, t. j. sredinu toga razreda, dakle 
10.25 cm. A će vjerojatno biti nešto veća ili nešto manja od 
prave srednje vrijednosti M. Između A i M bit će neka razlika, 
koju ćemo označiti sa b. Prema tome je: M = A ± b. Ako dakle 
označimo na pr. udaljenost varijante C od aproksimativne sred¬ 
nje vrijednosti A sa a, udaljenost A od M sa b i razmak između 
M i C sa a (alfa), onda je a = a 4- b. To se vidi i na slijedećoj 
skici: 


M b 



a 


Ako obje strane jednadžbe a — a + b pomnožimo sa 2p 
t. j. sa sumom (2) pojedinih razrednih frekvencija (p) izlazi: 

2pa = 5 p« t 2pb 

2p . « je zbroj udaljenosti pojedinih varijanata od srednje vri¬ 
jednosti. Ova je suma uvijek jednaka ništici. To se najbolje vidi 
na ovom primjeru: Srednja vrijednost (aritmetička sredina) 
brojeva: 21, 22, 25 i 28 jest 24. Ona je za — 3, — 2, + 1 i + 4 
manja odnosno veća od brojeva, koji čine aritmetičku sredinu. 
Suma ovih razlika je ništa. Iz toga slijedi, da je i u našem slu¬ 
čaju Ep« = 0, tako da je b = budući da je 2p suma 

sviju frekvencija u varijacionom redu ujedno zbroj svih vari¬ 
janata toga reda, dakle n, slijedi da je b — • Budući da 

se u praksi, zbog lakšeg računanja, radi najprije sa a kao jedi¬ 
nicom razrednog areala, mora se b pomnožiti još sa apsolut¬ 
nom vrijednošću tog areala. Dakle konačna formula, po kojoj 
se računa srednja vrijednost, glasi: 

M = A + a b = A -i- a 

n 
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a) Srednja vrijednost razrednih varijanata 

U našem primjeru za duljinu klasa pšenice izabrali smo kao 
polaznu točku A (aproksimativnu srednju vrijednost) apsolutnu 
vrijednost razreda 10 — 10.5 cm. Sredina toga razreda je 
A = 10.25. Ostali razredi imaju razredni razmak (areal = a) 
za 1 X 0.5, 2 X 0,5, 3 X 0.5 cm i t. d. pozitivno ili negativno 
od A. 

Sve brojke potrebne za izračunavanje srednje vrijednosti 
možemo poredati tako, kako je prikazano u tab. IV. 


Tab. IV. 



a = 0,5 
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Srednja je vrijednost: 

M = A + ab = 10.25 + 0.5(—0.4775) = 10.01125 cm ili 
okruglo 10.01 cm, t. j. prosječna je duljina klasa spomenute 
čiste linije Šomotske pšenice 10.01 cm. 

Kraći način za izračunavanje 2 (pa) je ovaj: (Tab. V.) 
(Tabela je proširena rubrikama, koje služe za izračunavanje 
i drugih varijaciono-statističkih vrijednosti). Za izračunavanje 
srednje vrijednosti odnosno 2 (pa) služe nam horizontalne ru¬ 
brike od 4 do 8). 


Tab. V. 



pozitivni produkti negativni produkti 

(rubrika 7 i 4) (rubrika 8 i 4) 


(P a) 

1-8 ~8 
2-9 18 

Ukupno 26 


2 (pa) = — 217 + 26 = — 191 


(pa) 


19 -1 

— 19 

22 -2 

— 44 

18 . - 3 

— 54 

9 • - 4 

— 36 

9 •• — 5 

— 45 

2—6 

- 12 

1—7 

— 7 

Ukupno 

— 217 





Kontrolu vršimo na različne načine: a) da ponovimo račun, 
b) da račun izvedemo s koje druge polazne točke za A (apro¬ 
ksimativne srednje vrijednosti) i c) Charlierovom me¬ 
todom. 


Ova je metoda radi svoje jednostavnosti naročito prikladna. 
Za novu polaznu točku Ai izabere se ona, koja je za jedan areal 

Sna. 

manja od pređašnje A. Polaznoj točci A odgovara b = - ^ , 

n 

a polaznoj točci Ai odgovara b x = 2a kontrolu je 

dovoljno da izračunamo 


2[p(a + l)] = 2pa -f 2p = Spa + n 


Ako je dakle nova polazna točka (aproksimativna srednja 
vrijednost) duljine klasa: A x — 9.75 cm, onda je: 


Tab. VI. 



Iz gornje tabele proizlazi: 
6-5 30 _ 3.1 _ 3 

6 - 6 36 2-2 — 4 

7 -7 49 _ 1.3 _ 3 

7- 8 56 _ 2.4 - 8 

4 . 9 36 - 

2-10 20 —18 

227 


Dakle: 2p(a+l) = 2pa + n = 227 — 18 = 209 = —191 + 400 
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Osim metode množenja, kojom se izračunava Spa postoji 
još i metoda zbrajanja, koja se jako malo upotrebljava, pa je 
zato ne ćemo ni primijeniti. 


b) Srednja vrijednost cijelih (diskretnih) varijanata 

Cijele ili diskretne varijante su, kako znamo, one 
kod kojih se vrijednost svojstva može izraziti cijelim brojem. 
Tako na pr. broj klasića na klasu spada u skupinu diskretnih 
varijanata. 

U spomenutoj već čistoj liniji Šomotske pšenice-golice istra¬ 
živano je osim duljine klasa i broj klasića na klasovima glavnih 
vlati. Najmanji broj klasića je u toj liniji 11, a najveći 25 (vari- 
jaciona je širina 25—11 = 14 klasića). Ako za razredni razmak 
(areal = a) uzmemo, da je a = 1 (klasić), onda je varijacioni 
red slijedeći: 


Tab. VII. 


Broj klasića 

11 

12 

13 

14 

F 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

Ukupno 

vari¬ 

janata 

Razredna 

frekvencija 

3 

10 

16 

28 

41 

59 

57 

49 

38 

_ 

29 

25 

19 

15 

7 

4 

400 


Srednju vrijednost broja klasića izračunamo na ovaj način: 

Za polaznu točku (aproksimativnu srednju vrijednost A) 
uzmemo na pr. razred, u kojemu 49 klasova ima po 18 klasića, 
dakle A — 18. (Kao A možemo dakako upotrebiti apsolutnu 
vrijednost ma kojega razreda.) 

Prema tome bi onda tab. VIII. za izračunavanje srednje vri¬ 
jednosti po skraćenoj metodi bila ovakova: (tabela VIII. je pro¬ 
širena za izračunavanje i drugih različnih statističkih vrijed¬ 
nosti). 


Tab. VIII. 


Rubrika 

i 


a 3 

i 

8 

27 

64 

125 

216 

343 

2 


a 

i 

2 

3 

4 

5 

6 

7 

3 

Frekvencija 

(P) 

+ 

38 

29 

25 

19 

15 

7 

4 

4 

— 

57 

59 

41 

28 

16 

10 

3 

5 

Diferencija 

(P) 

+ 

— 

— 

— 

— 

— 

— 

1 

6 

— 

19 

30 

16 

9 

1 

3 

— 

7 

Suma 

(P) 


95 

88 

66 

47 

31 

17 

7 

8 


a* 

1 

4 

9 

16 

25 

36 

49 

9 


a 4 

1 

16 

81 

256 

625 

1296 

2401 


pozitivni produkti 

negativni produkti 

(rubrike 5 i 2) 

(rubrike 6 i 2) 

(pa) 

(pa) 

1.7 = 7 

19.—1 = — 19 


30 . —2 = — 60 


16 . —3 = — 48 


9 . —4 = — 36 

= — 186 + 7 = — 179 

1 . — 5 = — 5 

“ 179 0447*5 

3 . — 6 = — 18 

= To--°- 4475 

Ukupno —186 


M = A + ab = 18 + l (—0.4475) = 17,5525 ili okruglo 17,6. 


Dakle srednja je vrijednost za broj klasića M = 17.6 kla¬ 
sića na klasu. 
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II. RAČUNANJE VARIJABILITETA 

Za točnu ocjenu nekog svojstva treba da znamo njegovu 
srednju vrijednost i njegov varijabilitet. Dvije sorte pšenice ili 
dvije pasmine goveda, koje imaju na pr. istu srednju vrijednost 
za neko svojstvo, nemaju jednaku gospodarsku vrijednost, ako 
je varijabilitet toga svojstva kod jedne sorte ili pasmine veći, 
a kod druge manji. Što je varijabilitet kod neke sorte gospo¬ 
darske biljke ili neke pasmine domaćih životinja i si. veći, to 
je ta sorta ili pasmina slabija. Za mjerenje varijabiliteta služe 
ove vrijednosti: 


1. Varijaciona širina — Vš 

Varijabilitet nekog svojstva možemo približno mjeriti po¬ 
moću varijacione širine. JPfiđ varijacionom ši¬ 
rinom razumijevamo razliku između apsolutne vrijednosti 
svojstva najveće i .najmanje varijante. Na pr. između 50 klipova 
gore spomenutog kukuruza Hrvatice, najmanji (min.) je klip 
bio dug 16.2 cm, a najdulji (max.) 23.7 cm. Varijaciona je širina 
prema tome Vš = max. — min. = 23.7 — 16.2 = 7.5 cm. Od 
400 klasova Šomodske pšenice-golice najmanji je bio dug 6.3 
cm, a najveći 14.7 cm, dakle je varijaciona širina za duljinu 
klasa Vš = 14.7 — 6.3 = 8.4 cm, a budući da je najmanji broj 
klasića između 400 klasova bio 11, a najveći 25, to je Vš = 
= 25 — 11 = 14 klasića varijaciona širina. (Str. 26, tab. VII.). 


2. Kvartil — Q 

Kvartil je vrijednost, koja rezultira iz diferencije između 
apsolutne vrijednosti (q 3 ) zadnje varijante u trećoj četvrtini 
varijacionog reda i (qj) apsolutne vrijednosti zadnje varijante 
prve četvrtine varijacionog reda, podijeljenom sa 2, dakle 

Q= ~ ~2 ^ 

U našem primjeru s duljinom klasa pšenice q 3 znači apso¬ 
lutnu vrijednost 300-te varijante, a q x apsolutnu vrijednost 
100-te varijante. 
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Iz tabele III. na str. 18 i si. 4 slijedi, da je 
q t = 8,5 4- 0,346 = 8,846 cm 

q 3 = 11 + 0,038 = 11,038 cm, tako da je Q = ’ — = 

= 1,096 cm. Kvartil je dakle apsolutno mje rilo za varijabilitet 
svojstva, ali nije točno, nego tek aproksimativno, jer obuhvaća 
svega 50% varijanata i to onih iz sredine varijacionog reda. 

U normalnoj je krivulji Q = 0.6745 a. U toj formuli je a 
standardna devijacija. 

3. Standardna devijacija — 0 

Varijabilitet nekog svojstva najbolje se mjeri pomoću 
standardne devijacije (a). Standardna devijaci ja se 
temelji na metodi najmanjih kvadrata. Principe, po kojima ovom 
metodom izračunavamo varijabilnost svojstva, upoznat ćemo 
najbolje na slijedećem primjeru: 

Dva lovca X i Y natječu se, koji će bolje pogoditi cilj. Svaki 
gađa dvaput. Lovac X je oba puta pogodio u 2. krug cilja, 
a lovac Y pogodio je najprije u 1., a onda u 3. krug. Pita se, 
koji je od ove dvojice bolji strijelac? (SI. 5.). 



Ako zbrojimo rezultate gađanja lovca X, dobit ćemo 
24-2 = 4, a rezultate lovca Y: 14-3=4. To bi značilo, da su obo¬ 
jica jednako dobri strijelci. Ipak je X trebao manju ukupnu po¬ 
vršinu da pogodi cilj nego li Y. 
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Površina kruga je r%: Lovac X pogodio je površinu od 
2 ji( 2) 2 , a lovac Y od Jt(l) 2 +jt(3) 2 . Prema tome su ove dvije 
površine u relaciji 8 :10. 

Još ćemo bolje ocijeniti streljačku vještinu obaju lovaca, 
ako se poslužimo brojkom, koju smo dobili iz drugog korijena 
sume kvadrata, sastavljenih od udaljenosti pogodaka od cilja 
prvog i drugog gađanja lovca X, odnosno lovca Y. 

Na taj način dobivamo slijedeće ocjene: 

lovac X .... y2?~+2? = 1/8“ = 2.83 

„ Y .... l/T* + 3* = 1/10 = 3.16 
a to znači, da je varijabilitet u gađanju kod lovca X manji (2.83) 
nego kod lovca Y (3.16), te je prema tome X bolji strijelac 
nego Y. 

Dakle varijabilitet pojedinih svojstava možemo najbolje 
izraziti pomoću drugog korijena iz sume kvadratnih udaljenosti 
(a ) 2 pojedinih varijanata od srednje vrijednosti, podijeljene s 
brojem varijanata (n). Prema tome je opća formula za odre¬ 
đivanje varijabiliteta svojstva ili standardne devijacije a ova: 



Iz jednadžbe: a=a + b upotrebljene već kod računanja 
srednje vrijednosti (vidi str. 22) izvodimo formulu za standardnu 
devijaciju, najpraktičniju za rad u ovom obliku: 

Ako jednadžbu a=a+b kvadriramo, dobit ćemo a 2 = a 2 + 
+ 2«b + b 2 , a pomnožimo li ovu jednadžbu sa sumom varijanata 
2p, dobit ćemo 2pa 2 = 2pcr + 22p«b + 2p . b 2 . 

Suma razmaka između pojedinih varijanata i srednje vri¬ 
jednosti, kako je već spomenuto, ravna je ništici. Zato je: 

2pa = 0 i 2p2ab = 0. 

Preostaje 2pa 2 = 2pa 2 + 2pb 2 
U iznosu 2pb 2 je 2p=n. 

Iz toga slijedi: 2pa 2 = 2pa 2 —nfc 2 ili 

0 2 = 2p<* 2 _ Upa 2 _ , 2 

n n 


Tamo gdje je areal a—1 računamo standardnu devijaciju po 



formuli: a = 
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Ako je areal (a) veći ili manji od 1, moramo o pomnožiti 
još s apsolutnom vrijednošću areala, dakle sa a. 

Opća formula za standardnu devijaciju, koja vrijedi za sve 


vrijednosti areala onda je slijedeća: o = a f— E-?-b 2 . 

Primjenu te formule prikazat ćemo'na našem već"spome¬ 


nutom uzorku pšenice. 


a) Standardna devijacija razrednih varijanata 

U varijacionom redu za duljinu klasa pšenice tab. III., str. 
18 pošli smo od aproksimativne srednje vrijednosti A= 10.25 
đm. Ovom vrijednošću služit ćemo se i kod računanja standardne 
devijacije. 

Vrijednosti a, b i n poznate su nam još od izračunavanja 
srednje vrijednosti. One su: a=0.5, b=—0.4774, n=400; 2pa 2 = 
=4711 izračunat ćemo iz tab. IV. ili V. na str. 23 i 24. 

Po metodi množenja tab. V. zbrajamo frekvencije (p) ra¬ 
zreda, koji su u pozitivnom i negativnom smjeru od točke A 
udaljeni za: a, 2a, 3a i t. d. U našem je slučaju na pr. frekven¬ 
cija (p) u razredu za a udaljenom od A u pozitivnom smjeru 37, 
a u negativnom 56, dakle 37+56=93. Zatim zbrajamo redom 
frekvencije u razredima udaljenim za 2a, 3a, 4a i t. d. od A. 
Svaku od ovih suma pomnožit ćemo onda s kvadratom udalje¬ 
nosti razreda od A. 

Na našem primjeru dobivamo tako slijedeće produkte i 
sume: 




Za duljinu klasa je standardna devijacija: 


a = a 



1/4711 ~~ 

400 ( 0 


• 4775) 2 = ± 1-699 


b) Standardna devijacija cijelih (diskretnih) varijanata 


Varijacioni red za cijele varijante, i to broi 
klasu, prikazan je u tab. VII. na str. 26. 


klasića na 


Za aproksimativnu srednju vrijednost (A) uzeli smo razred 
u kojemu ima 49 klasova, svaki sa 18 klasića, tako da je A = 18.’ 
Standardnu devijaciju izračunavamo i ovdje po formuli: 



S vrijednostima a, b i n smo se upoznali već kod izračuna¬ 
vanja srednje vrijednosti — M. (Vidi str. 26 i 27) i to: 


a— 1, b=0.4475 i n=400 

Spa” izračunavamo iz tab. VIII. ovako: (rubrike 7 i 8). 



Dakle broj klasića na klasu ima standardnu devijaciju: 



. . standardna devijacija izražava apsolutnu varijabilnost 
is raženog svojstva. Kod normalne krivulje ona označuje one 
točke krivulje projicirane na apscisu, kod kojih centralna kon- 
kavnost prelazi u Iateralnu konveksnost (si. 2.). 
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Računske operacije kontroliramo i ovdje najbolje Char- 
lierovom metodom 'na taj način, da za novu polaznu točku 
odaberemo A 1; koja je za jedan areal manja od A. 

Kad radimo s polaznom točkom A, onda je 



a kod operacija sa Ai je a = ± a 


Sp(a+ l) 2 


Za kontrolu nam je dovoljan izraz: 


2p(a + l) 2 =2pa 2 f 22pa + 2p= 2pa 2 4-22pa + n 


Ako je broj varijanata malen dobit ćemo točnije rezultate, 
poslužimo li se slijedećom formulom za izračunavanje standardne 
devijacije: 



Na veličinu standardne devijacije utječe i veličina razrednog 
razmaka (areala = a). Utvrđeno je, da je standardna devijacija 
nešto veća, ako su areali veliki, a manja, ako su maleni. To za¬ 
visi o podjeli apsolutnih vrijednosti pojedinih varijanata u ra¬ 
zredima. U praksi se, zbog jednostavnijeg postupka, uzima 
apsolutna vrijednost pojedinih razreda, kao da je u sredini 
razreda, iako to nije uvijek tako. Upravo poradi toga izradio 
je Američanin Sheppard za pravilniju standardnu devijaciju 
formulu, koja glasi: 


" /"Spa 2 _ 

“=± a |/-„—»■ 


Y2 odnosno a = ± a 



Od kolike je važnosti standardna devijacija za mjerenje 
varijabilnosti svojstva, vidi se iz slijedećeg primjera iz prakse: 

Pomoću sortnog pokusa ustanovili smo, da je na pr. pro¬ 
sječni prirod zrna sorte pšenice A, a isto tako i sorte B = 20.5q 
po ha. Standardna devijacija sorte A je o=±2.4, a sorte B pak 
0= ±0.7 q po ha. Iz toga slijedi, da prirod sorte A nije tako kon¬ 
stantan kao prirod sorte B, a to stoga, što sorta A ima veću 
varijabilnost u prirodu nego sorta B. Sorta B je prema tome 
s gospodarskog gledišta bolja od sorte A. 

Standardna devijacija služi dakako kao apsolutno mjerilo 
varijabiliteta. 


3 Tavčar: Biometrika u poljoprivredi 
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4. Varijacioni koeficijent — v 


Standardna devijacija je apsolutno mjerilo za varijabilitet 
nekog svojstva. 

Da usporedimo, kako jako variraju različna svojstva, mo¬ 
ramo upotrebiti neko relativno mjerilo. Takvo je mjerilo 
varijacioni koeficijent — v. Varijaci oni koeficijent 
izražava, standardnu devijaciju u .postocima srednje vrijednosti, 


dakle: v = 


100 o 
M 


v = 


U našem je primjeru za duljinu klasa varijacioni koeficijent: 

100 X 1,699 


10,01 


16,97, a za broj klasića na klasu 


v = 


100 X 2,934 
17,6 


= 16,67. 


Duljina klasova varira dakle nešto jače nego broj klasića. 


5. K v a r t i! n i koeficijent — Q k 

Tamo, gdje želimo usporediti varijabilitet različnih svoj¬ 
stava, na pr. duljinu klasa s brojem klasića na klasu pomoću 
kvartila — Q (vidi str. 28), valja nam upotrebiti neko rela¬ 
tivno mjerilo. 

U tu nam svrhu služi k v a r t i 1 n i koeficijent - Qk. 
Kvartilni koeficij enat izražava kvartil u postocima 

100 Q 


srednje vrijednosti. Dakle Qn = 


M 


Za duljinu klasa u našem je slučaju: (str. 29) 

Q _ 100 X 1,096 

io.oi lu,y4 - 


III. SREDNJA I VJEROJATNA POGREŠKA — m i E. P. 

Pretpostavimo li, da su srednja vrijednost, standardna de¬ 
vijacija i varijacioni koeficijent izračunati iz prosječnog uzorka 
i uza sve već spomenute opreze, ipak nas zanima, do koje mjere 
oni karakteriziraju istraživano svojstvo i do koje je granice 
moguće njihovom pomoću izvesti valjane zaključke za vrijed¬ 
nost toga svojstva. 


Sigurno je, da što više varira materijal koji istražujemo, to 
je manja vjerojatnost, da će uzorak iz toga materijala — na¬ 
ročito ako je manja količina (mali broj varijanata) — biti pro¬ 
sječan. Velika ili bar veća količina uzoraka (veliki broj vari¬ 
janata) mogla bi svakako bolje reprezentirati materijal, iz ko¬ 
jega je bila odvojena. 

Pomoću sta ndar d ne devijacije možemo ustanoviti varija¬ 
bilnost nekog svojstva. Na tab. IX. vidimo, koji se postotak va¬ 
rijanata nalazi u granicama između sredn je vrijednosti M i 
multiplog broja od a. 


Tab. IX. 


Razmak od 
srednje vrijednosti 

M ± X <3 

°/ 0 varijanata 
u granicama 

M ± x o 

°/o varijanata 
izvan granica 

M + xo 

M ± 0-5 a 

38-3 

61-7 

M + TO o 

68-3 

31-7 

M ± 1-5 o 

86-6 

13-4 

M ± 2'0 o 

95-5 

4"5 

M ± 25 a 

98-8 

1-2 

M ± 30 « 

99-7 

0"3 


Prema tome se u granicama M ± 3o nalazi 99.7% svih va¬ 
rijanata teoretskog varijacionog reda, a to znači, da možemo 
s vjeroja tnošću od 99.7% tvrditi, da je srednja vrijednost (M) 
najviše 3o udaljena od jedne slobodno izabrane varijante. 

Standardna devijacija ima dakle dvostruku važnost i to: 
L. kao mjerilo varijabilnosti i 2, kao indic ij. na osnovu koga 
možemo zaključiti iz koje mu drago varijante u varijacionom 
redu, kolika je srednja vrijednost toga reda. 

Da ustanovimo, do koje mjere neka biometrička vrijednost 
karakterizira neko svojstvo organizma ili mrtve tvari, valja 
nam izračunati njezinu srednju odnosno vjerojatnu pogrešku. 

Srednja i vjerojatna pogreška su mjerila, koja n am kazuju, 
možemo li dotičnu varijaciono-statističku vrijednost upotrijebiti 
za izvođenje valjanih zaključaka il i ne možemo. 

* 


i 







1. Srednja pogreška — m 


a) Srednja pogreška srednje vrijednosti — rn M 


Ako je oj pogreška jedne slobodno izabrane varijante, onda 
je srednja pogreška srednje vrijednosti od n varijanata: mw = 


= ■ ° t. j. standardna devijacija dividirana s drugim korije- 

V n 

nom od broja svih varijanata. 

U našem je primjeru za prosječnu duljinu klasa kod pšenice 
srednja vrijednost M = 10.01 cm (str. 24), a standardna devi¬ 
jacija a = ± 1,699 (str. 32). 

Prema tome je srednja pogreška srednje vrijednosti: 



1,699 

1/400 


= 0,085 


Da prosudimo valjanost srednje vrijednosti — M obzirom 
na njezinu srednju pogrešku, treba pisati: M±m. 

Prema tome je u našem slučaju M±m = 10.01±0.085. 

Srednja vrijednost za broj klasića je u našem primjeru 
(str. 27) M=17.6 i standardna devijacija o=±2,934, tako da je 
srednja pogreška srednje vrijednosti: 

2 Q34 

triM = T ~ — - = 0.1467 odnosno 0.147 
|/400 

Dakle je M±m = 17.6±0.147. 

Dokazano je da za stvaranje valjanih zaključaka srednja 
vrijednost M mora biti veća od trostruke srednje pogreške—3m. 


a) Granice bio m e t r i č k i opravdane srednje 
vrijednosti 


Srednja pogreška m označuje nam granice, u kojima se 
kreće srednja vrijednost M izračunata od n varijanata. 

Ako je broj varijanata n velik,""orida su'graniče točnosti, 
odnosno dopuštene pogreške srednje vrijednosti u pozitivnom 
i negativnom smjeru od »M« jednake trostrukoj srednjoj po¬ 
grešci (3m) srednje vrijednosti. Ako je broj varijanata malen 
netočnost se povećava. 

U donjoj su tabeli (tab. X.) prikazane granice točnosti od¬ 
nosno dopuštene pogreške u umnošku od m i to u varijacionom 


redu sa n—1 do 10 varijanata (n—1 uzima se radi toga, jer 
srednja vrijednost, izračunana od pojedinačnih mjerenja, ima 
n-^1 neovisnih varijanata, a samo je jedna varijanta ovisna). 


Tab. X. 


Broj neovisnih 
varijanata 
n —1 

Umnožak od m 
t 

Broj neovisnih 
varijanata 
n — 1 

Umnožak od m 
t 

1 

235.8 

6 

4.9 

2 

19.21 

7 

4.53 

3 

9.22 

8 

4.27 

4 

6,62 

9 

4.09 

5 

5,51 

10 

3.96 


U Kollerovu diagramu (si. 6) možemo naći granice točnosti 
za 10 do co broj varijanata. 

t 
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U našem primjeru za duljinu klasa pšenice, gdje je n—1 = 
= 400—1=399 slobodnih varijanata, srednja pogreška srednje 
vrijednosti (str. 36) m m =0.085 je t=3,02 (vidi Kollerov dia- 
gram si. 6). Prema tome su granice točnosti za duljinu klasa, 
gdje je srednja vrijednost M= 10,01 (str. 24) između 
M—tmM= 10,01—3,02X0,085=9,753 i 
M+tm M = 10,01 +3,02X0,085=10,267 

V 

b) Srednja pogreška standardne - devijacije 

— m 0 


Ovu vrijednost izračunavamo oo formuli: m 0 = , ?— • 

y 2 n 

U gore spomenutom našem primjeru s duljinom klasa je (vidi 

1 699 

str. 32): m 0 = 0.06 odnosno o±m 3 - 1,699±0,06, a 

za broj klasića na klasu je srednja pogreška standardne devija-' 

2 934 

čije: m 0 = p 7 ======= 0404. Prema tome je c±m„=2,934± 

±0.104 (vidi str. 32). 


Na temelju standardne devijacije možemo izvoditi pravilne 
zaključke samo onda, ako je ona veća od svoje trostruke srednje 
pogreške. 


c) Srednja pogreška varijacionog 
koeficijenta — m v 

Ako je varijacioni koeficijent v manji od 10, računamo nje- 

v \ 

govu srednju pogrešku po formuli: m T = __ , a ako je veći 

|/2 n 

od 10, po formuli m, = ±= ](1 +.2 (-i.)'. 

U našem primjeru za duljinu klasa (str. 34) varijacioni 
koeficijent v = 16.97, a u primjeru za broj klasića v = 16.67. 

Budući da je jedan i drugi varijacioni koeficijent veći od 
10, izračunavamo njihovu srednju pogrešku m v po drugoj for¬ 
muli. 
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Prema tome je u našem primjeru srednja pogreška vajdja- 
cionog koeficijenta za duljinu klasa: 

m v = t - 6,97 _ 1/7+2 = 0.6169 , 

1/2X400 1 \ 100 / 


a za broj klasića: m v = f1 72 (^)"= 0.6054. 

U jednom i drugom primjeru varijacioni koeficijent veći je 
od svoje trostruke srednje pogreške. Prema tome nam ovi 
varijacioni koeficijenti mogu poslužiti za pravovaljane za¬ 
ključke. 


2. Vjerojatna pogreška — E. P. (Error probabilis) 

Već je spomenuto, da k v a r t i 1 — Q označuje gra nice, u 
kojima se nalazi srednjih 50% svih varijanata izvjesnog"virT- 
jacionog reda. U takvom varijacionom redu, gdje su varijante 
razdijeljene točno po binominalnoj formuli, kvartil se i stan¬ 
dardna devijacija nalaze u jednom posve stalnom omjeru i to: 
Q : o=0.6745 . 

Dok evropski stručnjaci procjenjuju varijaciono-statističku 
opravdanost rezultata pomoću srednje pogreške m, engleski se 
i američki stručnjaci služe vjerojatnom pogreškom 
(E. P. = error probabilis). Nju dobivamo tako, da srednju po¬ 
grešku (m) pomnožimo s faktorom 0.6745. 

Prema tome glase formule za: 

a) vjerojatnu pogrešku srednje vrijednosti: 

EP m = ^ 0,6745, 

V n 

b) vjerojatnu pogrešku standardne devijacije: 

EP 0 = =-L= 0,6745 
y 2n 

c) vjerojatnu pogrešku varijacionog koeficijenta: 


EP v = =+= 0,6745, 

H? _ 


ako je v < 10 i 


-=|+r F 1 + 2 ("mo )"°- e745 
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Na osnovu M, a i v možemo-donijeti valjane zaključke samo 
onda, ako su te vrijednosti triput veće od njihovih vjerojatnih 
pogrešaka. 

Za binominalnu krivulju s 10.000 varijanata vjerojatnost, 
da prava srednja vrijednost (M) leži između multiplih brojeva 
od vjerojatne pogreške (E), prikazana je u tab. XI. 


Tati. XI. 


Granice za vjerojatnu 
pogrešku (EP) od M 

M ± E. P. 

Vjerojatnost 

M ± E. P. 

1 : 1 

M ± 2 E. P. 

4,64 : 1 

M ± 3 E. P. 

22,36 : 1 

M + 4 E. P. 

141,87 : 1 


Ako je dakle srednja vrijednost (M) tri puta veća od njene 
vjerojatne pogreške (EP), onda je vjerojatnost 1 :22,36 da je M 
varijaciono-statistički opravdana. 

Srednju i vjerojatnu pogrešku izražavamo vrlo često i u 
postocima. 

Izračunavanje srednje odnosno vjerojatne pogreške ne¬ 
ophodno je potrebno, da ocijenimo, jesu li rezultati, dobiveni 
različnim pokusima na pr. s prirodom različnih sorata poljo¬ 
privrednih biljaka, s gnojivima ili količinom nekog gnojiva, s 
dobom sjetve, različnim razmacima biljaka u polju, različnim 
gospodarskim životinjama i t. d. upotrebljivi za valjane za¬ 
ključke. 

Iz rezultata različnih pokusa, objelodanjenih u stručnim 
časopisima, možemo dobiti tek onda pravu sliku o njihovoj va¬ 
ljanosti, znamo li srednju, odnosno vjerojatnu pogrešku istra¬ 
ženih svojstava. 

IV. POREĐENJE VARIJACIONIH REDOVA S TEORETSKOM 
RASPODJELOM (DISTRIBUCIJOM) 

Za točniju ocjenu izvjesnog svojstva treba katkada da po¬ 
djelu varijanata u varijacionom redu usporedimo s teoretskom 
t. j. binominalnom raspodjelom varijanata. 


To možemo učiniti na dva načina i to: 

1. poređenjem brojaka i 2. grafičkim prikazom. 

1. Poređenje brojaka. Matematičari su izradili tabele, po¬ 
moću kojih možemo pronaći, koliko se varijanata od ukupnog 
broja 10.000 nalazi kod idealne ili binominalne podjele između 
idealne sredine binominalnog reda i ma kojega razmaka u po¬ 
zitivnom ili negativnom smjeru od te sredine. 

Udaljenost pojedinih varijanata od srednje vrijednosti izra¬ 
žena je u tom slučaju u jedinicama standardne devijacije. 

Ako sa D označimo apsolutni razmak izvjesne varijante od 
srednje vrijednosti, onda je D : o relativno mjerilo za razmake 
pojedinih varijanata (a = standardna devijacija). 

Na priloženoj tabeli (vidi tab. XII.) izračunate su vrijednosti 
za D : a. Upotrebu ove tabele prikazat ćemo na našem primjeru 
pšenice za duljinu klasa. 

Za duljinu klasa smo već prije izračunali srednju vrijednost: 
M = 10.01 cm (str. 24), isto tako i standardnu devijaciju 
a = + 1,699 (str. 32). Nadalje treba izračunati udaljenost (D) 
pojedinih razrednih granica od srednje vrijednosti M. Do ra¬ 
zredne granice 6 cm (str. 18, tab. III.) razmak je od siednje 
vrijednosti V — M = 6 — 10.01 = — 4.01; razredna granica 
6.5 udaljena je od srednje vrijednosti za 6.5 — 10.01 = —3.51 
i t. d. Prema tome je za razrednu granicu 6 omjer: D:o = 

=_4,01 : 1,699 = — 2,3602. Za razrednu granicu 6.5 je omjer 

D : a = — 3,51 : 1,699 = — 2,0659 i t. d. 

U tab. XII. potražit ćemo za D : a, koliko varijanata leži 
između srednje vrijednosti M (razmak je 0) i pojedinih razred¬ 
nih granica, kad bismo imali 10.000 varijanata. Ove su brojke 
unesene u rubriku 3 na našoj tabeli br. XIII. Na primjer, za ra¬ 
zrednu granicu 8 cm je D : a = 1,183. 

U tab. XII. naći ćemo izvan granica razreda 3749 varijanata 
za D : o = 1,15. Interpolacijom izračuna se onda za D : a= 1,183, 
a to u našem slučaju čini 3815 varijanata. 

U koloni 5 tab. XIII. izračunat je teoretski broj varijanata 
za svaki pojedini razred. Ove brojke dobijemo tako, da u 4. 
koloni suptrahiramo susjedne brojke u vertikalncmi smjeru,, 
dakle za razred 7.5—8 cm broj je varijanata: 3815-$400=485. 

Prve susjedne brojke do srednje vrijednosti (M), gdje je 
razmak 0, treba dakako zbrojiti. U 4. i 5. koloni su brojke pre¬ 
računate na 10.000 varijanata. 
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Tab. XII 


D: o za 

10.000 

■ 2 ° 

D: o za 

10.000 

cv 

P 
— o 

D:c 

za 10.000 

varijanata 

cu NJ 

a: N 

varijanata 

CO g 

Om N 

varijanata 

0,00.... 

.. 0 

40 

1,50.... 

..4332 


3,00. 

.4987 

0,05. 

0,10. 

.. 199 

.. 398 

40 

AC\ 

1.55.. .. 

1.60.. .. 

..4394 

.. 4452 

11 

3,05. 

3,10. 

.4989 

0,15.... 

.. 596 

39 

1,65.... 

..4505 

1 1 

10 

3,15. 

.4992 

0,20.... 

.. 793 

39 

1,70.... 

..4554 

9 

3,20. 


0,25. 

.. 987 

38 

1,75.... 

.. 4599 

Q 

3,25. 


0,30. 

..1179 

38 

1,80. 

.. 4641 

7 

3,30. 

.4995 

0,35. 

..1368 

37 

1,85. 

.. 4678 

7 

3,35. 


0,40. 

..1554 

37 

1,90. 

.. 4713 


3,40. 


0,45. 

.. 1736 

36 

1,95. 

..4744 


3,45. 


0,50. 

,.1915 

35 

2,00. 

.. 4773 

O 

5 

3,50. 

.4998 

0,55. 

.. 2088 

34 

2,05. 

.. 4798 


3,55. 


0,60. 

..2258 

33 

2,10. 

.. 4821 

D 

3,60. 


0,65. 

.. 2422 

32 

2,15. 

..4842 

4 

3,65. 


0,70. 

..2580 

31 

2,20. 

..4861 

4 

3,70. 

.4999 

0,75. 

..2734 

30 

2,25. 

..4878 


3,75. 

.4999 

0,80. 

..2881 

28 

2,30. 

..4893 

o 

3,80. 

.4999 

0,85. 

. 3023 

27 

2,35. 

..4906 

o 

3,85. 

.4999 

0,90. 

..3159 

26 

2,40. 

.. 4918 

2 

3,90. 

.5000 

0,95. 

.. 3289 

25 

2,45. 

.. 4929 

o 

CO . 

.5000 

1,00. 

..3413 

24 

2.50. 

.. 4938 

o 



1,05. 

. 3531 

22 

2,55. 

.. 4946 

1 



1,10. 

. 3643 

21 

2,60. 

.. 4953 

1 



1,15. 

. 3749 

20 

2,65. 

..4960 

1 



1,20. 

.3849 

19 

2,70. 

. 4965 

\ 



1,25. 

. 3944 

18 

2,75. 

. 4970 

l 



1,30. 

.4032 

17 

2,80. 

. 4974 

1 



1,35. 

.4115 

15 

2,85. 

.4978 

1 



1,40. 

.4192 

14 

2,90. 

..4981, 




1,45. 

.4265 

13 

2,95. 

.4984 





U 6. koloni je teoretska, a u 7. koloni pokusom dobivena 
podjela varijanata u varijacionom redu za n = 400 varijanata. 

Ako u našem slučaju usporedimo teoretsku podjelu vari¬ 
janata s podjelom dobivenom iz pokusa vidimo, da između njih 
nema većih razlika. 

2. Grafički prikaz. Do koje se mjere podjela varijanata u 
varijacionom redu eksperimenta slaže s teoretskom, binominal- 
nom podjelom, možemo vidjeti i tako, da tu podjelu prikažemo 
grafički. To se dade izvesti na nekoliko načina i to u našem 
primjeru za duljinu pšeničnog klasa: 

a) Na apscisu nanesemo granice pojedinih razreda a na 
ordinate onda pripadajuće frekvencije varijanata i to pokusom 
dobivene i teoretske (tab. XIII., kolona 7 i 6). Spojimo li linijom 
najprije jedne, a onda druge tako dobivene točke, vidjet ćemo, 
do koje se mjere one podudaraju (si. 7). 



SI. 7. 


b) Vrijednosti brojaka dobivenih zbrajanjem frekvencija po¬ 
jedinih razreda — (vidi str. 18, tab. III. i kolona 9 u tab. XIII.) 
— nanest ćemo na ordinate (razmak je na apscisi označen jedi¬ 
nicom za razredni areal). Spojimo li točke na ordinatama crtom, 
dobit ćemo sumarnu krivulju (si. 4) za duljinu klasa. 
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1155 46,20 56 54,00 


+ 0,6 40 


4' 



400,00 


Na isti se način konstruira iz teoretskih brojaka kolone 10, 
tab. XIII. (si. 4) i sumarna krivulja za duljinu klasa. 

c) x Uspoređivanje binominalne sumarne kri- 
v u 1 j e s eksperimentalnom sumarnom krivu¬ 
ljom pomoću Hazenova dijagrama 

Dokazano je, da vrijednosti integrala frekvencije za sumarnu 
binominalnu krivulju, izražene u postocima i prikazane gra¬ 
fički, daju pravac. Na osnovu toga izrađen je Hazenov dijagram, 
koji vidimo na si. 8. 

Da ustanovimo, do koje se mjere eksperimentalni podaci 
poklapaju s teoretskim, komparirat ćemo u Hazenovu papiru 
liniju konstruiranu od eksperimentalnih podataka s pravcem 
dobivenim od teoretskih podataka. 

Teoretski pra v ac konstruir a se iz podataka za srednju vri- 
jednost (M) i standardnu devijaciju (o). 

Srednja se vrijednost (M) nalazi kod binominalne krivulje 
na Hazenovu dijagramu na paraleli od 50%. U tab. XIV. vidimo, 
da se u granicama M -f- 3 a nalazi 99.7% svih varijanata, u gra¬ 
nicama M — 3 o samo preostalih 0.3%. 

Na apscisi ćemo potražiti vrijednost M + 3a, a na ordinati 
paralelu od 99.7%. Isto ćemo učiniti i za M — 3a za apscisu, 
a na ordinati ćemo potražiti paralelu od 0.3%. 

Točke M, m + 3o i M — 3a, koje smo tako dobili na para¬ 
lelama, leže u pravcu. 

U našem primjeru za duljinu klasa je M = 10.01 cm (str, 
24). Dakle na apscisi treba kod 10,01 cm povući ordinatu do 
paralele 50%. 

Na isti način valja tražiti i odgovarajuće točke za M + 3a 
i M —3a. U našem je primjeru a = 1,699 cm (str. 32). Dakle 
M + 3a = 10,01 + 3 X 1,699 = 15,107 cm (ordinatu treba povući 
do paralele 99.7%). 

M 3a = 10,01 + 3 X 1,699 = 4,913 cm (ordinatu treba po¬ 
vući do paralele 0.3%). 

Iz eksperimentalnih brojaka naći ćemo točke na parale¬ 
lama ovako: Na pr. u varijacionom redu za duljinu pšeničnog 
klasa u razredu 9.5—10 cm ima 56 biljaka. Sumarna je frekven¬ 
cija varijanata do granice toga razreda 216, a to je 54% od 
ukupnog broja svih 400 varijanata. Prema tome želimo li naći 
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traženu točku, na Hazenovu diagramu na si. 8. kod apscise 10 
cm povući ćemo ordinatu do paralele 54%. Na isti način dobit 
ćemo i točke na ordinatama za sve ostale razredne frekvencije. 
(Broj varijanata do granice razreda u % prikazan je u koloni 8 
tab. XIII.). Točke, dobivene tako na paralelama, spojit ćemo 
linijom. Što se bolje podudara eksperimentalna linija s teoret¬ 
skim pravcem, to se varijaciona krivulja bolje pokriva s bino- 
minalnom krivuljom t. j. slaže se s teoretskom pretpostavkom. 


Tab. XIV. 


L a t i t u d a 

% varijanata 

M + x o 

M — x a 

M ± 0-5 o 

383 

61-7 

M ± 1-0 a 

68-3 

31-7 

M -g 15 a 

866 

13 4 

M ± 2-0 a 

95-5 

45 

M ± 2-5 o 

98-8 

1-2 

M±30a 

99-7 

0-3 


V. ABNORMALNE VARIJACIONE KRIVULJE 

Varijante varijacionog reda, sastavljenog za izvjesno svoj¬ 
stvo od jednog istog genotipa, obično s_e_ prilično dobro podu¬ 
daraju s binominalnom podjelom, pa prema tome i s normal¬ 
nom krivuljom. Ali ima i svojstava s abnormalnom podjelom 
varijanata, prema tome i s abnormalnim varijacionim krivu¬ 
ljama. 

Takve su krivulje: a) asimetrične, b) simetrične s visokim 
ili niskim vrhom, c) asimetrične s visokim ili niskim vrhom ili 
d) bi-, tri- ili multimodaine t. j. sa 2, 3 vrha ili s više vrhova. 

Sve se vrste tih krivulja pojavljuju onda, ako varijacioni red 
sačinjavaju varijante različnih genotipova. U tom slučaju, da 
njihovu podjelu u varijacionom redu uzmognemo usporediti s 
binominalnom podjelom, valja nam varijaciono-statistički pro¬ 
učiti varijante svakog pojedinog genotipa. Abnormalna podjela 
varijanata može biti katkada i u genotipu. 
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Ovdje nas u prvom redu zanima abnormalna podjela vari¬ 
janata u varijacionom redu, sastavljenom od varijanata jednog 
istog genotipa. 

1. Asimetrične krivulje genotipa 

Asimetrija može biti dvovrsna i to: a) prividna ili /?) prava. 

ad a) prividna asimetrija pojavljuje se kod razrednih vari¬ 
janata, kad se srednja vrijednost ne nalazi u sredini varijacio¬ 
nog razreda t. j. kad leži na granici dvaju razreda. 

Kod cijelih varijanata pojavljuje se prividna asimetrija, ako 
srednja vrijednost nije cijeli broj ili ako srednja vrijednost leži 
između dva susjedna cijela broja. 

ad p) ima slučajeva i s pr avom asimetrijom u podjeli vari¬ 
janata jednog istog genotipa, iako su se varijante razvijale pod 
gotovo istim vanjskim uslovima. U takvim se slučajevima i 
srednja vrijednost ne nalazi u razredu s najvećim brojem vari¬ 
janata. 

Ovu je asimetriju Q u e t e 1 e t protumačio time, što faktori, 
koji uzrokuju atipičnu podjelu varijanata, djeluju jače u jednom 
smjeru nego u drugom. On je pošao od pretpostavke, da nor¬ 
malnoj podjeli varijanata odgovara formula (a + b)f|i to da 
je a — b, dok atipična podjela nastaje tek onda, kad je a > b 
ili a < b. 

No K a p t e y n je utvrdio, da ova pretpostavka vrijedi samo 
za male potencije od (a + b). Ako su potencije velike i ako 
razlika između a i b nije prevelika, ne može to utjecati na asi¬ 
metriju. 

K a p t e y n tumači asimetriju u podjeli varijanata jednog 
istog genotipa time, što organizmi aktivno reagiraju na vanjsk e 
utjecaje . To znači: ako je na pr. od dva organizma jednakog 
genotipa jedan iz nekog vanjskog razloga zaostao u rastu, a 
drugi se razvijao normalno, onda se i pod istim vanjskim uslo¬ 
vima zaostali organizam ne će jednako razviti kao normalni, 
nego tek proporcionalno prema svome rastu. 

Asimetrija može biti pozitivna ili negativna. 

Po J o h a n n s e n u je asimetrija pozitivna, ako se 
varijante u varijacionom redu protežu široko na desno od 
srednje vrijednosti. Protežu li se u svojoj širini lijevo od sred¬ 
nje vrijednosti, asimetrija je negativna. 


4 Tavčar: Biometrika u poljoprivredi 
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Jačinu asimetrije izračunavamo po raziičnim varijaciono- 
statističkim formulama. 

Za približnu orijentaciju jačine asimetrije služe nam m e- 
đ i j a n a (Med) (str. 20) i srednja vrijednost M (str. 21). 

Kod idealne podjele varijanata je Med = M, a kod asime¬ 
trije je Med > M ili Med <C M. Prema tome može kao mjerilo 
asimetrije poslužiti donekle razlika: M — Med (si. 9.). 



Ako AMć n 


SI 9. 

Točnije izračunamo jačinu asimetrije (S) pomoću Jo¬ 
ti a n nsenove formule : 

asimetrija: S = - ■: a 3 

Ovdje je a razmak između pojedinih varijanata i srednje 
vrijednosti (M). Budući da je račun jednostavniji, ako se slu¬ 
žimo aproksimativnom srednjom vrijednosti dakle polaznom 
točkom A (vidi str. 22) to treba, da u potonjem slučaju gornju 
formulu promijenimo ovako: 


asimetrija: S — (^ 1 _ 3b ^ !+2bl):pS 


Primjenu ove formule prikazat ćemo na našem primjeru 
s duljinom klasa kod pšenice. Varijacioni red za ovo svojstvo 
nalazi se na str. 18. Za aproksimativnu srednju vrijednost, koja 
nam služi kao polazna točka, uzeli smo već prije apsolutnu vri¬ 
jednost razreda s duljinom klasa od 10—10.5 cm t. j. sredinu 
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tog razreda tako, da je A = 10.25 cm. Iz toga smo izračunali, 
da je srednja vrijednost M = 10.01 cm. Budući da se ona ne 
nalazi u sredini razreda 10 i 10.25 cm, to treba da spriječimo 
prividnu asimetriju na taj način, da srednju vrijednost smje¬ 
stimo približno u sredinu razreda. Stoga treba promijeniti gra¬ 
nice razreda na pr. ovako : 6.25—6.75—7.25—7.75—8.25—8.75— 
— 9.25 i t. d. i nanovo raspodijeliti varijante. Kako nam je želja, 
da se uputimo u metodiku rada i da što preglednije i shvatljivije 
prikažemo sve operacije, to ćemo u našem slučaju izračunati 
asimetriju bez obzira na spomenuti nedostatak našeg pređašnjeg 
varijacionog reda. 

Neke su nam vrijednosti u formuli za asimetriju poznate 
već otprije. Pomoću skrižaljke na str. 23 i 24 tab. IV. i V. izraču¬ 
nali smo već b = 0.4775 prilikom računanja srednje vrijednosti, 

= SpaK _ 4,711 j lA _ b 2 = 3,398 kod računanja 
n 400 \ n 

standardne devijacije. Sada treba samo još da iz iste skrižaljke 
tab. V. i kolone 7, 8 i 2 izračunamo izraz 2pa 3 . U tu svrhu služe 
nam rubrike p i a 3 . Prema tome su: 



Ukupno | — 3157 

dakle: 2pa 3 = —3157 + 1970 = — 1187. 


Prema tome je asimetrija varijacione krivulje za duljinu 
klasa: 

s = [ - * - 3 <- °' 4775) m+ 2 ' ( “ °' 4775) '] : 

: 3,398 3 = 0,35 

* 



52 


53 


Asimetrija varijacione krivulje za broj kiasića na klasu je: 
(Tab. VIII., kolona 5, 6 i 1 na str. 27). 

S = [-' W - 3 ' ( ~ °’ 4475) • W + 2 ( ~ 0)08961467)] : 

: 25,2569 = 0,29 

Za jačinu asimetrije predlažem ovu klasifikaciju: 


Tab. XV. 


Asimetrija — S 

Jačina asimetrije 

0,0 -0,10 

0,1 -0,25 
0,25-0,5 

0,5 —više 

nema asimetrije 
mala asimetrija 
srednja asimetrija 
jaka asimetrija 


Kod nas je dakle u oba primjera pozitivna asimetrija. 

Treba konačno još napomenuti, da brojčana razlika između 
M — Med, kao aproksimativno mjerilo asimetrije nema uvijek 
isti predznak kao S. 

2. Ekscesivne krivulje (krivulje s visokim i 
krivulje s niskim vrhom) 

Kod genotipova sa simetričnom pa i kod onih s asimetrič¬ 
nom krivuljom vrh krivulje je katkad viši, a katkad niži od vrha 
idealne ili normalne krivulje. 

Takve se varijacione krivulje pojavljuju kod izvjesnih svoj¬ 
stava biljaka i životinja, te ih zovemo ekscesivnim krivuljama. 

Eksces može biti pozitivan ili negativan. U. prvom 
je slučaju vrh varijacione krivulje i z n a d, a u drugom ispod 
normalne krivulje (si. 10, krivulja A : eksces —, C : eksces -f, 
B : normalna krivulja). 

Koji sve činioci utječu na ekscesivnost pojedinih svojstava 
organizma do danas nije još pronađeno. Svakako vrijedi i ovdje 
ono, što smo naveli već kod asimetrije t. j. da dva organizma 
istog genotipa, od kojih se jedan jače razvio od drugog, ne će 
kod daljeg rasta razvijati svoja svojstva jednako. Organizmi, 



koji se razvijaju pod optimalnim prilikama rastu drugačije nego 
oni organizmi, za koje su uvjeti rasta nepovoljni. 

Za približnu orijentaciju, da li je neka varijaciona krivulja 
ekscesivna ili nije, služi nam omjer između kvartila (Q) i stan¬ 
dardne devijacije (o). Kod idealne binominalne podjele vari- 
j anata u varijacionom redu je: Q : a = 0.6745. Kad je Q:o> 
> 0.6745, . onda je eksces pozitivan, a kad je Q : o < 0.6745, 
eksces je negativan. 

Točnije možemo ocijeniti jačinu i smjer ekscesa (E) po¬ 
moću Johannsenove formule: 



I ovdje lakše i brže izvodimo računske operacije, ako gor¬ 
nju formulu prilagodimo za rad s aproksimativnom srednjom 
vrijednosti (polazna točka A). Gornja formula prelazi onda u 
ovu formulu: 

E= [(^_4b.^ + 6 l ^-3b*):,,*]-3 

Na naš primjer s duljinom klasa pšenice primijenit ćemo 
sada gornju formulu da vidimo, da li se vrh varijacione krivulje 
poklapa s vrhom normalne krivulje. 
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U gornjoj su nam formuli već iz računa srednje vrijednosti, 
standardne devijacije i asimetrije poznate slijedeće vrijednosti: 

b = — 0,4775; 2 pa 2 = 4711; 2pa 3 = — 1187 i p = 3,398 
Iz skrižaljke tab. V. na str. 24 izračunamo i 2pa 4 . U tu 
svrhu služimo se rubrikama p i a 4 , dakle: (kolone 9 i 1) 


p a< 

p a* 

93 ■ 1 

93 

74 16 

1.184 

60 • 81 

4.860 

45. 256 

11.520 

37 • 625 

23.125 

22 • 1296 

28.512 

15 • 2401 

36.015 

7 • 4096 

28.672 

2-6561 

13.122 

2 pa* 

147.103 


E “ {[-THn “ 4 X (-0.4775) X (- ^) + 

+ 6 X ( — 0,4775) 2 X ^ - 3 X ( - 0,4775)*] : 3,398*} - 3 = 

= - 0,16361 

za duljinu klasa. 

Eksces varijacione krivulje za broj klasića je pak: 
E = — 0,36043. 

Slabo pozitivno ekscesivne su one krivulje, kod kojih je 
eksces (E) manji od 0.4. Pozitivan eksces je karakterističan po 
tome, što ima u sredini varijacionog reda ili krivulje više vari- 
jana ta nego kod idealne binominalne podjele. 

ifod negativnog ek scesa ima, protivno pozitivnom, u sre¬ 
dini varijacionog reda premalo varijanata u usporedbi s 
normalnom krivuljom. 

^Varijacioni redovi s manje od 6 razreda nisu podesni za 
izračunavanje ekscesa. Kad varijacioni red sačinjava samo 6 
razreda, negativni eksces ne može biti manji od E — —2, jer u 
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tom slučaju varijacionu krivulju sačinjavaju dvije potpuno odi¬ 
jeljene krivulje. 

Eksces od E = — 2 već znači, da je varijaciona krivulja 
jako negativno ekscesivna. 

3. Varijaciona krivulja s dva ili više vrhova 

Kod bilja, a i kod životinja, ima genotipova, kojih su vari¬ 
jante za neka svojstva u varijacionom redu tako raspodijeljene, 
da čine dva, katkada i više maksimuma. Prema tome su i kri¬ 
vulje konstruirane na osnovi takvog varijacionog reda bi-, tri- 
odnosno multimodalne. Kako dolazi do takvih krivulja kod va¬ 
rijanata jednog istog genotipa, danas još ne znamo. Vjerojatno 
je ipak, da promjene u intenzitetu vanjskih faktora u pojedinim 
stadijima razvoja organizma mogu izazvati te komplicirane 
fiziološke procese. 


VI. »MODUS« — TEMELJNA TOČKA VRHA KRIVULJE 


Za karakteristiku nekog varijacionog reda ili krivulje upo¬ 
trebljava se katkada i modus. 

Potrebno je kadšto, da kod ocjene nekog svojstva znamo, 
do koje se mjere temeljna točka vrha varijacione krivulje — 
modus podudara s istom točkom u normalnoj krivulji. Modus 
je dakle najbrojnija vrijednost nekog svojstva t. j. ordinata~~š 
najviše VarijanataT " 

Kod normalne krivulje je modu s: Mo = M. Točno izra¬ 
čunavanje modusa prilično je komplicirano. 

Za monomodalne varijacione krivulje modus približno 
izračunamo pomoću Pearsonove formule : 

Mo = 3 Med —■ 2 M (M = srednja vrijednost a Med je medijana) 
Na si. 9. prikazana je asimetrična krivulja i njezine vari- 
jaciono-statističke vrijednosti srednja vrijednost (M), mediana 
(Med) i modus (Mo). 


VII. NEKE SPECIJALNE VARIJACIONO-STATISTIČKE 
METODE U POLJOPRIVREDI 

Na temelju već dosada opisanih varijaciono-statističkih 
vrijednosti izvedene su i neke specijalne statističke metode, s 
kojima se moramo upoznati, da bismo pomoću njih pravilno 
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ocijenili apsolutnu vrijednost i varijabilnost nekog svojstva kod 
čovjeka, biljke, životinje ili nežive stvari. 


I. Kako se može unaprijed izračunati broj 
varijanata potrebnih za pokus 


Kad bar približno znamo stepen varijabiliteta nekog svoj¬ 
stva, možemo već unaprijed izračunati, koliko puta treba pokus 
ponoviti t. j. koliko pokusnih životinja, parcelica ili ponavljanja 
(repeticija) treba za sortni pokus s gnojivom, dobom sjetve, 
načinom sjetve i t. d. da iz dobivenih podataka izvedemo točne 
zaključke. 

Srednja i vjerojatna pogreška srednje vrijednosti smanjuju 
se, ako je broj ponavljanja (varijanata) velik. Pokusima se usta¬ 
novilo, da je na velikoj parceli srednja i vjerojatna pogreška 
za prirod biljaka obično veća od 10% i da se ova pogreška u 
praksi ne može sniziti ispod 5%. 

Kod poljskog pokusa naročito je važno, da prije glavnog 
pokusa ustanovimo, nije li tlo pokusnog polja heterogeno. He¬ 
terogenost tla najbolje se uoči tako, da se pokusno polje zasije 
najprije jednom sortom poljoprivrednog bilja, koje jako reagira 
na razlike u kakvoći tla. 

Želimo li znati, koliko mora biti varijanata t. j. pokusnih 
parcela, životinja, imanja i t. d. da iz pokusa izvedemo ispravne 
zaključke, izračunat ćemo ili srednju ili vjerojatnu pogrešku 
jedne varijante. 

Kod žitarica je najmanja pogreška za prirod zrna po ha 
otprilike 5%. Broj varijanata, u našem slučaju pokusnih par- 
] ćela, koje trebamo, da pokus bude točan, izračunamo tako, da 
j vjerojatnu pogrešku (5%) podijelimo s drugim korijenom ne- 
| poznanice n, t. j. s brojem ponavljanja (varijanata). U drugi dio 
! jednadžbe dolazi procentualna razlika (a), koju želimo mjeriti, 
i između dvije sorte, a ovu razliku podijelimo sa omjerom između 
diferencije i vjerojatne pogreške te diferencije. Za spomenuti 
I omjer kod poljskog pokusa sa sortama uzimamo faktor 4.5. 


Prema tome je formula: , — 

V « 

Primjer: 


a 

4,5 


Želimo znati, koliko nam treba pokusnih parcelica (n) od 
svake sorte, da odredimo razliku od a = 10% u prirodu zrna 
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između sorte A i B. Dakle: 

dobijemo 


5 

1/Y 


10°/o 

4,5 ; 


10 V n = 22,5 


đaljnim računom 


n = (2,250) 2 = 5,0625 t. j. u pokusu treba svaku sortu sijati 
na 5 parcelica. U tab. XVI. izračunat je broj pokusnih parcelica 
potrebnih za signifikantnost pokusa, ako želimo mjeriti razliku 
u prirodu od a = 1 do 15%. 


Tab. XVI. 



Kod pokusa s krmljenjem goveda W ood i Stratton su 
izračunali, da je vjerojatna pogreška za prirast na težini jedne 
životinje (varijante) 14%, a 3.8 je faktor, koji odgovara vjero¬ 
jatnosti od 1 :40. 

Primjer: 

Neka se izračuna, koliko životinja treba za pokus krmlje¬ 
njem, da se još odredi razlika od 15% u prirastu na težini iz- 


i 



f 

l 
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među jedne i druge skupine pokusnih životinja. U tom se slu¬ 
čaju služimo ovom formiuom: 


14% 


dakle ^4% 


14 _ 
Vn 


]/ n 3,8 
15 \/~n = 14X3,8 

-m- 


=ja 

3,8 

12,5316 


t. j. za naš pokus trebali bismo 12 do 13 životinja. 

Mitchell i Grindley navode, da je varijacioni koefi¬ 
cijent (v) za prirast na težini u prosjeku za: ovce 21.63, svinje 
17.12 i goveda 16.73. Na temelju krmnih pokusa izvedenih na 
različnim institutima u Engleskoj i Americi, izradili su spome¬ 
nuti stručnjaci tabelu za najmanji broj životinja koji je potre¬ 
ban, da se ustanovi razlika u prirastu životinja uz dva različita 
načina krmljenja. (Vidi tab. XVII). 


Tab. XVII. 


Razlika u prirastu 
uz 2 različita obroka 

a % 

Fokus sa svinjama 

Pokus sa ovcama 

Broj životinja potreban za pokus | 

50 

i 

2 

40 

2 

3 

30 

3 

4 

20 

5 

8 

17.5 

7 

10 

15 

9 

14 

12 5 

13 

20 1 

10 

20 

31 

7.5 

36 

54 

5 

80 

121 

2.5 

317 

482 j 


2. Računsko izjednačenje varijacije priroda 
izazvane heterogenošću tla 

Usporedimo li prirode iste sorte posijane u sortnom po¬ 
kusu na 5 parcelica na jednakoj površini vidjet ćemo, da su 
prirodi jako različiti. U jednom slučaju kod sortnog pokusa s 


* 




ozimom zobi dobili smo ove prirode u zrnu: 21, 23, 20, 19, 
22 q/ha. 

Aritmetička sredina ili srednja vrijednost (M) dobivena od 
priroda sa svih 5 parcelica najbliža je pravom prirodu spome¬ 
nute pokusne sorte: 

M = 21+23 + a>+19 +j2 
5 

Prirod s pojedinih parcelica varira oko aritmetičke sredine 
M = 21 q/ha. 

Spomenutu varijaciju uvjetuju dvije skupine pogrešaka i 
to: a) slučajne pogreške ib) sistematske po¬ 
greške. 

Slučajne pogreške su ove: netočno odmjerivanje 
parcelica, netočnost kod sjetve i žetve, štete od miševa i ptica, 
ispadanje zrna kod žetve i t. d. 

Sistematske su pogreške u prvom redu pogreške 
uvjetovane heterogenošću tla. 

Najtočniji pokus je na izjednačenom tlu. Ali takova se po¬ 
kusna polja nađu vrlo rijetko. Bilo je stoga potrebno pronaći 
računske metode, koje će nam omogućiti, da na osnovu rezul¬ 
tata pokusa na heterogenim tlima donesemo pouzdane za¬ 
ključke. U tu svrhu služe ove najvažnije metode. 


a) Standardna metoda 

Kod ove metode smjestimo sorte, koje želimo istražiti na 
pr. na rodnost, između kontrolne ili standardne sorte. Isto tako 
kod pokusa s gnojivima, različnim načinom sjetve, okapanja 
i t. d. smjestimo pokusne parcele tako, da dolaze između stan¬ 
dardnih parcela. Standardna parcela može biti svaka treća- če¬ 
tvrta ili peta u redu. 

Evo primjene ove metode u praksi: 

Pravljen je sortni pokus s 4 oplemenjene sorte ozime zobi, 
posijane u jednom pojasu uz 4 repeticije. Ispred i iza svakog 
niza oplemenjenih sorti posijana je jedna domaća sorta zobi 
kao standardna ili kontrolna sorta. Oplemenjene sorte označit 
ćemo sa I, II, III, IV, a standardnu sortu sa K. 

Raspored parcela i prirode u zrnu (q/ha) ovih sorata vi¬ 
dimo u tab. XVIII. 
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Tab. XVIIf. 


Oznaka sorte i 
komparacije 

K, 

h 

ih 

III, 

IV, 

k 2 

U 

1I 2 

III 2 

iv 2 

Ka 

Parcelica broj 

i 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

a) Faktičan prirod 
q/ha 

21-3 

22-6 

197 

18-6 

23-8 

19-8 

21-3 

20-3 

19-3 

21 6 

20-7 

b) Vjerojatni pri¬ 
rod kontrolne 
l sorte 

21-3 

210 

20-7 

20-4 

20-1 

19-8 

19-98 

2016 

20-34 

20-52 

20-7 

c) Vjerojatni pri¬ 
rod pokusnih 
sorta,kad bi tlo 
bilo homogeno 

- 

22*104 

19-517 

18-727 

24-321 

- 

21-836 

20-682 

19-489 

21*621 

- 


Oznaka sorte i 
komparacije 

I. 

Ih 

IHa 

IV, 

K 4 

h 

»4 

Uh 

IV, 

Kc 

Parcelica broj 

12 

13 

14 

15 

16 

17 

18 

19 

10 

21 

a) Faktičan prirod 
! q/ha 

20-4 

13-8 

18-2 

209 

19-1 

19-9 

189 

19-7 

22-7 

21-8 

b) Vjerojatni pri¬ 
rod kontrolne 
! sorte 

20-38 

20 06 

19 74 

19-42 

19 1 

19-44 

19-78 

20-12 

20-46 

21-8 

c) Vjerojatni pri¬ 
rod pokusnih 
sorta, kad bi tlo 
bilo homogeno 

20-560 

20-273 

18-937 

22-105 

_ 

21-026 

19-626 

10-111 

22-788 

- 


Primjer: 

Vjerojatni se prirod sorte III na pokusnoj parcelici br. 9 uz 
pretpostavku, da je tlo »homogeno«, izračuna ovako: 

Najprije treba ustanoviti, kakav bi bio prirod kontrolne 
sorte (K) na parcelici br. 9. Budući da ta parcelica leži između 
parcelica br. 6 i br. 11, na kojima je rasla kontrolna sorta sa 
prirodima K 2 = 19,8 i K 3 = 20,7, znači, da se prirodi povisuju 


od parcelice br. 6 u smjeru parcelice br. 11 


za 


K, 




— 0,18 q/ha po svakoj parcelici. 


Iza parcelice br. 6 s prirodom K 2 = 19,8, parcelica br. 9 
je treća. Dakle bi prirod kontrolne sorte na parcelici br. 9 bio: 
Ko + 3 . 0,18 = 20,34 q/ha. 


Vjerojatni prirod poredbene sorte na homogenom tlu izra¬ 
čunatim tako, da faktični prirod pomnožimo s aritmetskom sre¬ 
dinom priroda svih kontrolnih parcelica i podijelimo s vjero¬ 
jatnim prirodom, koji bismo dobili, da je na odnosnoj parcelici 
rasla kontrolna sorta. 


Dakle: 


lilo x M 


19,3 X 20,54 


= 19,489 q/ha 


K 2 ^3X 0,18 19,8^-3X0,18 

Aritmetička je sredina svih' kontrolnih parcelica 

_ K,+K s -t-K s + K t +K 5 = 
m— . 5 

21,3 + 19,8 + 20,7 + 19,1 + 21,8 
— 5 

Uz pretpostavku, da je tlo homogeno, izračunat ćemo vje¬ 
rojatne prirode poredbenih sorta i na drugim parcelicama. 


20,54 q/ha 


Prosječne prirode neke sorte na nehomogenom i »ho¬ 
mogenom« tlu izračunamo tako, da sumiramo prirode s poje¬ 
dinih ^parcelica, U našem su primjeru pokusne sorte posijane na 
4 ,različite parcelice, a kontrolna sorta (K) na 5 parcelica. Dakle 
treba u svake pokusne sorte zbroj priroda na parcelicama divi- 
dirati sa 4, a zbroj priroda s parcelica, posijanih kontrolnom 
sortpm, sa 5. 

U tab. XIX. su sakupljeni prosječni prirodi zrna pojedinih 
sorata i njihove srednje pogreške i to na nehomogenom i »ho¬ 
mogenom« tlu. 


Tab. XIX. 


Oznaka sorte 

Prosječan prirod 
zrna q/ha 
nehomogeno tlo 

Prosječan prirod 
zrna q/ha 

„homogeno” tlo ■ 

Kontrolna K 

20,54 ± 0,49 

20,54 

Sorta I 

21,05 ±0,63 

21,40 ±0,23 

Sorta II 

19,68 ± 0,29 

20,03 ± 0,09 

Sorta III 

18,95 ± 0,34 

19,32 ± 0,3! 

Sorta IV 

22,25 ± 0,65 

22,71 ± 0,59 


r 
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b) Kako se kod sortnih pokusa s okopavinama računski korigira 
prirod modificiran uslijed praznih mjesta 

Kod sortnih pokusa s okopavinama posijanih ili posađenih 
u stanovitim pravilnim razmacima, neke biljke katkada uginu, • 
pa radi toga nastaju prazna mjesta. Zbog toga se biljke oko 
praznih mjesta bolje razviju, a to dakako utječe i na prirod. 

Ako želimo u takvom slučaju odrediti rodnost pokusne sorte, 
moramo korigirati vrijednost priroda, kojega dobijemo na ta¬ 
kvoj parceli. 

Ukupni prirod jedne pokusne parcele mogao bi se ^amo 
onda usporediti s prirodom druge poredbene parcele, kad ni na 
jednoj ne bi bilo praznih mjesta. 

Za izračunavanje korigirane vrijednosti priroda pokazalo 
se, da je vrlo dobra Henrichova formula. Korigiranom 
smatramo onu vrijednost priroda, koju bismo dobili, kad na 
pokusnim parcelama ne bi bilo praznih mjesta. Korigirani pri¬ 
rod p x dade se po Henrichovoj formuli izračunati ovako: 

p + b - v, . i 3 Jij jJj-- 

Pl ~ 2 * £ 4 a 

gdje je p faktični prirod, bi broj biljaka, koji bi bio na parceli, 
da nije uginula nijedna biljka, b 2 broj biljaka kod berbe. 

Za izračunavanje korigiranog priroda okopavina primje¬ 
njujem vlastitu formulu: pi = Najprije odredimo prirod 

p od prosječnih 25 biljaka, izvađenih s onih dijelova pokusne 
parcele, gdje nije bilo praznih mjesta. Prosječan prirod jedne 

biljke iznosi Korigirani prirod, koji bismo dobili, da na po¬ 
kusnoj parceli nije bilo praznih mjesta, izračunamo dakle tako, 

da prosječan prirod jedne biljke pomnožimo s normalnim 
brojem biljaka. 

3. Peterova i Besselova metoda za izračuna¬ 
vanje vjerojatne pogreške 

Kod pokusa je u poljoprivredi vrlo često već iz tehničkih 
razloga ograničen broj pokusnih varijanata. Obično moramo 
raditi s malim brojem biljaka, pokusnih parcela, životinja i si.. 
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tako da kod sortnih, gnojidbenih i drugih pokusa imamo 3, 4, 5 
a najviše deset repeticija. 

U takvim slučajevima izračunavamo vjerojatnu pogrešku 
(E. P.) jedne varijante i pogrešku srednje vrijednosti po Pete- 
rovoj ili po Besselovoj formuli. 

Peterova formula za vjerojatnu pogrešku 
jedne varijante (na pr. za jednu pokusnu parcelu ili 
životinju) glasi: 

S(±D) 

“) E,Rse ~ 0,8453 V'n (n—1) 

a za vjerojatnu pogrešku srednje vrijednosti 
(M): 

2(±D) 

P) E.P.m = 0.8453 n j /-^—y 

gdje je n broj varijanata, d zbroj diferencija pojedinih varija¬ 
nata srednje vrijednosti bez obzira na predznak, a 0,8453 je 
konstanta. 

Besselova formula za vjerojatnu pogrešku 
jedne varijante glasi: 

7) E.P. S g = 0,6745 j/ 

a za vjerojatnu pogrešku srednje vrijednosti: 

5) E.P. m = 0,6745 y n( ^ 1} 

Kako se po Peterovim i Besselovim formulama izračuna¬ 
vaju vjerojatne pogreške, pokazat će nam slijedeći 
primjer: 

Od jedne linije pšenice Prolifik, koja je s drugim različnim 
sortama pšenice uvrštena u sortni pokus na 10 parcelica, dobi¬ 
veni su ovi prinosi zrna preračunati na hektar: 
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^6 ± E. D . sg — 23,5 ± 0,676 (Besselova formula za jednu 

varijantu.) 

Za srednju vrijednost M je vjerojatna pogreška po Pete- 
rovoj formuli: 


E-P-m ■= ± 0,8453 = + 0,214, 


a po Besselovoj: 


E.P.m = ± 0,6745 [. _? I 06_ = + 0 214 

I 10 ( 10 — 1 ) - 


Prosječan je dakle prinos zrna naše Prolifik-pšenice 
M± E.P.m —23 ±0,214 (Po Peterovoj i po Besselovoj for¬ 
muli). 

U našem se slučaju vjerojatne pogreške dobivene po Pete- 
rovim formulama dobro slažu s pogreškama, koje smo izra¬ 
čunali po Besselovim formulama. Vjerojatne su pogreške izra¬ 
čunate po B e s s e 1 u točnije nego po P e t e r u. Računanje 
je po Besselovoj formuli kompliciranije i dulje nego po 
Peterovoj formuli. 


Možemo uštediti nešto vremena, ako se kod računa po 
Besselu služimo Z 6 11 e r o v i m tabelama. 


4. Modificirana Peterova metoda za srednju 
pogrešku — m M 


Spomenuta izračunavanja iziskuju u svakom slučaju dosta 
vremena. Da olakšam i ubrzam te operacije barem kod pokusa 
s 2 do 10 ponavljanja — a taj je broj čest kod poljskih pokusa 
izradio sam na osnovu Peterove formule za vjerojatnu po¬ 
grešku srednje vrijednosti tabelu za srednju pogrešku te vri¬ 
jednosti. 


Za srednju pogrešku srednje vrijednosti 

tmram frvrrnuln ♦-*-» s (+ D) 


upotrebljavam formulu m M = 


n 1/ n — 


U tab. XXI. od ove je formule izračunan izraz - 

nj/ n — 1 

za 2 do 10 varijanata (n). Prema tome treba izračunati 2 (±D) 
t- j- vrijednost, koju dobijemo, ako bez obzira na pred¬ 
znak zbrojimo razlike između srednje vrijednosti (M) 


5 Tavčar: Biometrika u poljoprivredi 
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i vrijednosti pojedinih varijanata. Upotrebu te tabele prikazat 
ćemo najbolje na ovom primjeru: 

Prirode zrna po hektaru pšenice Prolifik, posijane u sort¬ 
nom pokusu na 10 parcelica — dakle 10 varijanata — vidimo 
na str. 64. - 

Ako prirode pojedinih parcelica zbrojimo i zbroj podije- 

230 

Urno sa n = 10, dobit ćemo srednju vrijednost M = —= 
= 23 q/ha zrna. 

Odbijemo li od ove srednje vrijednosti prirod svake parce- 
lice (varijante) posebno, i ove razlike označimo sa D 1( Do do 
Dio, pa ih onda zbrojimo bez obzira na predznak (+ ili —), 
dobit ćemo zbroj 2 (± D). U našem je primjeru 2 (± D) = 7,6. 
U tab. XXI. naći ćemo za n = 10 vrijednost 0,03333; dakle je za 
M=23 srednja pogreška: m M = 7,6.0,3333 = 0,25333. 

Ako je n = 5 t. j. ako imamo u sortnom pokusu od iste 
sorte 5 parcelica (repeticija), onda je srednja pogreška srednje 
vrijednosti t. j. prosječnog priroda: him = 0,1 2 (± D). 

U tom slučaju ne trebamo tabelu, nego 2 (± D) pomno¬ 
žimo s 0,1, i to je već tražena srednja pogreška. 
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5. Studentova metoda, kojom se određuje 
opravdanost razlike između dviju srednjih 

vrijednosti 

Ako je broj varijanata u nekom pokusu malen, određujemo 
vjerojatne srednje pogreške E. P m za srednju vrijednost (M) 
najčešće po Pete r ovoj ili Besselovoj formuli. Da li je 
diferencija (D) između srednjih vrijednosti (Mi i M 2 ) varija- 

ciono-statistički opravdana, vidi se iz omjera — 

Ondje, gdje treba usporediti varijantu po varijantu jedne 
serije pokusa s odgovarajućim varijantama druge serije, da 
bismo pronašli razliku između ovih serija, služimo se često 
Studentovom metodom. Pomoću ove metode račun¬ 
ski ustanovimo, da li razlika između dvije srednje vrijednosti 
leži u granicama vjerojatnosti ili probabiliteta. 

Ovu ćemo metodu upotrebiti na pr. kod sortnih pokusa, 
koje provodimo kroz više godina; kod pokusa o djelovanju 
dvaju različnih krmiva na iste životinje; kod pokusa o djelo¬ 
vanju paše ili stajskog krmljenja na sastav mlijeka istih živo¬ 
tinja i si. 

Studentova je metoda prvi put objelodanjena god. 
1908. u časopisu »Biometrika«, a na njenu važnost i upotreblji¬ 
vost u eksperimentalnoj poljoprivredi upozorio je god. 1923. 
američki učenjak prof. dr. Love. U Americi se ova metoda 
mnogo upotrebljava na poljoprivrednim pokusnim stanicama, 
dok je u Evropi još i dan danas slabo poznata. 

Da je što bolje upoznamo, prikazat ćemo njenu primjenu 
kod sortnog pokusa sa sojom. Kroz pet godina vršili smo sortne 
pokuse, kako bismo pronašli, koja od 4 sorte soje oplemenjene 
u našem Zavodu za oplemenjivanje bilja sadržava najviše masti. 
U pokusu je svake godine svaka od sorata A, B, C i D posijana 
u 6 repeticija na parcelicama od 5 m duljine i u 5 redova širine. 
Nacrt sijanja tih sorata, te kontrolne sorte K bio je ovaj: (si. 11) 





Od svake sorte uzet je sa svake od 6 parcelica prosječan 
uzorak sjemena, i nakon toga se određivao sadržaj masti u sva¬ 
kom uzorku 3 puta. Interesiralo nas je, koja od sorata A i B, 
sadržava nasljedno uvjetovan veći postotak masti. 

Sadržaj masti u tim sortama u pojedinim godinama pokusa 
prikazan je u donjoj tabeli: 



Godina 

Sadržaj masti u sjemenu 

% 

A —B = D 

D - Md 


Sorta A 

Sorta B 



1934 

19,25 

18,36 

0,89 

0,41 

1935 

18,77 

17,84 

0,93 

0,45 

1936 

17,58 

17,67 

-0,09 

-0,57 

1937 

20,36 

19,65 

0,71 

0,23 

1938 

18,09 

17,34 

0,75 

0,27 

1939 

18,45 

18,76 

-0,31 

-0,79 


2č_ 109,62 

_ n 6 

= 18,75% = 18,27 o/o 



S(D—Md) 2 
= 1,4454 


= 0,49081 


z= Mp_ = Md: 

x 


a odgovarajuće su vrijednosti izračunate po ovim formulama: 
M d = M==^---=0,48, 

x = y S(D ~ Mp)2 -= = 0,49081 

v Md m 1/-S(D-MdV _ 0,48 —p 0770 

Z = ^T = Md - |/-n-~ 0,49081 - U ’ y//9 

Najprije se izračuna srednja vrijednost (Md) diferencije 
A — B = D; iza toga treba za svaku godinu pronaći razliku 
između D i Md, a onda (D — Md) 2 . Sumu (D — Md) 2 podije¬ 
limo s brojem svih varijanata, u našem slučaju s brojem godina 
(n), te iz konačnog broja izvadimo drugi korijen, a to daje X. 

Pomoću Md i X izračunamo Z = —1 a iz njega u tabeli XXII. 

X 

pronalazimo vjerojatnost P. Vjerojatnost, koja odgovara 
vrijednosti Z, naći ćemo u tab. XXII., u rubrici Z i n (broj vari¬ 
janata). 


1.45 

4.20 

1.5 

4.34 

1.55 

4.48 

1.6 

4.62 

1.65 

4.76 

1.7 

4.91 

1.75 

5.05 

1.8 

5.20 

1.85 

5.34 

1.9 

5.49 

1.95 

5.63 

2.0 

5.78 

205 

5.92 

2.1 

6.07 

2.15 

6.21 

2.2 

6.36 

2.25 

651 

2.3 

6.66 

2.35 

6.81 

2.4 

6.96 

2.45 

7.10 

2.5 

7.26 

2.55 

7.40 

2.6 

7.55 

2.65 

7.70 

2.7 

7.86 


> 













Tab. XXII a 


z. 

n = 2 

n = 3 

n = 4 

n = 5 

n*= 6 

n = l 

ti —8 

H = 9 

' 71 = 10 

1 71 = 11 

71 = 12 71 = 13 71= 14 71= 15 71= 16 71= 17 

.1 

1.14 

1.22 

1.29 

1.35 

1.40 

1.46 

150 

1.55 

1.59 

1.64 

1.68 

1 72 

1 76 

1.80 

1.84 

1.88 

.15 

1.21 

1.35 

1.46 

1.56 

1.66 

1.75 

183 

1.92 

2.00 

2.08 

2.16 

2.24 

2.31 

2.39 

2.47 

254 

.2 

1.29 

1.49 

1.60 

1.82 

1.97 

2 12 

2 26 

2.41 

2.55 

2.70 

2 84 

2.99 

3.14 

3.29 

3.44 

3 60 

.25 

1.37 

1.64 

1.88 

2.10 

2.32 

2.54 

2.75 

2.97 

3.19 

341 

3.64 

3.87 

4.11 

4.36 

4.60 

4.86 

.3 

1.46 

1.81 

2.13 

2.44 

2.76 

3.08 

3.41 

3.75 

4.11 

4.48 

4.86 

5.27 

5.69 

6.13 

6 59 

7.08 

.35 

1.54 

1.98 

2.40 

2.81 

3 24 

3 68 

4 14 

462 

5.13 

567 

6 24 

6.84 

7.47 

8.15 

8.86 

9.62 

4 

1.64 

2.18 

2.72 

3.27 

3.85 

4.48 

5 15 

5.88 

6.67 

7.53 

8.45 

9.47 

10.6 

11.8 

13.1 

14.5 

.45 

1.74 

2.39 

3 05 

3.75 

4.51 

5.33 

6.24 

7.24 

8.34 

9.54 

10.9 

12.3 

13.9 

15 7 

17.7 

19.8 

.5 

1.83 

2 62 

3.44 

4.35 

5.36 

6 50 

7.80 

9.26 

10.9 

128 

149 

17.3 

20.1 

23.3 

26.8 

30.8 

.55 

1.94 

2.85 

3.85 

4.97 

625 

7.72 

9 42 

114 

13.6 

16 2 

19.2 

22.7 

26 7 

31.3 

36.5 

42.7 

6 

2.05 

3.12 

4 33 

5.75 

7.42 

9 42 

11.8 

14.6 

18.0 

22.0 

26.8 

324 

39.3 

47.3 

56.8 

68.4 

.65 

2.16 

3.39 

4.82 

6.54 

8.62 

11.2 

14.2 

17.9 

22.5 

27.9 

34.6 

42.5 

52.2 

63.9 

77.7 

95 2 

.7 

2.27 

3.69 

5.41 

7.55 

10.2 

13.6 

17.8 

23.1 

29.8 

38.1 

48 5 

61.5 

77.7 

99.0 

124. 

155. 

.75 

2.39 

3 99 

5.99 

8.55 

11.8 

16.0 

21.4 

28.3 

37.2 

48.3 

62.7 

81.0 

104. 

134. 

171. 

216. 

.8 

2.51 

4.33 

6.70 

9.82 

14.0 

19.5 

26 8 

36 5 

49.3 

66.1 

88.3 

118. 

158. 

207. 

277. 

356. 

*85 

2.62 

4.66 

7 39 

11.1 

16.1 

22.9 

32.1 

44 5 

61.1 

83.7 

114. 

155. 

207. 

277. 

369. 

499. 

.9 

2.75 

5.04. 

8.22 

12.7 

18.9 

27.7 

40.0 

57.1 

81.0 

114. 

160. 

226. 

311. 

434. 

587. 

832. 

95 

2.87 

5.41 

9.03 

142 

21 7 

32.4 

47.8 

69.4 

100. 

144. 

207. 

293. 

416. 

587. 

832. 

1110. 

1.0 

3.00 

5.83 

10.0 

16.2 

25 5 

39.2 

59.2 

89.1 

132. 

195. 

293. 

434. 

624. 

908. 

1428. 

1999. 

1.05 

3.12 

6.24 

11.0 

18.2 

29.1 

45.7 

70 4 

108. 

163. 

243. 

369. 

555. 

832. 

1249. 

1999. 

2499. 

1.1 

3.26 

6.69 

12.1 

20.6 

34.0 

54.9 

87.5 

138. 

216. 

332. 

525. 

832. 

1249. 

1999. 

3332. 

4999. 

1.15 

3.39 

7.13 

13.2 

22.9 

38.7 

63.5 

103. 

166. 

262. 

416. 

666. 

939. 

1666. 

2499. 

3332. 

4999. 

1.2 

3.52 

7.63 

14.5 

25.9 

44.9 

75 9 

127. 

212. 

344. 

555. 

908. 

1428. 

2499. 

3332. 

4999. 

9999. 

1.25 

3.65 

8.11 

15.7 

28.8 

50.8 

87.5 

151. 

255. 

416. 

713. 

1110. 

1666. 

3332. 

4999. 

4999. 


1.3 

3.79 

8.64 

17.2 

32.3 

58 5 

104. 

184. 

322. 

555. 

999. 

1666. 

2499. 

4999. 

9999. 

9999. 


1.35 

3.92 

9.16 

18.6 

35.8 

66.1 

119. 

216. 

384. 

666. 

1249. 

1999. 

3332. 

4999. 




1.4 

4.07 

9.74 

20.3 

40.0 

75.9 

142. 

262. 

475. 

908. 

1666. 

3332. 

4999. 

9999. 




1.45 

4.20 

10.3 

21.9 

44.0 

85.2 

163. 

302. 

555. 

1110. 

1999. 

3332. 

4999. 





1.5 

4.34 

10.9 

23.9 

49.0 

98.0 

191. 

369. 

713. 

1428. 

2499. 

4999. 

9999. 





1.55 

4.48 

11.5 

25 7 

53.9 

109. 

216. 

434. 

832. 

1666. 

2499. 

4999. 






1.6 

4.62 

12.2 

27.7 

60.0 

124. 

255. 

525. 

9:>9. 

1999. 

3332. 

9999. 






1.65 

4.76 

12.8 

29.8 

65.7 

138. 

285. 

587. 

1110. 

2499. 

3332. 







1.7 

4.91 

13.5 

32.2 

72.5 

158. 

332. 

713. 

1428. 

3332. 

4999. 







1.75 

5.05 

14.2 

34.5 

79.0 

174. 

369. 

832. 

1666. 

3332. 

4999. 







1.8 

5.20 

14.9 

37.0 

86.7 

199. 

434. 

999. 

1999. 

4999. 

9999. 







1.85 

5.34 

15.6 

39.5 

94 2 

216. 

499. 

1110. 

2499. 

4999. 








1.9 

5.49 

16.4 

42.5 

103. 

243. 

587. 

1249 

3332. 

9999. 








1.95 

5.63 

17.1 

45.1 

111. 

269. 

624. 

1428. 

3332. 









2.0 

5.78 

17.9 

48.3 

122. 

302. 

713. 

1666. 

4999. 









2.05 

5.92 

18.7 

51.4 

132. 

332. 

768. 

1999. 

4999, 









2.1 

6.07 

19.6 

54.9 

144. 

369. 

908. 

2499. 

4999. 









2.15 

6.21 

20.4 

58.2 

155. 

399. 

999. 

2499. 

4999. 









2.2 

6.36 

21.3 

61.9 

168. 

454. 

1249. 

3332. 

9999. 









2.25 

651 

22.1 

65.2 

181. 

499. 

1249. 

3332. 










2.3 

6.66 

23.1 

69.4 

199. 

555. 

1428. 

4999. 










2.35 

6.81 

24.0 

73.1 

212. 

587. 

1666. 

4999. 










2.4 

6.96 

25.0 

77.7 

232. 

666. 

1999. 

4999. 










2.45 

7.10 

25.9 

81.6 

249. 

713. 

1999. 

4999. 










2.5 

7.26 

26.9 

86.7 

269. 

768. 

2499. 

4999. 










2.55 

7.40 

27.9 

91.6 

285. 

832. 

2499. 

4999. 










2.6 

7.55 

29.0 

97.0 

302. 

908. 

2499. 

9999. 










2.65 

7.70 

30.1 

102. 

322. 

999. 

2499. 











2.7 

7.86 

31.2 

108. 

356. 

1110. 

3332. 
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Tab. XXII c. 


z. 

n = 2 

n = 3 

II 

c 

n =5 

n = 6 

n = 7 

2.75 

8.00 

32.2 

113. 

369. 

1110. 

3332. | 

28 

8.16 

33.4 

118. 

399. 

1249. 

4999. 

2.85 

8.30 

34.5 

124. 

416. 

1249. 

4999. 

2.9 

8.46 

35.6 

131. 

454. 

1428. 

4999. 

2.95 

8.61 

36.7 

136. 

475. 

1428. 

4999. 

3.0 

8.77 

37.9 

144. 

525. 

1666. 

4999. 

3.05 

8.91 

39.2 

151. 

555. 

1P66. 

4999. 

3.1 

9.07 

40.5 

158. 

587. 

1999. 

9999. 

3.15 

9.21 

41.6 

163. 

587. 

1999. 


3.2 

9.37 

42.9 

171. 

624. 

2499. 


3.25 

9.52 

44.0 

178. 

666. 

2499. 


3.3 

9.67 

45.5 

188. 

713. 

2499. 


3.35 

9.82 

46.8 

195. 

768. 

2499. 


3.4 

9.98 

48.3 

203. 

832. 

3332. 


3.45 

10.1 

49.5 

212. 

832. 

3332. 


3.5 

10.3 

51.1 

221. 

908. 

3332. 

- • 

4.0 

11.8 

66.1 

322. 

1666. 

4999. 

* 

4.5 

13.3 

83.0 

434. 

2499. 

4999. 

V 

5.0 

14.9 

102. 

624. 

3332. 

9999. 

ft 

5.5 

16.5 

122. 

832. 

4999. 



6.0 

18.0 

146. 

999. 

9999. 


> 

65 

19.6 

171. 

1249. 




7.0 

21.1 

199. 

1666. 




7.5 

22.7 

226. 

1999. 




8.0 

24.3 

255. 

2499. 




8.5 

25.8 

293. 

2499. 



1 

9.0 

27.4 

322. 

3332. 



1 

9.5 

28.9 

369. 

3332. 




10.0 

30.5 

399. 

4999. 




15.0 

46.2 

908. 

9999. 




20.0 

61.9 

1666. 





25.0 

77.7 

2499. 





30.0 

' 93.3 

3332. 





35.0 

109. 

4999. 





40.0 

124. 

4999. 




j 

45.0 

140. 

9999. 





50.0 

153. 






60.0 

188. 





* 

70.0 

221. 





I 

80.0 

249. 





1 

90.0 

285. 






100.0 

311. 





1 

120.0 

369. 





1 

140.0 

434. 






150.0 

475. 






160.0 

499. 





t 

180.0 

555. 






200.0 

624. 






250.0 

768. 





; 

300.0 

908. 






350.0 

1110. 





4? 

400.0 

1249. 






450.0 

1428. 






500.0 

1666. 






600.0 

1999. 






700.0 

1999. 





\ 

1000.0 

3332. 






1500.0 

4999. 






2000.0 

4999. 






3000.0 

9999. 
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U našem je primjeru Z = 0,9779, an = 6. Vjerojat¬ 
nost — P — se nađe interpretacijom u tabeli XXII. iz Z = 0,95, 
dakle 21,7 i za Z = 1,0... 25,5. 

Prema tome za Z = 0,9779 vjerojatnost — P — je 1 :23,2, 
da sorta A sadrži više masti od sorte B. 

U praksi se drži, da je već vjerojatnost P ... 1 :40 stati¬ 
stički opravdana. Prema tome u našem primjeru nema oprav¬ 
dane razlike u sadržaju masti u sortama A i B. 

Vidimo dakle, da tamo, gdje treba varijante jedne serije 
pokusa usporediti s odgovarajućim varijantama druge serije 
pokusa, najbolje odgovara vjerojatnost P, koja se izračuna po 
Studentovoj metodi. Ona nam je u tom slučaju ujedno 
vrlo dobra podloga, pomoću koje izvodimo pravilne zaključke. 

6. Mjerenje točnosti neke razlike pomoću 
Kollerove tabele i dijagrama 

Jedno isto svojstvo očituje se u različitim genotipovima 
vrlo različito. Opravdanost razlike između dva genotipa za 
jedno isto svojstvo možemo ustanoviti i pomoću Kollerovih ta¬ 
bela (vidi tab. X., str. 37 i dijagram si. 6). 

Postupak kod izračunavanja najbolje će nam ilustrirati ovaj 
primjer: 

Pokusom smo htjeli ustanoviti, koja od dvije sorte ozime 
repice proizvedene u 5 različnih klimatskih i edafskih područja 
(dakle n = 5) daje veći prirod sjemena. Od sorte A dobiveni 
su prirodi i to u kraju: a ... 14,1; b ... 13,2; c ... 12,5; d ... 15,1 
i e... 13,6 q/ha, a od sorte B prirodi u kraju a... 13,8; 
b ... 13,5; c ... 12,1; d ... 14,3 i e ... 13,0 q/ha sjemena. 

Da li sorta A bolje rodi od sorte B? 

Da se uoči utjecaj klime i tla svakog proizvodnog područja, 
ustanovljena je za svako takvo područje razlika u prirodu sje¬ 
mena između sorte A i B. Razlika u prinosu sjemena jedne i 
druge sorte, dakle A — B, je za proizvodno područje a... 14,1 
— 13,8 = 0,3 q/ha; za područje b... 13,2 — 13,5 = — 0,3; 
c ... 12,5 — 12,1 = 0,4; d ... 15,1 — 14,3 = 0,8 i područje 
e ... 13,6 — 13,0 = 0,6 q/ha. 

Srednja je vrijednost gornjih razlika: M D — ^ ^ = -^p = 

= 0,36 u korist sorte A. Da li je ova razlika statistički oprav¬ 
dana? 
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Razlika 

A - B = D 

D2 

a 0,3 

0,09 

b -0,3 

0,09 

c 0,4 

0,16 

d 0,8 

0,64 

e 0,6 

0.36 

SD = 1,8 

SD 2 = 1,34 


Treba izračunati srednju 


pogrešku: m = 


o = 


m = 


SD 2 
n - 

= 0,517 
0,517 


1/T 


1,34 


0,231 


Za n — 1 = 5 — 1 —4 slobodna varijanta je t = 6,62 
(vidi tab. X.). Dakle za m = 0,231 i srednju vrijednost razlike 
(Md) = 0,36 su granice točnosti mXt=0,231X6,62 = 1,52922, a 
to znači, da je ova vrijednost veća od srednje vrijednosti. Prema 
tome razlika od 0,36 nije statistički opravdana, i nema dokaza, 
da je sorta A bolja od sorte B. 

Ako je broj slobodnih varijanata veći od 10, ustanovlju¬ 
jemo t pomoću Kollerova dijagrama (si. 6, str. 37). 

7. Srednja i vjerojatna pogreška diferencije 

i s u m e 

U poljoprivrednoj nam je praksi često od važnosti da 
znamo, jesu li brojčane razlike, a katkada i zbrojevi (sume) 
srednjih vrijednosti nekog svojstva varijaciono-statistički oprav¬ 
dane ili nisu. 

Da bismo to utvrdili, poslužit ćemo se zakonom o preno¬ 
šenju pogreške. Po tomu je zakonu srednja pogreška di¬ 
ferencije (mp) ili sume (m?) jednaka drugom korijenu iz sume 
kvadrata pojedinih srednjih vrijednosti, dakle: 

m D = l/m 1 2 + m 2 2 

rn s = |/m 1 2 -t-m i , 2 -|-m s 2 -K.m x 2 

a) Srednja i vjerojatna pogreška diferencije 

Kod uspoređivanja rezultata različitih pokusa najčešće 
treba ustanoviti, da li je razlika između dvije srednje vrijednosti 
nekog svojstva varijaciono-statistički opravdana. To možemo 


i 
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odrediti na različite načine, a jedan od njih je pomoću srednje 
pogreške diferencije. 

Srednja pogreška diferencije je drugi ko¬ 
rijen iz zbroja kvadrata srednjih pogrešaka 
(m x i m 2 ) srednjih vrijednosti (Mi i M 2 ). 

Dakle: m D = ]/m, 2 -f m 2 2 

Primjer: a) U sortnom pokusu bile su uz druge sorte 
pšenice i sorte A i B, od kojih je svaka uzgajana na 10 pokusnih 
parcelica iste veličine. Prosječni prirodi zrna (M) po hektaru 
pšenice A i B s njihovim srednjim pogreškama bili su ovi: 

Sorta A.Mi ± mi_22,5 q/ha ± 0,659 

Sorta B .M 2 ± m 2 _20,1 q/ha ± 0,732 

Razlika: Mi — M 2 = D .... 2,4 q/ha 

Dakle je razlika u prirodu između jedne i druge sorte 
2,4 q/ha. 

Srednja je pogreška ove razlike: 

m D = l/m, 2 + m 2 2 = ]/0,659 2 + 0,732 2 = 0,9944 

Kad je omjer —— = 3 ili veći od 3, onda smatramo, da je 

mp 

razlika (D) varijaciono-statistički opravdana, jer je u tom slu¬ 
čaju vjerojatnost 1 :22,5 da je ona sorta, od koje smo dobili 
veći prirod zrna, uistinu rodnija. Spomenuti omjer vrijedi kao 
mjerilo i kod različnih drugih pokusa u poljoprivredi na pr. kod 
pokusa s gnojivima, razmacima biljaka u polju, dobom sjetve, 
mliječnošću krava i si. 

D 24 

U našem je primjeru ~^j _ = ' 09^44 ~ 2,4135. Kako je ovaj 

omjer manji od 3 slijedi, da je razlika između srednjih vrijed¬ 
nosti Mj i M 2 premalena, a njihove srednje pogreške mi i m 2 
prevelike, a da bi se moglo ustvrditi, da je sorta pšenice A 
uistinu bolja od sorte B. 

Primjer: b) Prof. dr. Stj. Filipović je proučavao sadržaj 
vode u mlijeku u zimsko doba i uspoređivao ga sa sadržajem 
vode u mlijeku u ljetno doba. Pokuse je vršio s mlijekom krava 
fakultetskog dobra Maksimir. 

U zimsko doba mlijeko je sadržavalo u prosjeku 86,62% ± 
± 0,13 vode, a u ljetno doba 85,72% ± 0,1 vode. Dakle je razlika 
u sadržaju vode ljetnog i zimskog mlijeka 0,9%. 
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Zanima nas, je li ta razlika varijaciono-statistički oprav¬ 
dana. To ćemo ustanoviti na ovaj način: 

Zimsko mlijeko 86,62% ± 0,13 vode 
Ljetno mlijeko 85,72% ±0,1 vode 

Razlika D.0,9% vode 


mo = |/ m i 2 + m a 2 = j/0,13 2 + 0,1 2 = 0,164 


_D^ 

mo 


0,9 

0,164 


= 5,48 


Budući da je taj omjer veći od 3, možemo sa sigurnošću 
ustvrditi, da je zimsko mlijeko u konkretnom slučaju sadržalo 
za 0,9% više vode od ljetnog mlijeka. 

Ako mjesto srednjih pogrešaka m upotrebimo vjerojatne 
pogreške E. P., dakle Mi ± E. P x i M 2 ± E. P 2 , onda računamo 
vjerojatnu pogrešku (E. P.d) razlike u prirodu (D = 
= M x — M 2 ) po ovoj formuli: 

E.P. d = J/E.P/+EP., 2 

Pokusima je ustanovljeno, da je kod sitnih žitarica s razli¬ 
kom u prirodu zrna od 2—3 q po hektaru između jedne i druge 

sorte omjer manji od 3. To znači, da je razlika premalena, 
ma 

a da bi bila varijaciono-statistički opravdana. 

Kod poljskih je pokusa srednja odnosno vjerojatna po¬ 
greška srednjih vrijednosti obično veća od 5%. 

Gore spomenute formule za srednju i vjerojatnu pogrešku 
diferencije vrijede samo onda, ako istraživana svojstva nisu u 
međusobnoj povezanosti (korelaciji). Ako su ova svojstva kore- 
lativno vezana, eliminiramo pogreške pomoću metode najma¬ 
njih kvadrata. Formula za srednju pogrešku diferen¬ 
cije glasi za taj slučaj ovako: 


m D = J/m, 2 + m 2 2 —2r m, • m 2 , a za vjerojatnu pogrešku dife¬ 


rencije: 

E.P.d = 1/E.P., 2 + E.P., 2 - 2rE.P,-E.P 2 
r je korelacioni koeficijent (vidi str. 100). 


b) Srednja i vjerojatna pogreška prosjeka 
srednjih vrijednosti m p i E. P. p 
Srednju pogrešku prosjeka od više srednjih vrijednosti 
izračunamo tako, da drugi korijen iz zbroja kvadrata pojedinih 
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srednjih pogrešaka podijelimo s ukupnim brojem svih srednjih 
vrijednosti, dakle 

_ j/m , 2 + m 2 2 -f m 3 2 +~.. ■ m x 2 
p x 

Primjenu gornje formule pokazat ćemo najbolje na slije¬ 
dećem primjeru: 

Proveden je sortni pokus sa 6 različnih sorta ozimog ječma 
(Hordeum vulgare). Svaka je sorta u pokusu posijana na 5 par- 
celica. Iz priroda zrna preračunatog na hektar, izračunan je 
prosječan prinos zrna i srednja pogreška za svaku sortu napose. 
(Tab. XXIII.). Zbroj srednjih vrijednosti je 2 M = 127,26. 

Prosječni prirod zrna svih sorata: =.—g—• = 

= 21,21 q/ha. 

Srednja pogreška prosjeka srednjih vrijednosti je prema 
tome: 


„ V Em 2 1/60,8143 7,798 onn 

m P x — 6 — 6 _ ’ 

Na temelju tih brojeva možemo usporediti srednje vri¬ 
jednosti pojedinih sorata s prosjekom sred¬ 
njih vrijednosti svih sorata u pokusu. Razlike 
između priroda svake pojedine sorte i prosjeka priroda svih 
pokusnih sorata zajedno opterećene su jednom izvjesnom sred¬ 
njom pogreškom. Nju ćemo izračunati po formuli mD == 


= j/m 2 + m 2 p (vidi rubriku 6 u tab. XXIII.). 

Razlike u prirodu pojedinih sorata u nekom sortnom po¬ 
kusu obično nisu velike. Dakako da valjane zaključke možemo 
izvoditi samo onda, kad su srednje vrijednosti (M) veće od 
2 m %. Na temelju ovog postotka možemo pravilnije ocijeniti 
vrijednost pojedinih sorata i pomoću spomenute metode gru¬ 
pirati ih po visini prosječnog priroda. 

Iz tab. XXIII. vidi se, da je prosječni prirod sorte br. 1 naj¬ 
veći; iza njega dolazi sorta br. 2. Kako je sorta br. 1 opterećena 
dosta velikom srednjom pogreškom, nije sigurno, da li ona tu 
rodnost dalje nasljeđuje. 

Uzimajući u obzir nasljednu vrijednost u našem pokusu 
istraženih sorata, grupirat ćemo ih po visini njihova priroda 
najbolje tako, da izračunamo t. zv. koeficijent kole¬ 
banja (k). 




77 

Taj koeficijent izražava omjer između razlike (D) u pri¬ 
rodu jedne izvjesne sorte i prosjeka priroda svih sorata, te 

srednje pogreške niD spomenute razlike, dakle — = k. 

m d 

Iz toga nađemo slijed pokusnih sorata, prikazan u tab. 
XXIII. 

8, V a r i j a n c a 

a) Općenito 

Po dosadašnjim se principima eksperimentalnog rada sma¬ 
tralo, da je pokus točno proveden samo onda, ako je na nj isto¬ 
vremeno utjecao što manji broj činilaca. Nastojalo se urediti 
pokus tako, da budu osim onog faktora, koji se želio proučiti, 
svi ostali faktori po mogućnosti konstantni. Ali u prirodi na 
neko svojstvo organizma utječe redovito vrlo velik broj čini¬ 
laca. Njihovo je djelovanje veoma komplicirano. Oni utječu na 
svojstvo paralelno i komplementarno. Stoga je važno, da neko 
svojstvo pokusnog organizma istražimo pod utjecajem što ve¬ 
ćeg broja tih činilaca, jer su u tom slučaju zaključci mnogo 
točniji i općenitiji. To se najbolje provodi kompleksnom ili 
faktorijalnom podjelom pokusa. Ako želimo na pr. ustanoviti, 
koje sorte pšenice od onih, što se u većem broju proizvode u 
izvjesnom kraju, tom kraju najbolje i odgovaraju, onda je slika 
točnija, ako u sortni pokus uklopimo i vrstu tla (lagano i teško), 
gustoću sjetve i dobu sjetve. 

Kod sortnih pokusa, pokusa s gnojenjem, ili pokusa s dubi¬ 
nom sijanja, gustoćom sjetve, dobom sjetve i t. d., pokusne 
smo parcelice dosada obično sistematski poredali u skupine u 
obliku šahovskog polja. 

Na pr. trebalo je načiniti sortni pokus s 10 sorata pšenice 
(I—X) u peterostrukoj komparaciji i u 5 pojaseva. Sistematski 
raspored bio bi obično ovakav: (SI. 12). 
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U tom se slučaju za svaku sortu posebno izračunala sred¬ 
nja vrijednost za prirod i srednja pogreška svake srednje vri¬ 
jednosti napose. Na temelju tako dobivenih pogrešaka izraču¬ 
nali smo srednje pogreške diferencije između dvije i dvije sorte. 
Ove srednje pogreške diferencije služile su nam kao mjerilo da 
ustanovimo, je li razlika u prosječnom prirodu (srednja vri¬ 
jednost) ma koje dvije sorte sortnog pokusa varijaciono-stati- 
stički opravdana ili nije. 

Kako su u gornjem slučaju pokusne sorte bile sistematski 
poredane, srednje se pogreške srednjih vrijednosti dobivenih 
priroda ne će posve podudarati sa srednjim pogreškama onih 
srednjih vrijednosti priroda, koje bismo dobili, da je naših 10 
sorata na parceli bilo raspodijeljeno bez sistema. 


F i s h e r je ustanovio, da nam pokus, u kojemu su sorte 
poredane bez sistema, daje točniju sliku o gospodarskoj vrijed¬ 
nosti pojedine sorte, negoli onaj sa sistematskom podjelom 
parcelica. 

Za proučavanje pojedinih svojstava kod nesistemat¬ 
ske podjele pokusnih parcelica najbolja je Fisherova 
analiza varijance, koja se u novije vrijeme mnogo upo¬ 
trebljava. 

Pod varijancom razumijevamo kvadratnu vrijednost stan¬ 
dardne devijacije: 

Ako na pr. želimo proučiti jedno svojstvo nekog određe¬ 
nog uzgoja goveda na osnovu istraživanja pojedinih individua 
(varijanata) kojih su vrijednosti mjerenog svojstva Vi, V 2 , 
V 3 ... V„, onda je prosječna vrijednost ili srednja vrijednost 
toga svojstva jednaka zbroju apsolutnih vrijednosti pojedinih 
varijanata (Vi + V 2 + V 3 ... V„) dividiranom s brojem svih 
istraženih individua — varijanata (n). Dakle: 


M = —, dok je: a 2 = 
n J 


S (V — M) 2 


= kvadrat standardne devijacije. 

Ako raspolažemo samo malim brojem varijanata, onda se 


služimo formulom: a 2 — 


_ S (V— M) 2 
n — 1 


Varijancu analiziramo. 

Kod toga izraz 2 (V — M) 2 razdijelimo na pojedine sku¬ 
pine, kojih varijaciju uvjetuju različni činioci. 
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Izraz za broj varijanata (n) smanjenih za 1, dakle (n_1), 

razdijelimo već prema broju skupina na stupnjeve neovisnih 
varijanata, (degrees of freedom). To će reći, da možemo od n 
varijanata, koje se nalaze u zbroju, slobodno izabrati samo 
(n—1), jer n-ta varijanta nije više slobodna, nego se ravna po 
ostalim (n—1) varijantama. 

'Treba voditi računa i o tome, da svaka vrijednost koja ka¬ 
rakterizira neki uzorak i koja se upotrebljava za izračunavanje 
nekih drugih vrijednosti, smanjuje i broj neovisnih varijanata. 

Tako je na pr. za varijacioni koeficijent v = ^ ° broj 

neovisnih varijanata 2n — n — 2, jer prije toga treba izraču¬ 
nati M i a. 

Da u sortnom pokusu imamo 10 sorata, koje smo uz 5-ero- 
struku repeticiju smjestili u 5 pojaseva, i bez sistema raspore- 
dali na 50 parcelica, onda bi slučajna podjela ovih 10 sorata 
mogla biti i ovakova: (si. 13). 


SI. 13 



Kako je n — 50 parcelica ili varijanata, to imamo ukupno 
n — 1 = 50 — 1 = 49 slobodnih varijanata. 

Kod ove se podjele pogreška ukupnog pokusa sastoji od 
pogrešaka u prirodu pojedinih parcela. 

U analizi varijance broj se neovisnih varijanata rasporedi 
ovako: (Tab. XXIV.) 

Premda je princip, da svaki materijal, koji želimo varija- 
ciono-statistički istražiti, mora predstavljati slučajni i neizabrani 
uzorak jednog varijacionog reda, ipak treba po mogućnosti 
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Tab. XXIV. 



Broj neovisnih 
varijanata j 

a) broj parcela: n a = 50 

b) broj sorata: nb =10 

n a —1=50—1= 49 

n b -1 = 10-1= 9 

Razlika = pogreška: (ti a — 1) — (nb — 1) — 40 


ograničiti pogreške nastale heterogenošću tla. To možemo uči¬ 
niti na taj način, da pokusne parcele svrstamo u skupine poja- 
seva (blocks) i to tako, da se u svakom pojasu nađe svih 10 
sorata, ali bez sistema. Na pr. (si. 14). 


SI. 14. 
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Ovom smo prilikom ukupnu varijancu razdijelili na sortnu 
varijancu i na pojasnu (blockvariance). 

Pojasna je varijanca uvjetovana heterogenošću tla. 

Kod izvjesne je ukupne varijance pogreška u pokusu to 
manja, što je veća pogreška u pojasu. Pogreška za sortu mora 
biti dakako što manja. Prema tome treba rasporediti parcelice 
tako, da može varijanca između pojasa biti veća, a u samom 
pojasu što manja. 

Gdje je 10 sorata u 5-erostrukoj komparaciji smješteno u 
5 pojaseva, a u svakom se pojasu nalazi svih deset sorata pore¬ 
danih bez sistema (si. 14), pogreška je u pokusu i broj neovisnih 
varijanata ovaj: (Tab. XXV.). 


Tab. XXV. 



Broj neovisnih varijanata 

a) broj parcela: n a = 50 

b) broj sorata : nb = 10 

c) broj pojasa : nc = 5 

n- 1=50-1 =49 

n — 1 = 10—1= 9 

n — 1 = 5—1= 4 



Razlika = pogreška 

(n a — 1) — [(nb — 1) + (nc — 1)] = 36 


Kod ove je raspodjele parcela i sorata pogreška manja (36) 
nego u prvoj raspodjeli (40) (vidi str. 80). 

Pogreška u pokusu još je manja, ako pokusne parcelice 
rasporedimo u obliku kvadrata (latin square). Najbolji su pak 
oni pokusi, kod kojih je broj sorata odnosno pokusa na pr. s 
različitim gnojenjem, dubinom sjetve, različnim razmacima među 
biljkama, gustoćom sjetve i t. d. jednak broju komparacija. 

Primjer: Treba izvršiti pokus sa gustoćom sjetve kod 
lana. U tu svrhu posijali bismo 5 različnih količina sjemena na 
jedinicu površine, a svaku količinu u pokusu komparirali bismo 
5 puta. Parcelice bismo smjestili u kvadrat i to tako, da je u 
svakom horizontalnom i u svakom vertikalnom pojasu poredano 
bez sistema svih 5 količina sjemena (si. 15). 


SI. 15. 


11 

IV 

V 

III 

I 

III 

II 

IV 

I 

V 

IV 

V 

I 

11 

III 

I 

III 

II 

V 

VI 

V 

I 

III 

IV 

11 


(Rimskim su brojevima označene parcelice s pojedinim količinama sje¬ 
mena po hektaru.) 


6 Tavčar: Biometrika u poljoprivredi 
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Takva se podjela može u praksi upotrebiti onda, ako je 
broj pokusnih članova malen. Tom se podjelom dadu jako eli¬ 
minirati pogreške nastale heterogenošću tla. 

Broj neovisnih varijanata i pogrešku pokusa s gustoćom 
sjetve lana ustanovljujemo ovako: (Tab. XXVI.). 


Tab. XXVI. 



Broj neovisnih 
varijanata 

a) parcelica: na = 25 

b) različitih gustoća sjetve: nb = 5 

c) horizontalnih pojasa: n c = 5 

d) vertikalnih pojasa: na = 5 

na _ i = 25 - 1 = 24 

n b — 1 = 5-1= 4 

nc — 1=5—1= 4 

na — 1 = 5-1= 4 

Razlika = pogreška 

(na — 1) — [(n b — l)-f(nc-l) +(nd-l)] = 12 


b) Primjenu Fisherove metode na nekvadratnoj podjeli par- 
celica prikazat ćemo na slijedećim primjerima. 

Postavljen je sortni pokus od 10 različnih sorata kukuruza 
(koje su označene rimskim brojevima) u 5-erostrukom ponav¬ 
ljanju. Sorte su raspodijeljene u 5 pojasa tako, da je u sva „ k ° 
pojedinom pojasu slobodno smješteno svih 10 sorata. Slučajna 
podjela tih sorata u pojedinim pojasima prikazana je na - 

broj 16. 


SI. 16. 
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Prirodi pojedinih sorata na ovim parcelicama (preračunati 
na hektar — q/ha) bili su slijedeći: (Tab. XXVII.) 
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Analiza varijance našeg-sortnog pokusa je slijedeća: 

(X x) 2 

1. Ukupna (totalna) varijanca: 2(x 2 ) — 

S(x 2 b ) (S x) a 

2. Skupinska (blocks) (pojasna), varijanca: —A-- 

X(x 2 s ) (Xx) ž 

3. Sortna varijanca:--— 

nB °t 

Oznake se odnose na ova svojstva: 

‘ x = prirodi pojedinih parcela preračunati na hektar, 

X B = ukupni prirod jedne skupine (bloka-pojasa), 
n s = broj pokusnih sorata (= 10), 
x s = ukupni prirod jedne sorte sa 5 parcelica 
n T = broj svih pokusnih parcela (= 50), 
n B = broj skupina (blokova) koje kompariramo (— 5). 
Za računanje varijance trebamo još i slijedeće vrijednosti: 
Zbroj kvadrata od: a ) priroda pojedinih parcela preiaču- 
natih na ha, /?) sorata i y) skupina (blokova — pojasa). 

Kvadrati priroda po hektaru s pojedinih parcelica i za po¬ 
jedine sorte sakupljeni su u tabeli XXVIII. 

(- x) 2 

1. Totalna varijanca: 2 (x 2 )-= 


= 39.444,15-^^=137,6778 
(Sx) 2 _ 1 401,9-' 


= 39.306,4722 


---- - 

nx 50 

2. Skupinska (blokna) varijanca: 

S(x 2 b) (Sx) 2 __ 393.356,79 
n s n T 10 


X B 

xu 

286,8 

82254,24 

276,9 

76673,61 

291,5 

84972,29 

270,5 

73170,25 

276,2 

76286,44 


393.356,79 = 2 (X 2 b) 


1 — 39.306,4722 = 

= 29,2068 = skupinska 
varijanca. 


Tab. XXVIII. 
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2(x 4 ) = 39.444,15 3893,06 4396,66 4122,75 3653,87 3972,91 3455,54 3303,91 4371,94 3736,26 4537,25 
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S(x 2 s ) (Sx) 2 197.029,69 

3. Sortna varijanca: ——- ~ 5 

— 39.306,4722 = 99,4658 =-- sortna varijanca. 


X s 

X 2 s 

139,4 

19.432,84 

148,2 

21.963,24 

143,5 

20.592,25 

135,1 

18.252,01 

140,9 

19.852,81 

131,4 

17.269,64 

128,5 

16.512,25 

147,8 

21.844,84 

136,6 

18.659,56 

150,5 

22.650,25 

Z(X's) = 

197.029,69 


/ 


Dakle analiza varijance našeg sortnog pokusa s kukuruzom 
je ovakva: (Tab. XXIX.) 

Pomoću varijance pogreške ustanovit ćemo, kakova je va¬ 
rijanca u pojedinim komponentama. 

Kad bi na pr. sortna varijanca bila približno jednaka vari- 
janci pogreške, onda između sorata ne bi bilo statistički oprav¬ 
danih razlika. 

Isto vrijedi i za varijancu skupina (blokova — pojasa). 

Jesu li spomenute razlike statistički opravdane ili nisu, naj- 
bolie možemo ustanoviti pomoću F i s h e r o v e vrijednos . 
z ii Wl _ 1 Mna 2 2 , gdje su a 2 i i o 2 2 prosjeci kvadrata onih 

komponenata, koje uspoređujemo. 

Prema tome Z je razlika od prirodnih logaritama 
(ln) od obaju o 2 . 

Po Fisherovim tabelama ustanovimo, za koju^je vri¬ 
jednost z = 5o/o (P = 0,05), l«/o (P = 0,01) ih 0,1 / ; ^ °’ 001) ’ 
ako imamo 2 komponente (svojstva) od ^ 1 n 2 individu . 
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Tab. XXIX. 


Neovisnih 

varijanata 


varijance 


Prosječna var. 


a) ukupno 
parcelica 
n r =50 


n-r — 1 = 

= 50 — 1 = 49 


Totalna 

varijanca 

137,6778 


b) skupina bio na — j_ Skupinska var. 29,2067 : 4 = 

kova-pojasa 0,99406 

K = 5- 1 =4 29,2068 = 7,3017 ’ 

ns = 5 


c) sorata 
ns = 10 


n s — 1 = 
= 10-1 = 9 


Sortna var. I 99,4658 : 9 = 


99,4658 


= 11,0517 


1,20129 


Razlika = 
= pogreška 


t 1 a 

o u in 

a 5 « o 
> . o 

5* *s Ć3 

> II 


9,0052 : 36 = 
= 0,25 


- 0,69315 


Za sortnu je varijancu Z= 1,20129 — (—0,69315) = 1,89444. 

Opravdanost ove vrijednosti utvrđujemo po Fishero- 
voj tabeli (Tab. XXX.) za 5% (str. 88). U našem je slučaju 
n i = 4, jer ima (5 — 1 sorta) i n 2 = 36 (razlika = pogreška), 
dakle z = 0,4852, a to znači, da u 5% od svih slučajeva može 
biti vrijednost veća od 0,4852. Prema tome je naša razlika od 
z = 1,89444 opravdana. 

Tek kad ustanovimo, da su razlike između sorata biome- 
trički opravdane, treba izračunati opću pogrešku za prirod 
sorata. 
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Raspodjela z za 5% 


Tab. XXX. 
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Ova se pogreška izračuna ovako: 

varijanca jedne parcelice (Tab. XXIX.).0,25 

varijanca jedne sorte u peterostrukoj komparaciji 5X0,25 = 1,25 

pogreška: m = ]/ varijanca jedne sorte = ]/1,25 = 1,5, a u 
postocima: 


m% = 


m X 100 
Sx 
n s 


1,5 X 100 
140,19 


1,69 


Sx 

ns 


1401,9 

10 


140,19 


Pogreška diferencije: 

m D = m % j/lT = 1,69 X 1.4124 = 2,389 

Signifikantna je razlika u prirodu tek kod 2X2,389 = 
= 4,778%). 


Ako potonju vrijednost upotrebimo kao mjerilo razlike 
između pojedinih sorata i ako uzmemo prirod sorte I kao 100, 
onda su prirodi ostalih 9 pokusnih sorata ovi: (Tab. XXXI.) 


Tab. XXXI. 


Sorta 

apsolutni prirodi 
(Tab. XXVII.) 

Relativni prirod sjemena 

1 

27,88 

100,00 

II 

29,64 

106,31 

III 

28,70 

102,94 

IV 

27,02 

96,91 

V 

28,18 

101,07 

VI 

26,28 

94,26 

VII 

25,70 

92,18 

VIII 

29,56 

106,02 

IX 

27,32 

97,99 

X 

30,10 

107,96 


Signifikantnu razliku od 4,778% obzirom na sortu I imaju 
prema tome samo sorte II, VIII i X, a sorte IV, VI i VII izrazito 
su loše. 
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c) Kvadratna podjela pokusnih parcela 

Želimo ustanoviti, kako utječu na prinos sjemena soje tri 
različne vrste bakterijskog cjepiva. U tu svrhu postavit ćemo 
pokus u četverostrukoj komparaciji. Na kontrolnim parcelama 
(K) posijat ćemo necijepljeno sjeme. 

Parcelice zasijane cijepljenim sjemenom označene su s. A, 
B i C; one su nesistematski poredane u obliku kvadrata (latin 

square): (SI. 17). 


SI. 17. 


K 

B 

C 

A 

A 

K 

B 

C 

c 

A 

K 

B 

B 

C 

A 

K 


Prirod sjemena na pojedinim parcelicama, preračunat na 
hektar, prikazan je u tabeli XXXII. 


Tab. XXXII. 




Xft 

X 2 h 


K 

18.73 

A 

18,93 

C 

19,16 

B 

19,42 

76,24 

5.812,5376 


B 

19,68 

C 

18,95 

K 

17,26 

A 

18,51 

74,40 

5.535,3600 


A 

18,86 

K 

16,38 

B 

20,03 

C 

19,28 

74,50 

5.557,7025 


C 

18,76 

B 

19,35 

A 

19,32 

K. 

17,84 

75.27 

5.665,5729 

Xv 

76,03 

73,61 

75,77 

75,05 

2x h = 300,46 

£ X 2 h = 
22.571,173 

X 2 v 

5.780,5609 

5.418,4321 

5.741,0929 

5.632,5025 




| £x\ = 22.572,5884 




Za izračunavanje varijance i za analizu varijance treba 
odrediti slijedeće vrijednosti: 

a) sumu priroda u sjemenu u svakom horizontalnom redu 

(h), 

b) sumu svih horizontalnih redova, 

c) kvadrate ukupne sume priroda u svakom horizontalnom 
redu zasebno, 

d) sumu tih kvadrata, 

e) sumu priroda u svakom vertikalnom redu (v) zasebno, 

f) kvadrate tih zbrojeva, 

g) sumu kvadrata priroda pojedinih parcelica (tab. XXXIII). 


Tab. XXXIII. KVADRATI: 



K 

350,8129 

A 

358,3449 

C 

367,1056 

B 

377,1364 


B 

387,3024 

C 

359,1025 

K 

297,9076 

A 

342,6201 

A 

355,6996 

K 

268,3044 

B 

401,2009 

C 

371,7184 

C 

351,9376 

B 

374,4225 

A 

373,2624 

K 

318,2656 

X 2 

1.445,7525 

1.360,1743 

1.439,4765 

1.409,7405 = 

X 2 


£(x 2 t) = 5.655,1438 


Nadalje moramo sakupiti u posebnu tabelu prirode svake 
cijepljene i necijepljene parcelice (Tab. XXXIV.) 

Tab. XXXIV. 



K 

A 

B 

C 



18,73 

18,93 

19,42 

19,16 



17,26 

18,51 

19,63 

18,95 



16,38 

18,86 

20,03 

19,28 



17,84 

19,32 

19,35 

18,76 


X 

70,21 

75,62 

78,48 

76,15 

£ x = 300,46 

Prosjek 1 

17,5525 

18,9050 

19,6200 

19,0375 

£1=75,115 

X2 

4.929,4441 

5.718,3844 

6.159,1104 

5.798,8225 

V 

221505/7614 


T=10° 

4 

107,7054 

111,7789 

108,4603 
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n x = 16 (pokusnih parcelica) 
n v = 4 (broj vertikalnih pojasa) 
n h = 4 (broj horizontalnih pojasa) 

Nakon toga sve te sume zbrojimo. 

Iz navedenih podataka izračunamo pojedine varijance po 
ovim formulama: 


Ukupna ili totalna varijanca (Tab. XXXIII. i 
XXXII.): 

2 (x 2 t ) — = 5.655,1431 — = 12,8831 

(S x) 2 = 300.4 6 2 = 5 ^2 26 (Tab. XXXII.) 
n T 16 v 

Varijanca horizontalnih redova (h): (Tab. 
XXXII) 

_ lSx)i = 22.571J73 _ ^ = 0 53 

n v n T 4 

Varijanca vertikalnih redova (v): Tab. (XXXII.) 


S(x 2 v ) (2 x) 2 
n h n T 


22.572,5884 

4 


5.642,26 = 0,8871 


Varijanca pokusa (cijepljena i necijeplje- 
n a soja): (Tab. XXXIV. i XXXII.) 


S ( X V _ = ?2 : 605Jbl4 _ 5 g 4 2 26 = 9,1803 

n rep . n T 4 


Nadalje treba izračunati varijancu pogreške i vrijednost Z. 
(vidi tabelu XXXV.). 


Da li je ova razlika statistički opravdana, ustanovljujemo 
po Fisherovoj tabeli (Tab. XXX.) za Z od 5%. Kod nas je ni—3, 
jer je za pokus (4—1) = 3 neovisnih varijanata, a za n 2 = 6, a 
to izlazi iz neovisnih varijanata za razliku = pogrešku. Prema 
tome je Z = 0,7798; to znači, da samo u 5% može biti ta vrijed¬ 
nost slučajno prekoračena. 


Kako je u našem slučaju Z = 1,05507, to je ta razlika bio- 
metrički posve opravdana. Dakle ima opravdanih razlika u pri¬ 
nosu sjemena između necijepljene soje i soje cijepljene različnim 
bakterijskim cjepivima. 


93 





94 


Pogrešku diferencije izračunat ćemo na slijedeći način: 


prosječna varijanca jedne parcelice (pogreške) 
varijanca jednog cijepljenja 

= 0,37095 

ili necijepljenja (4 komparacije) 4 X 0,37095 

=--1,4838 

pogreška: m = |/varijanca = |/1,4838 

= 1,2183 

1,2183X100 
m /“ — 75,115 

= 1,6221 

Pogreška diferencije =m]/2 = 1,6221 X 1,4142 

= 2,29397 

Dakle je signifikantna razlika u prinosu 
2.2,29397 = 4,58794%. 

sjemena kod 


Kad s tim mjerilom (4,58794%) usporedimo prirode neci- 
jepljene soje (K) = 100 s cijepljenim A, B i C izlazi: 

(Tab. XXXVI.) 


Tab. XXXVI. 


Soja 

Relativni 

prirodi 

Apsolutni 

prirodi 

K 

100 

17,5525 

A 

107,71 

18,9050 

B 

111,78 

19,6200 

C 

108,46 

19,0375 


Vidimo, da su sa svim cjepivima postignuti opravdano bolji 
prirodi u poredbi s prinosom necijepljene soje. 


' B — 2. BIOMETRIKA KORELACIJE 

' 

J Općenito 

Do sada smo obrađivali varijaciono-statistički svako poje¬ 
dino svojstvo organizma zasebno. Istraživajući različna svoj¬ 
stva jednog istog organizma opažamo, da neka od njih uvijek 
dolaze u zajednici. Kažemo, da između takvih svojstava ima 
neka ovisnost ili korelacija. 

Na pr. brašnasto zrno pšenice sadržava manje bjelančevina 
nego caklavo. Zrno zobi s mnogo masti lakše je od zrna s manje 
masti. Zrno ječma s više bjelančevina teže je od zrna s manjom 
količinom bjelančevina. Sjeme soje sa više masti sadržava ma¬ 
nje bjelančevina nego zrno s manje masti i t. d. 

Za oplemenjivanje gospodarskih biljaka i životinja vrlo je 
važno da doznamo, koja su svojstva organizma međusobno 
ovisna, a koja su neovisna. Kad to doznamo, treba nam proučiti 
samo ono svojstvo, koje se dade lakše i brže istražiti, a onda 
po njemu zaključujemo, kakva je vrijednost ostalih svojstava. 

Po Johannsenu su organizmi »sistemi u dinamičkoj ravno¬ 
teži«. Razlikujemo: a) fiziološke ko r e 1 a c i j e t. j. fizio¬ 
loške povezanosti pojedinih organa u organizmu kao i među- 
I sobni utjecaj gena, koji putem spolnih stanica (gameta) dolaze 

u novi organizam, b) korelativni varijabilitet t. j. 

I varijabilitet u povezanosti svojstava jednog organizma u po- 

i redbi s drugim organizmom. 

Ovdje nas specijalno interesira korelativni varija- 
b i 1 i t e t t. j. slučajevi, gdje treba ustanoviti, da li pojedina 
svojstva organizma variraju ovisno ili neovisno o drugim svoj¬ 
stvima. 

Konstantnu vezu dviju veličina zovemo u matematici funkci¬ 
jom. U matematičkoj funkciji možemo prvu veličinu izračunati 
iz druge, i obratno. Određenoj nekoj vrijednosti jedne veličine 
odgovara samo jedna određena vrijednost druge veličine. Na pr. 
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ako je y = 4, i ako je x = 3, onda je xy = 12. Ovdje su po¬ 
znate uzročne veze. 

Kod korelativnog varijabiliteta osim poznatih uzroka ima 
i nepoznatih tako, da se iz jednog svojstva jedne varijante ne 
može točno unaprijed odrediti i drugo njezino svojstvo. Dakle, 
jedno će se svojstvo izvjesnog organizma razviti do izvjesnih 
granica različito, ako se razvija pod drugim intenzitetom vanj¬ 
skih faktora. Prema tome ono varira. 

Što je više poznatih činilaca, koji djeluju na istraživana 
svojstva, to se točnije može unaprijed odrediti veza među tim 
svojstvima. Istražujemo li samo jednu ili vrlo mali broj vari¬ 
janata (organizama), onda radi event. jačeg utjecaja kojeg ne¬ 
poznatog činioca na promjenu nekog svojstva ne možemo naći 
povezanost s kojim drugim svojstvom, makar ona i postojala. 
Tek kod velikog broja varijanata može se točnije odrediti, jesu 
li neka svojstva organizma povezana (u korelaciji) ili nisu. 

Naučenjaci su nastojali da izrade matematičku metodu, ko¬ 
jom bi se dao mjeriti stupanj korelacije. 

Jednu je takvu metodu izradio poznati biolog G a 11 o n. 
Kasnije su ovu metodu usavršili i neki drugi naučenjaci. 

Gdje su dva svojstva međusobno linearno povezana, istra¬ 
žujemo jačinu korelacije najčešće po Bravaisovoj me¬ 
todi. Ovom se metodom izračunava t. zv. korelacioni 
koeficijent, a označujemo ga sa r. 

U pređašnjim poglavljima izračunali smo samo neke vari- 
jaciono-statističke vrijednosti za duljinu klasa i broj klasića kod 
pšenice (tab. III. i tab. VII.). Zanima nas sada na pr., da li su 
duljina klasa i broj klasića na klasu međusobno povezana svoj¬ 
stva, a ako jesu, kakva je jačina njihove povezanosti. 

Da se pronađe event. zavisnost između tih dvaju svojstava, 
treba najprije izraditi t. zv. korelacion u-tabelu. U toj 
se tabeli svaka varijanta smjesti u onaj varijacioni red, koji joj 
pripada po njezinoj duljini klasa i broju klasića. 

Jedno svojstvo t. zv. suponirano označit ćemo sa x, a 
drugo »relativno« sa y. Posve je svejedno, koje ćemo svoj¬ 
stvo uzeti za suponirano, a koje za relativno. Ako su dva svoj¬ 
stva u međusobnoj ovisnosti (korelaciji), onda varijanta V 
obzirom na svojstvo x može varirati u izvjesnim granicama, koje 
određuje svojstvo y. 


Raspored svih varijanataV x za svojstvo x zavisi o svojstvu 

x KoTeTf SVih Vadj f ata V y za svo J' s tvo y : zavisi o svojstvu 
, , 0 r a c 1 o n a t a b e 1 a za duljinu klasa i broj klasića na 

klasu je ovakova: (Tab. XXXVII). 

Za različite parove svojstava korelacione su tabele razli- 

tome R tT 1 frekvenci j a varijanata u tim tabelama zavisi o 
tome, kako su jako međusobno povezana svojstva, koja sači¬ 
njavaju neki istraživani par. 

Općeniti oblici takvih tabela mogu biti i ovakvi: 1 
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a ^ 

SI. 18. f 

. ; t ■ 

S1 - br - 18 a ) pozitivna korelacija, b) negativna korelacija i 
c) korelacije nema. J 

v P [ e i ma l . t0me zay isnost ili korelacija između dva svojstva 
može biti ili pozitivna (si. a) ili negativna (si. b). 

Ako su svojstva međusobno nezavisna, onda su varijante 
poredane po čitavoj korelacionoj tabeli tako, da posve zapre¬ 
maju površinu kruga, kako se vidi na si. pod c) 

karJt rStU k ° relaC J je ( + ) ’ (-) iH (0 > 1 °j ezinu jačinu najbolje 
produkt razmaka pojedinih varijanata od srednje 
rijednosti M x za suponovano i M y za relativno svojstvo. 

Ako u našoj korelacionoj tabeli za svojstvo klasa pšenice 

SZ Tm l Crt “ “ raZredu ' “ « nalazi srednja 

vrijednost (M x t za suponovano svojstvo (x), i horizon- 

M.u crtu u razredu u kojemu se nalazi srednja vrijednost 

M y )^ za svojstvo (y), razdijeli se korelaciona tabela na četiri 

' gdj ' 86 ° Ve crte križ “j». označit čemo sa 0. 

(Slika 18). 

Na horizontalnoj su crti desno od 0 pozitivne, a lijevo ne¬ 
gativne varijante. Na vertikalnoj, crti pozitivne su varijante ispod 

7 Tavčai*: Biometrika a poljoprivredi 
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horizontalne linije a negativne su iznad nje. 

Tako nastaju pozitivni i negativni razmaci pojedinih vari¬ 
janata od srednjih vrijednosti. Produkti tih razmaka su ili po¬ 
zitivni ili negativni. 

Ako spomenute kvadrante označimo s I, II, III, i IV (vidi 
korelacionu tab. XXXVII. i sliku br. 18), onda su produkti raz¬ 
maka pojedinih varijanata 2pa x a y obzirom na srednju vrijed¬ 
nost za I. i IV. kvadrant pozitivni, a za II. i III. kvadrant su 
negativni. 

Ako je suma spomenutih produkata iz I. i IV. kvadranta 
veća od sume iz II. i III. kvadranta, onda je korelacija pozitivna 

^P a x a y = + t. j. ako se jedno svojstvo promijeni u jednom 
smjeru, mijenja se i drugo svojstvo u istom smjeru. 

Kad je suma spomenutih produkata iz II. i III. kvadranta 
veća od sume iz I. i IV., onda je korelacija negativna t. j, 
2 P a x a y = (—); to znači: kad se jedno svojstvo mijenja u jed¬ 
nom smjeru, drugo se svojstvo mijenja u protivnom smjeru. 

U onom slučaju, kad su varijante u jednakom broju i jedna¬ 
kom poretku porazmještene u svim kvadrantima, suma je pozi¬ 
tivnih produkata jednaka sumi negativnih produkata, dakle 
^P*x a y = 0 i istraživana svojstva nisu međusobno povezana. 

Korelacija je to jača, što je razlika između sume produ¬ 
kata iz I + IV i II + III kvadranta veća. Jačina pak zavisnosti 
jednog svojstva od drugog ne zavisi samo o veličini spomenute 
razlike, već također i o broju svih varijanata (n) i o varijabilitetu 
( a x > a y) suponovanog i relativnog svojstva. Dakle, pravi se stu¬ 
panj korelacije može izračunati tako, da se razlika spomenutih 
produkata iz I + IV i iz II + III kvadranta podijeli s produ¬ 
ktom dobivenim od broja svih varijanata (n) i standardne devi¬ 
jacije a x suponovanog i a y od relativnog svojstva. 

LINEARNA KORELACIJA 

L Bravaisov korelacioni koeficijent za 
mnogo varijanata 

Linearna korelacija je ona, kod koje se u korelacionoj ta¬ 
beli varijante raspodjeljuju diagonalno u obliku elipse. 

Linearnu korelaciju izrazujemo pomoću Bravaisovog 
korelacionog koeficijenta tako, da zbroj produkata svake po- 
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jedine frekvencije i njezine udaljenosti od srednje vrijednosti x 
svojstva i srednje vrijednosti y svojstva, podijelimo s produktom 
broja svih varijanata i standardne devijacije x i y svojstva. 
Dakle općenita formula za Bravaisov korelacioni koeficijent 
glasi: 

r ^ S p a x a y 

II G x Oy 

Kad je korelacija između dva svojstva potpuna t. j. ako je 

a„ = a i a = a , onda je r = =1. U toj jednadžbi 

’ y no" 

je = ° 2 (vidi str. 30). 

Iz toga izlazi, da se vrijednost korelacionog koeficijenta 
kod pozitivne korelacije kreće između 0 i 1, a kod negativne 
između 0 i —1. 

Kod računske se operacije ne služimo razmacima od srednje 
vrijednosti M, nego od aproksimativne srednje vrijednosti A, 
s kojom smo se već upoznali (vidi str. 22). S A x označit ćemo 
aproksimativnu srednju vrijednost svojstva x, a s A y svojstva y. 

Iz toga izvodimo' slijedeću formulu, kojom se služimo kod 

izračunavanja Bravaisovog korelacionog koefi- 

c i j e n t a, a koja glasi: 

S pa x 3y nb x by i* 

r — -—--5-i- U toi le formuli: 

np x Py 


p .frekvencija 

a.areal 

n . ..ukupni broj varijanata 

b.(vidi str. 22) 


Ova formula vrijedi samo za linearne korelacije t. j. kad 
su varijante poredane tako, kao što je prikazano u našoj ko- 
relacionoj tabeli br. XXXVII. ili na slici br. 18 a i b. Ako kore¬ 
lacija nije linearna nego na pr. parabolična t. j. kad su vari¬ 
jante poredane kao u tab. XLV., str. 121, onda je računanje 
korelacije veoma komplicirano. Za poljoprivrednu praksu posve 
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Kako se to radi, vidjet ćettro na 'Sje^ 

' a ) Pozitivna korelacija: 

dulii^ "i?” 1 Primj f U sa P šeni com želimo ustanoviti, da li je 
w2^fliSfc a Ut,r 3u Il SObnoj P° vezanosti » brojem klasića na 

Trelat l P ° ređamo vari J a nte za suponovano (x) 
i relativno (y) svojstvo u korelacionu tabelu. ' 

Iz podjele se varijanata u korelacionoj tabeli vidi da broi 

SSS Pretr,' đuljina - k ' asa sloie “ 

vaisovoj fZuH ' ‘ 2raCUnat Ćem ° jafl ”“ korelac ti e PO Bra- 

y —— pa x 3y — nb x by 

flpx Py 

jednosti:° j ^ f ° rmUli P ° Znate Već 0tprije sli jedeće vri- 

• / ** = -0,4475; b y = _ 0,4775; n = 400 

^ ■ ' Px= 2,934; p y = 3 ; 398 

Redovi u smjeru od A x i Ay dijele korelacionu tabelu na 
cetin kvadranta. U ovim je redovima 2pa x a y = 0, jer su a x ili 

% ih i a x . a y jednaki 0. Produkti a x a y su u I. i IV. kvadrantu 
pozitivni, a u II. i III. negativni. 

U našem slučaju je u IV. kvadrantu i to u najdonjem redu, 

gdje je frekvencija p = 2, produkt pa x a y = 112 jer je a - 7 
i a y = 8 . . j j a x — / 

a, mn ° Žim ° (P) S 


I- kvadrant 
(pozitivni produkti) 

P a y 

18 . 1 . 1 18 

11.1.2 22 

5.1.3 15 

2.1.4 8 

1.1.5 5 

I’ 4 * 2 . 1 28 

18 .2 . 2 72 

9.2.3 54 


P 

7 

6 

6 

2 

5 

3 

6 

4 


IV. kvadrant 
(pozitivni produkti) 

• 1.1 7 

• 1.2 12 

• 1.3 18 

1.4 8 

- 2.1 io 

• 2.2 12 

■ 2.3 36 

• 2.4 32 
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Ukupno 1553 


II. kvadrant 


(negativni produkti) 

p a x • 3y 


2 .2.2 

7 

1.3.1 

8 

4.1.2 

6 

7.1.1 

8 

3 .2.1 

3 


Ukupno — 32 


III. kvadrant 
(negativni produkti) 


P 

• a x 

a y 


5 

. i . 

1 

5 

3 

. i . 

2 

6 

3 

2 

1 

6 

1 

. 2 . 

2 

4 

1 

. 3 - 

1 

3 


Ukupno 


— 24 
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2p a x a y = (2062 + 1553) - (24 + 32) = 3559; 

2 p a x a y - n b x b v 3559 - (400 • - 0,4475 • - 0,4775) 

n pxpy — 400 • 2,934 • 3^98 —0,871 

Prema tome je korelacioni koeficijent za duljinu klasa i 
broj klasića na klasu r = 0,871 a to znači, da su oba svojstva 
u vrlo jakoj međusobnoj vezi t. j. između njih postoji vrlo jaka 
pozitivna korelacija. 

Iz korelacionog koeficijenta (r) možemo odrediti jačinu 
korelacije poRoemerOrphalovoj tabeli: (Tab. XXXVIII.) 


Tab. XXXVIII. 


Korelacioni 
koeficijent :r 

Korelacija 

0,0 -01 

nema | 

0,1 —0,25 

jako slaba 

0,25 - 0,4 

slaba 

O 

1 

O 

Ul 

f srednja j 

0,5 -0,75 

jaka 

0,75 - 0,9 

vrla jaka 

0,9 -1 

potpuna 


Evo još jednog primjera iz živinogojstva za pozitivnu ko¬ 
relaciju: 

U korelacionoj tabeli XXXIX. na str. 104 prikazana je vari¬ 
jacija za nosivost i debljinu niti vune kod solčavske ovce. Htjeli 
bismo znati, kakva je veza između nosivosti i debljine niti vune; 
važno je, da to znaju stočari, koji oplemenjuju ovce na kvali¬ 
tetu vune. 



Ova je svojstva istražio prof. dr. Ogrizek. U korelacionoj 
tabeli je xs vojstvo nosivost n iti u gramima, a y svojstvo dej- 
bljin a niti u mikronima. 


liOŽr'-yi 




J ^ U ISIS 

T /l{j 



I 
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Tab. XXXIX. 


1 Nosivost 
niti u gr 
(x) 23- 

1 

Debljina niti u ji (y) 

25 26-35 29- 
1 

45 32- 

55 35' 

65 38-75 41-85 44‘95 48- 

lili 

05 

P 

10-5 

15-5 

20-5 

25-5 

30-5 

35-5 

40-5 

45-5 

P 

1 








1 | 

m 

i 







1 

u 

6 

9 

24 

1 

3 

4 


47 

m 


1 

30 

21 

3 

1 


56 

■ 




43 

10 



53 




4 

1 

9 

5 

1 

19 







8 

7 

10 








7 

7 

i 

7 

10 

58 

66 

25 

18 

15 

200 


Kako se iz korelacione tabele vidi, i ovdje je korelacija 
linearna, tako da možemo izračunati korelacioni koeficijent po 
Bravaisovoj formuli. Po toj formuli naći ćemo, da je r=0,779, 
a to znači, da je između debljine i nosivosti niti vrlo jaka po¬ 
zitivna korelacija. 


b) Negativna korelacija 

Istraživali smo, u kakvoj je vezi stupanj caklavosti 
zrna pšenice s duljinom vremena potrebnog za meljavu. U tu 
smo svrhu proučili 360 različnih sorata pšenice. 

U korelacionoj tabeli br. XL. prikazan je raspored istraže¬ 
nih varijanata-sorata. Iz te se tabele vidi, da i ovdje postoji 
linearna, ali negativna korelacija. 

Korelacioni koeficijent izračunan po Bravaisovoj for¬ 
muli je r = — 0,736 t. j. caklavost zrna i trajanje meljave u 
jakoj su negativnoj zavisnosti. 



♦ 


Prosječno trajanje meljave — u sekundama (") 
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2 . Izračunavanje korelacionog koeficijenta 
za mali broj varijanata 


Kad je broj varijanata, za koje treba izračunati korelacioni 
koeficijent, razmjerno malen, ne možemo se služiti korelacio- 
nom tabelom. U tom slučaju izračunavamo eventualnu ovisnost 
između dva svojstva pomoću ove formule za korelacioni koe¬ 
ficijent: 


r _ SV x V y — n M x M y 
n o x Oy 

Primjenu ove formule prikazat ćemo na slijedećem pri¬ 
mjeru: 

12 različnih sorata soje istražili smo sadržinu sirovih 
bjelančevina (sirovi protein) i sirove masti u sjemenu. Slijede¬ 
ćom i acunskom operacijom nastojali smo pronaći, da li postoji 
neka zavisnost među spomenutim tvarima. (Tab. XLI.). 

Sa x označit ćemo sirovu mast a sa y sirovi protein. Prema 
tome su V x varijante sirove masu, a V y varijante sirovog pro¬ 
teina. M x je srednja vrijednost za sirovu mast, M y za sirovi 
protein, n je ukupni broj varijanata, a o x i o y su odgovarajuće 
standardne devijacije. 


Pojedine vrijednosti korelacionog koeficijenta izračunat 
ćemo iz podataka u Tab. XLI. , 


2V Xr= 224 1 19 _ _2V y ... 45771 

12 WUU5, M 






38-1425; V 

"4225,^891 
/ 12 


18-6825 = 1-7546 


38-1425 = 2-2472 


a 1/2 V 2 y ~\[ 17518-8063 “ 

7 y = j ~n ~ M '>- = /- 12 -38-1425 = 2-2472 

r = 2v xV y -nM x M v 8511-5882 — 12 X 18-6825 X 38-1425 

n °x a T 12 X 1 -7546 X 2 2472 — 

= — 0-8004 

Kako je r = —0,8004, znači, da su u našem slučaju koli¬ 
čina sirove masti i sirovog proteina u sjemenu soje u vrlo jakoj 
negativnoj korelaciji. 
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Tab. XLI. 



3. Srednja i vjerojatna pogreška 
korelacionog koeficijenta 



Da li je korelacioni koeficijent, izračunat na temelju neke 
korelacione tabele, varijaciono-statistički opravdan, najbolje 
nam pokazuje njegova srednja odnosno vjerojatna pogreška. 
Srednju pogrešku korelacionog koeficijenta 
izračunatćemo po Pearso n-F ilonovoj foi muli. 

mr — 1 , a vjerojatnu pogrešku po formuli 

l/n 

E P. r = 0,6745 • gdje je n broj varijanata, a r korela- 

V n ... 

cioni koeficijent. Pokusima je ustanovljeno, da se iz korelacio¬ 
nog koeficijenta može zaključiti na korelaciju između svojstava 
samo onda, ako je on barem tri puta veći od svoje srednje 
odnosno vjerojatne pogreške. 

Izračunali smo, da je korelacioni koeficijent za duljinu 
klasa i broj klasića r = 0,871. Prema tome je srednja po- 
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1 _o 0712 

g r e s k a tog koeficijenta: m r = — = 0,012, a vjero- 

V400 J 

1 _n 071 2 

jatna pogreška E. P r = 0,6745 . - =0,0081 

J/400 

Dakle: r ± m r .... 0,871 ± 0,012 i 
r±E. P r . . . 0,871 ±0,0081 


Vidimo, da je trostruka srednja, odnosno vjerojatna 
pogreška manja od korelacionog koeficijenta, a to znači, da je 
korelacija duljine klasa i broja klasića varijaciono-statistički 
opravdana. 


Ako je broj varijanata malen, srednja i vjerojatna pogreška 
korelacionog koeficijenta je točnija, ako se mjesto divizora 
]/ n služimo J/n — 1 


4. Regresija 

Za određivanje intenziteta, kojim jedno svojstvo reagira na 
promjene nekog drugog svojstva, služi regresija t. j. ko- 
relacioni koeficijent (r) pomnožen kvocijentom standardne de¬ 
vijacije X i Y svojstva. 

Rx/y — r — odnosno R y/X = r • — 

°y a x 

U našem primjeru za pšenicu mogli bismo izračunati re¬ 
gresiju na ovaj način: 

a) Da doznamo, koliko se produlji ili skrati klas (relativno 
ili y svojstvo), kad se broj klasića (suponovano ili x svojstvo) 
poveća ili smanji za jedinicu moramo upotrebiti formulu: 

Ry/x = r ^ =0,871 • 0,504 cm, jer je a x = 2,934, a o y = 

= 1,699 (vidi str. 32). 

Dakle, ako se broj klasića na klasu poveća ili smanji za 1, 
klas se produlji, odnosno skrati za 0,504 cm. 

b) Da doznamo, za koliko se poveća ili smanji broj klasića 
na klasu (suponovano ili x svojstvo), ako se duljina klasa (rela¬ 
tivno ili y svojstvo) promijeni za 1 cm, onda se služimo for¬ 
mulom: 


R . == r — =0,871 • ^’cnn '~ 1,503, a to znači kad se broj klasića 
x/y a y l,b99 

na klasu promijeni za 1,503, onda se duljina klasa mijenja za 1 cm. 

5. Srednja i vjerojatna pogreška regresije 

Srednju pogrešku regresije izračunamo po formuli: 

m R = m r ako je R x/y i 
°y 

po formuli: 

m R = m r • —, ako je R y/X . 
o x 

Vjerojatna pogreška regresije je E. P. R = 0,6745 m R . 

6. Korelacija između više svojstava — 
parcijalna korelacija 

Proučili smo međusobnu ovisnost dvaju svojstava bez obzira 
na druga svojstva, s kojima ova dva svojstva koreliraju. Ima 
slučajeva, kad treba istražiti korelaciju dvaju svojstava elimi¬ 
niravši utjecaj trećeg ili četvrtog i t. d. svojstva, koja su u 
međusobnoj vezi s prva dva svojstva. Na pr. biometrički smo 
istraživali duljinu, širinu i debljinu sjemena gra¬ 
ha i po Bravaisovoj formuli izračunali korelacioni koeficijent 
za: a) duljinu i širinu, b) duljinu i debljinu i c) širinu i debljinu. 
Ovdje nam korelacioni koeficijent izrazuje stupanj korelacije 
između dva svojstva uz sudjelovanje i trećeg, dakle korelacioni 
koeficijent za duljinu i širinu označuje jačinu povezanosti tih 
dvaju svojstava uz sudjelovanje i trećeg svojstva t. j. debljine. 
Isto tako korelacioni koeficijent za duljinu i debljinu označuje 
međusobnu ovisnost tih dvaju svojstava uz sudjelovanje i širine, 
a korelacioni koeficijent za širinu i debljinu mjeri povezanost 
ovih dvaju svojstava uz sudjelovanje duljine. 

Kad treba izračunati korelacioni koeficijent za dva svoj¬ 
stva eliminirajući pri tom treće svojstvo, služimo se parcijal¬ 
nim korelacionim koeficijentima, a izračunavamo ih po ovim 
formulama: 
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II. r ,3 


(1 — r 12 2 ) (1 — r 23 2 ) ’ 


III r _ *23 1 12 M3 _ 

‘ 23/1 l/(l-r,?)(l-r^> 

r i 2/3 j e parcijalni korelacioni koeficijent za svojstvo 1 
i 2, kad je eliminiran utjecaj 3. svojstva; 

r, 3 / 2 znači korelacioni koeficijent za svojstva 1 i 3, ako 
je eliminiran utjecaj svojstva 2 i 

r 23 /, je korelacioni koeficijent za svojstva 2 i 3 elimini¬ 
ravši utjecaj svojstva 1 , 

r 12 je korelacioni koeficijent izračunat po Bravaisovoj 
formuli 


S pa x - a y — nb x b y 

n Px py 

kad i svojstvo 3 sudjeluje u razvoju svojstva 1 i 2, ri 3 je kore¬ 
lacioni koeficijent za svojstvo 1 i 3 uz sudjelovanje svojstva 2, 
u I 23 je korelacioni koeficijent za svojstva 2 i 3, kad i svojstvo 
1 utječe na razvoj svojstva 2 i 3. 

Primjenu tih formula najbolje ćemo vidjeti na slijedećem 
primjeru s grahom: 

Označimo s 1 duljinu, s 2 širinu i s 3 debljinu zrna graha. 
Ako po spomenutoj Bravaisovoj formuli izračunamo korelacione 
koeficijente za pojedine parove dimenzija, naći ćemo za njih 
slijedeće vrijednosti: r i2 — 0,733; ri 3 = 0,691; r 2 3 = 0,628. 


Parcijalni korelacioni koeficijent između duljine (1) i širine 
(2), ako je debljina (3) sjemena konstantna, izračuna se po for¬ 
muli I.: 


r _ r i 2 — r ia • r 23 - _ 0,733-(0,691 X0,628) _ n 

12/3 ]/(l — r 13 2 ) (1 — r 28 2 )" 1/(1 - 0,691 2 ) (1 -0,628 a ) 

Iz toga slijedi, da je parcijalni korelacioni koeficijent za 
•" 12/3 manji od Bravaisova ri 2 . 

Po formuli II izračuna se parcijalni korelacioni koeficijent 
P13/2 = 0,435 za duljinu (1) i debljinu (3), a po formuli III par¬ 
cijalni korelacioni koeficijent r 23 /x = 0,258 za širinu (2) i de¬ 
bljinu (3). 

Ako kod četiri svojstva u korelaciji želimo istražiti jačinu 
zavisnosti između jednog relativnog i jednog od triju supono- 
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vanih svojstava, uz pretpostavku, da su druga dva suponirana 
svojstva konstantna, onda izračunamo parcijalni korelacioni 
koeficijent po slijedećim formulama: 

I\ _ f 12/4 r !3/4 ' r 23/4 

odnosno po formulama: 

jt\ _ r i3/4 r !2'4 ’ r 23/4 

13/24 1/(1 -r 2 12 / 4 )(l-r 2 23 / + ) 

im -_ r 23/4 r i2 ; 4 ‘ f 13/4 

T\/\ r _ f 24/8 r i 2/3 ' r i4/3 


\n r _ r 34 / 2 r !3/2 * f 14/s 

84/12 1/(1— r2 13 / 2 )(l — r 2 14/2 ) 

\T|\ -_ r i4/3 r iS/3 ‘ r 24/3 


NELINEARNA KORELACIJA 

Osim linearnih korelacija ima i različnih nelinearnih. Ne¬ 
linearne korelacije pojavljuju se, kad u korelacionoj tabeli va¬ 
rijante nisu poredane ni u elipsi, ni u krugu nego u kojem dru¬ 
gom obliku. 

Česta nelinearna korelacija kod poljoprivrednih biometrič- 
kih istraživanja je parabolična korelacija. Takva je korelacija 
na pr. između caklavosti zrna pšenice i konzistencije tijesta do¬ 
bivenog iz prekrupe (Pelshenkeova metoda). 

Konzistenciju tijesta ustanovljujemo na taj način, da 5 
prekrupe pšeničnog zrna zamijesimo s 0,25 g kvasca i 2,5 do 
ccm vode od 32° C. Od toga tijesta načinimo kuglicu i ubacimo 
je u staklenu čašu, koja je do 3 A napunjena vodom tempe¬ 
rature 31° do 33° C. Vrijeme, otkad kuglica u čaši ispliva na 
površinu vode, pa do časa, kada se raspane, uzima se kao mje¬ 
rilo konzistencije tijesta. Što se dulje kuglica u vodi raspada, 
to je bolja i konzistencija. 


i 


bJD co 
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U korelađonoj tabeli br. XLV. str. 121 prikazana je podjela 
360 sorata pšenice istraženih na caklavost i konzistenciju tije¬ 
sta. Iz te se tabele vidi, da su spomenuta svojstva od 0 pa do 
75% caklavosti u pozitivnoj, a iznad 75% u negativnoj korela¬ 
ciji. Korelacija ovdje nije linearna, pa stoga ne možemo upo¬ 
trijebiti formulu za Bravaisov korelacioni koeficijent. Točniju 
sliku raspodjele varijanata i njihove korelacije možemo u ovom 
slučaju dobiti grafičkim prikazom. 

7. Grafički prikaz korelacije 

Radi bolje ocjene treba katkada korelaciju 1 prikazati i gra¬ 
fički pomoću regresivnih linija. U tu svrhu treba naj¬ 
prije izračunati za relativno i suponovano svpjstvo prosječne 
razredne vrijednosti, i to za svaki razred napose. Na temelju 
tih brojaka konstruiraju se faktične regresione linije. Teoretske 
regresione linije za linearne i parabolične korelacije ižračunat 
ćemo najbolje pomoću metode najmanjih kvadrata uz uporabu 
Mikičevihformulai tabelaitoza: (str. 113). 

a) linearnu podjelu varijanata po formuli za pravac . . . . 
y = a + bx. 

Tab. XLII. 


Konstante u linearnoj jednadžbi za: 


j a) parni broj razreda 


p) neparni broj razreda | 

n 

k* 

2 <n) 


n 

k,* 

*<2 1 

2 

1 

3 


i 

0 

1 

4 

10 

10 


3 

2 

2 

6 

35 

21 


5 

10 

3 

8 

84 

36 


7 

28 

4 

10 

165 

55 


9 

60 

5 

12 

286 

78 


11 

110 

6 

14 

455 

105 


13 

182 

7 

16 

680 

136 


15 

280 

8 

18 

969 

171 


17 

408 

9 

20 

1330 

210 


19 

570 

10 

22 

1771 

253 


21 

770 

11 

: 24 

2300 

300 


23 

1012 

12 


* Vrijede i za parabolu 
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a) parni broj razreda 

b = i[2 2(xy) -(n + l) 2y] 

a = {[2(y)-2(n) b] 

/?) neparni broj razreda 

b=i-[Z(xy)-k 2 Sy] 

a — — • 2 v — k Q b 
n J 

b) Paraboličnu podjelu varijanata: y = a + bx + cx 2 
c = jy-[« s (y) — p2(xy) + ?2 (x 2 y)] 

b =-^-[2 2 (xy) — (n + 1) -Z(y)] — (n -f 1) c— (parni broj vari¬ 
janata) 

b = ~ [2 (xy) — k 2 2 y] — (n-fl)-c-(neparni broj vari- 

K 1 

janata) 

a = i-[2(y)-S (n) b — 2 (n 2 ) ■ c] 

U ovim su formulama: 

n ... broj razreda; vrijednosti n, k, kj, k 2 , k 3 a, ji i y nađemo 
u tabelama XLII i XLIII. 

a) Linearna korelacija ' 

Primjer: U korelacionoj tabeli XXXVII. str. 98 prika¬ 
zana je raspodjela varijanata za duljinu klasa i za broj klasića 
pšenice. Da ocijenimo jačinu povezanosti istraženih svojstava 
treba grafički prikazati korelaciju između spomenutih svojstava 
pomoću faktične i teoretske regresione linije. Brojčane podatke 
za konstruiranje eksperimentalne regresione linije dobijemo, 
ako za relativno i suponovano svojstvo najprije izračunamo 
prosječne razredne vrijednosti. Na pr. prosječan broj klasića u 


8 Tavčar: Biometrika u poljoprivredi 
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Tab. XLIII. Konstante u jednadžbi za parabolu: 


n 

ka 

a 

P 

V 

■S(n) 

2 (n 2 ) 

3 

2 

10 

12 

3 

6 

14 

4 

4 

5 

5 

1 

10 

30 

5 

14 

7 

6 

1 

15 

55 

6 

112 

28 

21 

3 

21 

91 

7 

84 

12 

8 

1 

28 

140 

8 

168 

15 

9 

1 

36 

204 

9 

924 

55 

30 

3 

45 

285 

10 

528 

22 

11 

1 

55 

385 

11 

858 

26 

12 

1 

66 

506 

12 

4.004 

91 

.39 

3 

78 

650 

13 

2.002 

35 

14 

1 

91 

819 

14 

2.912 

40 

15 

1 

105 

1.015 

15 

12.376 

36 

48 

3 

120 

1.240 

16 

5.712 

51 

17 

1 

136 

1.496 

17 

7.752 

57 

18 

1 

153 

1.785 

18 

31.008 

190 

57 

3 

171; 

2.109 

19 

13.566 

70 

20 

1 

190 

2.470 

20 

17.556 

77 

21 

1 

210 

2.870 

21 

67.298 

253 

66, 

3 

.231 

3.311 

22 

28.336 

92 

23 

1 

253 

3.795 

23 

35.420 

100 

24 

1 

176 

4.324 

24 

131.560 

325 

75 

3 

300 

4.900 

25 

53.820 

117 

26 

1 

325 

5.525 


razredu s duljinom klasa od 6.5 —7 cm izračunamo tako, da 
od 8 varijanata tog razreda potražimo prosječan broj klasića 
na klasu. Dakle: 


Broj 

klasića 

Broj varijanata 
u razredu 

V • p 

11 

1 

11 

12 

3 

36 

13 

2 

26 

14 

2 

28 


00 

II 

a 

101 = 2 Vp 


Prema tome je prosječan broj klasića na klasu u razredu 
sa 6,5—7 cm duljine klasa i s 8 varijanata ovaj: 


S v- p 

2 p 


101 

8 


= 12,625 


Na isti se način izračuna i prosječan broj klasića u svim 
ostalim razredima. Postupak kod izračunavanja prosječne du¬ 
ljine klasa u pojedinim razredima je dakako isti. Na pr. za ra¬ 
zred s 12 klasića na klasu izračunamo prosječnu duljinu klasa 
na ovaj način: 


Duljina klasa 
j V u cm 

Broj varijanata 
u razredu p 

V-p 

/ ». e* 

6 —6,5 

2 

12,50 

6,5-7 

3 

20,25 

7 -7,5 

2 

14,50 

7,5-8 

2 

15,50 

8 -8,5 

1 

8,25 


-£p = 10 

71,00 = 2 Vp 



Prosječna je duljina klasa za razred s 12 klasića na klasu 

ova: 


SVp 71 

2 p “10 


7,1 


Eksperimentalne regresione linije nađemo tako, da na pr, 
za vrijednosti x uzmemo razmake od 1 cm, a onda na ordinate 
nanesemo prosječne razredne vrijednosti za y svojstvo i obrnuto 
I. j. ako na apscisi za y svojstvo nanesemo razmake na pr. od 
1 cm, onda dolaze na ordinate prosječne razredne vrijednosti za 
x svojstvo. (Vidi sliku br. 19). 


Budući da smo za duljinu klasa i broj klasića na klasu pro¬ 
našli linearnu korelaciju, treba za računanje njihovih teoretskih 
regresionih linija upotrijebiti Mikičeve formule i tabele za 
linearnu jednadžbu (vidi str. 112, tab. XLII.). 



* 
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Regresionu liniju za duljinu klasa i z r a č Ut 
namo ovako: 

Duljina klasa uvrštena je u 18 razreda, dakle moramo se 
poslužiti Mikičevom formulom i tabelom za parni 
broj varijanata — r a z r e d a t. j. 

/i- ; . b = i-[2S(xy)-(n+l)2(y)] ' 

... U ovu formulu uvrstimo podatke iz naše tabele na :str ; 112 
i str. 117 za 18 varijanata-razreda: ' '! • 'V 

• b i JL- [2 • (3.493,74) — (18+ 1)- 328,96] = 0,7608 ., .,, 

. • j | *• ; » i ; " .u:. 

a = * jS (yi — 2 (n) • b] = ^ [328,96 — (71 • 0,7608J = 11,048 


.• Iz a i b izračuna se y t. j. ordinate za odgovarajuće x. Kako 
je kod nas regresiona linija pravac, to treba pronaći samo dvije 
točke na ordinati i to na pr. za x = 1 i x = 18. 

Formula za linearnu jednadžbu glasi: y — a 4- bx. 

. Ako je x = 1, onda je ■ ; - 

y = 11,048 + 0,7608 • 1 = 11,8088. 

Za x — 18 je 

y = 11,048 + 0,7608 -18 = 24,7424. 

Iz ovih podataka konstruiramo teoretsku regresionu liniju 
(vidi Si. 19). 


Tab. XL1V. 

Klasovi pšenice 


Razredi za 
duljinu klasa 

Prosječni 
broj klasića 
na klasu 

y 

Broj 

razreda 

x • y 

Teoretski 

y 

6 - 6,5 

11,67 

1 

11,67 

11,8088 ; 

6,5— 7 

12,62 

2 

25,24 

i 

7 — 7,5 

13,25 

3 

39,75 


7,5— 8 

14,17 

4 

56,68 . 


8 — 8,5 

14,88 

5 

74,40 


8,5— 9 

15,51 

6 

93,06 


9 - 9,5 

16,72 

7 

117,04 


9,5 - 10 

17,63 

8 

136,24 


10 -10,5 

18,20 

9 

163,80 


10,5—11 

18,86 

10 

188,60 


11 —11,5 

19,03 

11 

209,33 


11,5— 12 

19,16 

12 

237,12 


12 —12,5 

20,72 

13 

269,36 


12,5-13 

22,07 

14 

308,98 


13 -13,5 

22,30 

15 

334,50 


13,5 — 14 

23,42 

16 

374,72 


14 —14,5 

24,25 

17 

412,25 


14,5 - 15 

24,50 

18 

441,00 

24,7424 


-2y = 328,96 


2J(xy) = 3493,74 



y = a + bx; a = 11,048; b = 0,7608. 
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Na slici br. 19 vidimo teoretsku regresionu liniju i liniju 
dobivenu kod eksperimenata za duljinu klasa. 

Brojčani podaci za regresionu liniju za broj klasića. nalaze 
se u tabeli XLIVa., str. 119. Oni se izračunaju na isti način kao i 
za duljinu klasa samo s tom razlikom, da se ovdje služimo M i- 
kičevom formulom za podatke o neparnom broju razreda 
(str. 112, tab. XLII.), jer je ovdje samo 15 razreda. Za naš pri¬ 
mjer je: 

b = -1 [2 (xy) - k ž 2 y ] = 2^0 C 1378 - 24 ~ 8 X 153 >"1 = °- 5225 

a = — 2 y - k. 2 b =4 X 153,99 - 8 X 0,5225 = 6,086 7- 
n 15 

Po formuli: y = a + bx izračunamo y za x = 1 

y = 6,086 + 0,5225 • 1 = 6,6085. 

Za x — 15 je 

y = 6,086 + 0,5225 • 15 = 13,9235. 

Na apscisi za x = 1 nanesemo ordinatu y = 6,6085, a za 
x — 15 ordinata je y = 13,9235. Ako ove dvije točke na ord'i- 


CaJdavosis 

$ 5 IO 15 IO 25 30 35 W h5 90 59 60 65 X> 75 SO 35 90 99 Uđ 



Tab. XLIV a. 


Razredi za 
broj klasića 
na klasu 

Prosječna 
duljina klasa 

y ' 

Broj 

razreda 

X 

x y 

Teoretski 

y 

11 

6,75 

1 

6,75 

6,6085 

12 

7,10 

2 

14,20 


13 

7,71 

3 

23,13 


14 

8.10 

4 

32,40 


15 

8,79 

5 

43,95 


16 

9,36 

6 

56,16 


17 

9,78 

7 

68,46 


18 

1013 

8 

81,04 


19 

10,69 

9 

96,21 


20 

11,02 

10 

110,20 


21 

11,71 

11 

128,81 


22 

12,22 

12 

146,64 


23 

12,78 

13 

166,14 


24 

13,60 

14 

190,40 


25 

14,25 

*7 15 

213,75 

13,9235 


2 y = 153,99 2 (x y) = 1378,24 


nati spojimo pravcem, dobijemo regresionu liniju (vidi sliku 
broj 19). 

Na si. br. 19 prikazana je teoretska regresiona linija i re- 
gresiona linija konstruirana iz eksperimentalnih data. 

Između caklavosti zrna pšenice i trajanja meljave postoji 
negativna linearna korelacija (Tab. XL. str. 105). Kod grafičkog 
prikaza te korelacije (slika 20) služimo se isto tako Mikičevim 
tabelama i formulama za linearnu jednadžbu. U tom se slučaju re- 
gresione linije protežu u drugom smjeru, nego kod pozitivne 
linearne korelacije. 


SI. 20. 



b) Parabolićna korelacija 

U korelacionoj tabeli XLV., str. 121 prikazana je nelinearna 
korelacija za caklavost zrna i konzistencu tijesta. Sudeći po 
razdiobi varijanata u toj tabeli, korelacija je ovdje približno 
parabolićna. Treba grafički prikazati, do koje se mjere ekspe¬ 
rimentalni podaci podudaraju s pretpostavkom o paraboličnoj 
korelaciji pa prema tomu i o paraboličnoj regresionoj liniji. U 
tu svrhu treba pronaći najprije faktičnu, a onda teoretsku re- 
gresionu liniju. 

Podatke za faktičnu regresionu liniju dobivamo iz pro¬ 
sječnih razrednih vrijednosti za suponovano i za relativno svoj¬ 
stvo. Prosječne razredne vrijednosti izračunamo i ovdje na isti 
način kao i za linearnu korelaciju. 

U tab. XLV, na str. 121 sakupljeni su brojčani podaci za pro¬ 
sječne razredne vrijednosti suponovanog svojstva (x ... konsi- 
stenciju tijesta) i relativnog svojstva (y... caklavost zrna). 

Teoretsku regresionu liniju za suponovano i relativno svoj¬ 
stvo konstruiramo iz ordinata dobivenih po jednadžbi za para¬ 
bolu, koja glasi: 

y = a + bx + cx 2 

Vrijednosti a, b i c računamo po Mikičevim formu¬ 
lama i tabelama (vidi str. 112 i 113) za parabolu i to za 
parni broj podataka, jer smo caklavost i konstituciju uvrstili 
u 20 razreda. 

Za caklavost računamo vrijednosti a, b i c iz podataka u 
tab. XLIII. i XLV. na str. 114 i 121 za 20 razreda i iz tab. XLVI. 
na str. 123. 

Po tome je: ' 

c = [77 X 940,6 — 21 X 11285,4 + 1 X 154.681,6] = 

17,556 

= — 0,56308 

b = l 2 X 11285,4 — 21 X 940,6] — 21 X (— 0,56308) = 

= 13-94362 

a = ^ [940,6 — 210 X 13 94362 — 2870 X ( - 0,56308 =) 

-18 576 


f j 

! 

Prosječna 

caklavost 

u % 

kO 

CM 

0-01 

17-5 

15.833 

40.833 

IO 

TK 

IO 

-rr 

ip 

42-794 

51-5 

49-5 

46-3 

47-368 

52-743 

54-012 

601-85 

CO 

00 

p 

kb 

kO 
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Iz a, b i c izračunamo ordinatne točke y za pojedine vri¬ 
jednosti od x = 1 do x = 20 po formuli za parabolu y = a + 
+ bx + cx 2 . (Vidi kolonu 6 str. 123). 

Ako te točke na ordinatama spojimo, nastaje parabolicna 

regresiona linija (slika 21). 

Caklavost zrna u % 



Usporedimo li ovu teoretsku liniju s regresionom linijom 
dobivenom iz eksperimentalnih podataka (ordinate za prosječne 
razredne vrijednosti), naći ćemo, da se ove donekle po u ara J u - 
Na isti način izračunat ćemo i ordinate za regresionu liniju 
za konzistencu tijesta. U tu svrhu nam služi tabe a na 


Tab. XLVI. 


1 

Caklavost 

u °/o 

2 

Prosječna 
konzistenca 
tijesta 
u minutama 

3 

Broj 

razreda 

X 

4 

x-y 

■ 

5 

x 2 y 

6 

Teoretski : 

(y) 

■ minute 

0-5 

2,5 

1 

2,5 

2.5 

-5-2 

5 — 10 

7,5 

2 

15,0 

30,0 

6-1 

10 — 15 

12,5 

3 

37,5 

112,5 

18-2 

15-20 

24,4 

4 

97,6 

390.4 

28-2 

20 — 25 

36,6 

5 

183,0 

915,0 

371 

25-30 

49,1 

6 

294,6 

1767,6 

44-8 

30 — 35 

53,4 

7 

373,8 

2616,6 

51-4 

35—40 

54,5 

8 

436,0 

3488,0 

56-9 

40 — 45 

55,5 

9 

499,5 

4495,5 

61-3 

45 — 50 

59,3 

10 

593.0 

5930,0 

64-6 

50-55 

63,6 

11 

699,6 

7695,6 

667 

55-60 

66,4 

12 

796,8 

9561,6 

67-7 

60-65 

68,6 

13 

891,8 

115934,5 

67-5 

65 — 70 

70,5 

14 

987,0 

13818,0 

66-3 

70-75 

71,1 

15 

1066,5 

15997,5 

639 

75 — 80 

65,5 

16 

1048,0 

16768,0 

60-4 

80 — 85 

64,6 

17 

1098,2 

18669,4 

55-7 

85-90 

47,5 

18 

855,0 

15390,0 

499 

90-95 

40,0 

19 

760,0 

14440,0 

431 

95 - ICO 

27,5 

20 

550,0 

11000,0 

351 


940,6 = 

210 = 

11285,4 = 

154681,6 = 



= ^(y) 

= JS(x) 

= -S(xy) 

= £(x’y) 



str. 124. Pokusom dobivene ordinate sadržava kolona 2, a teoret¬ 
ske pak kolona 6. 

Spojimo li (slika 21) ordinatne točke iz kolone 2, nastaje 
eksperimentalna regresiona linija. Do koje se mjere ona podu¬ 
dara s teoretskom, vidimo, ako točke na ordinatama iz kolone 
6 spojimo linijom. 

U našem slučaju je bila opravdana pretpostavka, da su kon- 
zistenca tijesta i caklavost zrna u paraboličnoj korelaciji. 
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1 1 

2 

Konzistenca 
tijesta u 
minutama 

0 

o 

© 

cd ^ 

KJ -*-< 

0) (0 

V g * 

£ 5 

0- as 

cO 

O 

0'— 5’ 

2'500 

5’ — 10’ 

lO’OOO 

10'- 15' 

17-500 

15' - 20' 

15-833 

20' - 25' 

40-833 

25' — 30' 

41-500 

30' — 1 35'. 

41'545 

35' - 40' 

42-794 

40’ — 45' 

51-500 

45' — 50' 

49-500 

50' — 55' 

46-300 

55' - 60' 

47-368 

60' — 65' 

52-743 

65' — 70' 

54-012 

70 — 75' 

58-109 

75' - 80' 

55"086 

80' — 85' 

59500 

85' — 90' 

56-730 

90' — 95' 

58-928 

95' — 100' 

65-500 



2’500 

20-000 

52’SOO 

63‘332 

244-165 

249-000 

290-815 

342-352 

463-500 

495-000 

509-300 

568-416 

685-659 

756-168 

871-635 

881-376 

IOH’500 

1021-140 

1119-632 

1310-000 


2T(xy) = 

= 10917-99 



y = a-(- bx + Cx2 

C = -i. [«(Xy) - ps (xy) + 7-S (x 2 y)] = 

_1_ [77 + 867-781 — 21 X 10917-99 + 1 X 159.012,954] 

17.556 1 

== — 0*19626 


f 
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b = ^[22(xy)-(n + l)2y]-(n + l)c = 

== im [2X 10-917-99-(2° +1) X 867781]- (20 + 1) X 
X -0-19626 = 6*83766 
a = -J- [2 (y) - 2 (n) • b — 2 (n 2 ) c] = 

~ [867-781 —210 X 6-83766 — 2870X-019626] = — 0-24306 




C. BIOMETRIĆKA ISTRAŽIVANJA 
KVALITATIVNIH IU 
ALTERNATIVNIH SVOJSTAVA 


0—1. BIOMETR1KA VARIJACIJE 
Općenito 

Do sada smo proučavali varijaciju svojstava cijelih (di¬ 
skretnih) i razrednih kvantitativnih varijanata. 
U varijacionom redu takvih svojstava razlike su izražene kao 
kvantitete. 

U eksperimentalnoj poljoprivredi imamo vrlo često posla 
i s kvalitativnim svojstvima. Stoga je potrebno, da 
i ova svojstva varijaciono-statistički obradimo, kako bismo iz 
izračunatih vrijednosti stvorili pravilne zaključke. 

Kvalitativna svojstva imaju alternativnu 
varijaciju. Na pr. jednospolni su organizmi ili ženski ili 
muški. 

Naročito mnogo kvalitativnih varijacija pojavljuje se u po¬ 
tomstvu križanaca. Na pr. u drugoj filijalnoj generaciji kod 
križanja homocigotnog (samooplodnja) graška žutog sjemena 
s homocigotnim graškom zelenog sjemena, razviju se u F 2 ge¬ 
neraciji biljke sa žutim i s zelenim sjemenom. Dakle imademo 
dvije alternative. 

Križanjem pšenice-golice žutih pljevica s pšenicom brku- 
Ijom smeđih pljevica izrasle su u F 2 generaciji biljke u četiri 
tipa i to: a) smeđa golica, b) žuta golica, c) smeđa osata i d) 
žuta osata pšenica — dakle četiri alternative. 

Isto je tako i kod životinja: 

Na pr. u F 2 generaciji križanjem homocigotnih bijelih 
Leghorn kokoši golih nogu s drugima crne boje i obraslih nogu 
razviju se životinje: 
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a) bijele boje, obraslih nogu, b) bijele boje, golih nogu, 
c) crne boje obraslih nogu i d) crne boje, golih nogu. I ovdje 
su četiri alternative. 

Nadalje nalazimo alternativna svojstva kod di- i polimor- 
fizma, kao i u svim onim slučajevima, gdje treba u populaciji 
ustanoviti omjer pojedinih kvalitativnih svojstava, koja sači¬ 
njavaju smjesu. ...... .v 

Prema tome imamo kod alternativnih svojstava dvije ili vise 

alternativa. .... 

I alternativna svojstva dadu se i treba da budu varijaciono- 

statistički obrađena. 

Primjenu varijacione statistike na kvalitativna svojstva pn- 
kazat ću na nekoliko primjera iz svojih genetskih istraživanja. 

Iz genetike (nauke o nasljeđivanju) poznato je, da su geni 
t j. osnove za razvitak izvjesnog svojstva, lokalizirani na hro- 
mosomima, te da se sa spolnim stanicama prenose od roditelja 
na potomstvo. 

Koliko ima izvjesna odlika (varijetet) organizma različitih 
hromosoma u spolnim stanicama (haploidni broj hromosoma), 
toliko će biti i skupina gena, koji se križanjem organizma ne¬ 
ovisno kombiniraju s genima drugog organizma iste odlike 

Na pr. dominantan se gen Ž za žutu boju sjemena graška 
nalazi na jednom kromosomu, a dominantan gen O za okrug i 
oblik sjemena na nekom drugom hromosomu u spolnoj stanici 
roditelja žutog i okruglog sjemena. Genetska je konstitucija 
takvog homocigotnog roditelja ŽŽOO. Ako ga križamo s ro¬ 
diteljem genetske konstitucije žžoo t. j. s recesivmm svojstvima 
sjemena zelene boje i naboranog oblika, bit će F x gen. genetske 
konstitucije ŽžOo. U ovoj se generaciji stvaraju 4 različite ga- 
mete i to ŽO, Žo, ŽO i žo u omjeru 1 :1 : 1 : 1. 

U Fo generaciji izrast će ćetiri fenotipa i to konstitucije 
ŽO, Žo, žO i žo u omjeru 9:3 : 3 :1. 

Jednaki se fenotipovi i u istom omjeru kao gore stvaraju 
u F 2 generaciji, kad križamo homocigotne roditelje žutog i na¬ 
boranog sjemena, dakle genetske konstitucije ŽŽoo, s ro jte 
ljima zelenog i okruglog sjemena t. j. ŽŽOO. U tom je s ucaju 
u svakog roditelja po jedno dominantno i jedno lecesivno svoj 
stvo. Dakle je svejedno, da li povratno križamo F^ generaciju, 
nastalu križanjem roditelja genetske konstitucije Z ^ .zzoo 
ili križanjem roditelja genetske konstitucije ZZoo. zzUV, Fi 
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generacija konstitucije ŽžOo stvara gamete ŽO, Žo, ŽO i žo u 
omjeru 1 :1 :1 : 1. 

Za svojstva, kojih se geni nalaze u istom hromosomu, gore 
spomenuti omjeri ne vrijede. 

Ako križamo homocigotnog roditelja, kod kojega su oba 
dominantna gena lokalizirana u istom hromosomu (AB/AB), s 
drugim roditeljem, u kojega se oba recesivna gena nalaze za- 
jedno-u homolognom hromosomu (ab/ab), onda je Fi genera¬ 
cija konstitucije AB/ab. 

Ona stvara gamete AB, Ab, aB i ab u omjeru n : 1 :1 : n i 
to tako, da je n > 1. Ovakvo križanje spada u seriju spa¬ 
janja (coupling serie). 

Fi generacija Ab/aB, koja se razvila od križanja roditelja, 
U kojima je po jedno dominantno i po jedno recesivno svoj¬ 
stvo, dakle Ab . Ab X aB . aB, stvara gamete AB, Ab, aB i ab 
u omjeru 1 : n : n : 1. 

Takva križanja zovemo razdvajanjem, a skupinu, u koju 
spada ovo križanje, zovemo serijom razdvajanja (re- 
pulsion serie). 

Poradi razlike u omjeru različitih gameta dakako da su i 
omjeri fenotipova u F 2 generaciji u seriji spajanja dru¬ 
gačiji od omjera u seriji razdvajanja. 

Ima više metoda za varijaciono-statističko proučavanje po¬ 
dataka o svojstvima organizma u različitim generacijama nakon 
križanja i to: 

A) uspoređivanje eksperimentalnih i teoretskih brojaka i 
mjerenje razlike pomoću srednjih odnosno vjerojatnih pogre¬ 
šaka i to 1. za relativne i 2. za apsolutne brojke, 

B) izračunavanje, da li su razlike između eksperimentalno 
dobivenih i teoretskih brojaka u granicama vjerojatnosti P 
(metoda x 2 ). 

ad A) Primjenu metode A—1 i A—2 prikazat ćemo na sli¬ 
jedećem primjeru: 

Jari rijetko-klasi šesteroredni ječam golog zrna križali smo 
s ozimim dvorednim ječmom obuvenog zrna da u F 2 i daljnim 
generacijama dobijemo po koju biljku ozimog dvorednog ječma 
golog zrna. 


1 


9 TavEar: Biometrika u poljoprivredi 
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U F 2 generaciji izraslo je u konkretnom slučaju 1020 bi¬ 


ljaka ovih svojstava: 

a) dvoredac obuvenog zrna.583 

b) dvoredac golog zrna .197 

c) rijetkoklasi šesteroredac obuvenog zrna . 184 

d) rijetkoklasi šesteroredac golog zrna ... 56 


1020 


I. DVIJE ALTERNATIVE 


Treba da istražimo, kako se nasljeđuje broj redova 
na klasu i u kojem su odnosu obuvena i gola 
zrna, i kako ta svojstva variraju. 

Najprije nam valja ustanoviti broj biljaka s jednim svoj¬ 
stvom prema broju biljaka s drugim svojstvima. 

Kod alternativnih svojstava izaberemo jednu alternativu 
— svejedno koju — kao mjernu jedinicu-. Omjer pojedinih alter¬ 
nativa izražavamo ili u postocima ili u dijelovima jedinice. 

Treba izračunati najprije omjer pojedinih alternativa, da¬ 
kle srednje vrjednosti M, a onda varijabilitet pojedinih svoj¬ 
stava t. j. standardnu devijaciju o kao temelj srednje pogreške. 

Zanima nas najprije, kako se nasljeđuje broj redova na 
klasu ječma. Za ova svojstva postoje dvije alternative i to: 
a) biljke sa dvorednim i b) biljke sa rijetkokla- 
sim šesterorednim klasovima. Jednu alternativu uz¬ 
memo kao mjernu jedinicu. Prema tome ima druga alternativa 
u donjem varijacionom redu vrijednost ništice. Onda je vari- 
jacioni red ovakav: 



1. Srednja vrijednost 


Iz toga izračunamo najprije srednju vrijednost 
(M), a onda pomoću standardne devijacije (u) još i srednju po¬ 
grešku (m) srednje vrijednosti. 
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Srednju vrijednost kvantitativnih svojstava izračunamo po 

21 D cl 

spomenutoj već formuli M = A + B; b =—-—Iz toga pro¬ 
izlazi, da je M = — formula za srednju vrijednost kvalitativnih 
n 

svojstava. 

240 

U našem je slučaju: M = jq 2 q — 0,2353, a to znaci, da ima 

0,2353 ili 23,53% biljaka rijetkoklasnog šesteroreca. 

Uzmemo li za mjernu jedinicu razred sa biljkama dvored- 
nih klasova, onda je vrijednost druge alternative ravna ništici 

JOf) 

i M = = 0.7647 ili 76.47% biljaka s dvorednim klasovima. 

Budući da je kod križanja tip biljaka dvorednog klasa 
dominantan nad tipom biljaka rijetkoklasnog šestero¬ 
reca, razvijaju se u F 2 generaciji dvored i rijetkoklasni še- 
steroreci u teoretskom omjeru 3 : 1. 

Kombinacioni je broj: K = 3 + 1 = 4. 

Zanima nas, da li je u našem pokusu dobiveni omjer bio- 
metrički opravdan ili nije. 

Ima dva načina da to utvrdimo, i to tako da: 
a) kompariramo teoretske relativne "brojke s relativnim 
brojkama, koje smo izračunali iz podataka dobivenih u pokusu, 
ili 

ft) kompariramo teoretske apsolutne brojke s brojkama 
iz pokusa. ■, 

ad a) Relativne brojke eksperimenta izračunavamo tako, da 

K 

apsolutne brojke, koje smo dobili pokusom, pomnožimo sa —, 

gdje je n apsolutni broj svih varijanata u pokusu, a K kombi¬ 
nacioni broj izveden iz formule (1 + 2 + 1)", odnosno (3 + l) n 
za omjer fenotipova u F^ generaciji, ako je n broj hetero- 
cigotnih genetskih parova t. j. 

K = 4 n ; dakle 3:1 = (3+1) 1 = 4; 9:3:3:1 = (3 + 1) 2 = 16; 
’ 27 : 9 : 9 : 9 : 3 : 3 : 3 : 1 = (3 + l) 3 = 64 i t. d. 
ad ft) Teoretske apsolutne brojke izračunat ćemo tako, da 

v. n 

teoretske relativne brojke na pr. 3 : 1 pomnožimo sa 

ad a) Apsolutne brojke fenotipova naše F 2 generacije naj¬ 
bolje odgovaraju omjeru 3:1 t. j. kombinacionom broju K—4. 





Prema tome su relativne brojke fenotipova našeg eksperi¬ 
menta ove: 

780 x mo = =3,059 biljaka sa dvorednim kIasom > 

240 x t4o = m 6 =0,941 sa šesteroredn5m klasom 

3,059 odgovara teoretskoj relaciji 3, a 

0,941 teoretskoj relaciji 1. 

ad jS) Teoretske apsolutne brojke su slijedeće: 

3 1020 ^ 3060 = 765 biljaka dvorednog klasa, 

4 4 

1 X 1^2 = 255 biljaka šesterorednog klasa. 

780 odgovara teoretskom broju 765 
240 odgovara teoretskom broju 255. 


2. Standardna devijacija 


Želimo li izračunati srednju pogrešku srednjih vrijednosti 
za relativne odnosno apsolutne brojke, valja nam najprije izra¬ 
čunati standardnu devijaciju o. 

Standardnu devijaciju za kvantitativna svojstva naći ćemo 
po formuli: 




Kod kvalitativnih svojstava je razredni areal: a — 1. Već 
smo prije spomenuli, da je kod kvalitativnih svojstava b = 0l - 

Iz toga slijedi, da je b 2 = tako da je onda: 


~/Pi_Pi 2 V'PoPi 
1 n n 2 n ’ 


jer je po + Pi = n. 

„ . 1 / 780 X 240 nA9A17 

U našem je slučaju o-1020~ 
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Dakle: 

Mi = 3,059 .... dvoredac 

M 2 = 0,941 .... šesteroredac 

o = 0,4241 .... dvoreca ili šesteroreca. 

Standardnom devijacijom mjerimo varijabilitet svojstva i 
kod kvalitativnih svojstava. Najveći je varijabilitet onda, kad 
je kod dviju alternativa a = ± 0,5. Dakle kreće se od 0 — 0,5. 

Srednje vrijednosti izražavamo u postocima, a isto tako i 
standardnu devijaciju. Formule za izračunavanje standardne 
devijacije u postocima su ove: 


_ IPOj/poP, 


, odnosno a = 


100 po 100 p, 

n X n 


ili a = |/ % P 0 X %P| 


U našem slučaju: 


a = 1/76,47 X 25,53 = 44,183% 
Prema tome: 


Mi — 76,47% dvoreca 
M 2 = 25,53% šesteroreca 
a = 44,18% dvoreca ili šesteroreca. 


Kod dviju je alternativa standardna devijacija u granicama 
od 0 — 50%. Ovdje se varijacioni koeficijent ne upotrebljava 
kao relativno mjerilo, jer je već standardna devijacija relativna 
vrijednost. 


3. Srednja pogreška srednje vrijednosti 

Na temelju srednje pogreške srednje vrijednosti kvantita¬ 
tivnih varijanata izvedene su formule za srednje pogreške kva¬ 
litativnih ili alternativnih svojstava. 


m — p ,— je općenita formula za srednju pogrešku srednje 


vrijednosti. Srednja pogreška srednjih vrijednosti M = 0.7647 
dvoreca odnosno M = 0,2553 šesteroreca izračunamo po 


formuli: 



1 1 j 780 X 240 

1020 I 1020 


= 0,0132 


za dvorece odnosno šesterorece. 
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Kako je trostruka srednja pogreška manja od naših srednjih 
vrijednosti Mi = 76,47, odnosno M 2 = 25,53, te su vrijednosti 

varijaciono-statistički opravdane. 

B) Ako srednje vrijednosti i standardnu devijaciju izrazimo 
u postocima, onda treba izraziti u postocima i srednju pogrešku. 
Formula tada glasi: 

m „y %Pox«£pr 

U našem primjeru: 

76,47 X 25^53 = 13834% 

1020 

Mi = 76.47°/o dvoreca 

M 2 = 25.53% šesteroreca 

m = ± 1.3834% dvoreca ili šesteroreca. 

I ovom je metodom potvrđeno, da su u našem primjeiu 
srednje vrijednosti M t i M 2 varijaciono-statistički opravdane, 
jer je trostruka srednja pogreška manja od svake ove vrijed¬ 
nosti. . 

«) Kad proučavamo omjer različitih fenotipova u h 2 gene¬ 
raciji, najbolje je, ako nam za temelj služi kombinaciom broj K. 

U našem primjeru ima 780 biljaka s dvorednim klasom i 
240 biljaka sa šesterorednim klasovima. 

Omjer 780 : 240 se najbolje slaže s teoretskim omjerom 3.1 
kombinacionog broja K = 3 + 1 = 4. 

Valja nam istražiti, da li pokusom dobiveni omjer varija¬ 
ciono-statistički odgovara teoretskom omjeru 3 : 1. 

Moramo izračunati srednju pogrešku spomenutog omjera. 
Srednju pogrešku izračunamo pomoću standardne devijacije a. 

U ovom slučaju upotrebit ćemo formulu: o = 100 
i izračunati a obzirom na kombinacioni broj K, -dakle. 

„ _ K VIv pT 
n 

Ako uzmemo K za bazu, onda treba kod računanja uvijek 
upotrebiti samo dvije alternative. Na pr. kod kombinacionog 
broja K = 4 radimo s omjerom 3:1 ili 2 : 2; kod K 16, ako 
je cjepanje 9 :3 :3 : 1 s omjerom 9:7, 13:3 ili 15 : 1, ako je 
cjepanje 9:3:4 operiramo s omjerom 9 . 7 i 12 i . • , 




Jasno je, da je za pojedine alternativne relacije i standardna 
devijacija različita. 

Ako je K kombinacioni broj, a jedna alternativa označena 
s N, onda je druga (K — N). 


Prema tome je: 

K ]/Po Pi _ K ]/N (K — N) __ j/n _N) 

n K 

Kako je općenita formula za srednju pogrešku srednje vri¬ 
jednosti: m = ~=, to je u našem slučaju m = f N ( K ~ 

Treba sada istražiti, u koliko se omjer biljaka 780 dvore- 
daca : 240 šesteroredaca podudara s idealnim omjerom 3:1. 


Za 3 je 0,00294 

Za t je m — j/—J (|2 jh " °' 00294 

I ovdje vrijedi pravilo, da je srednja vrijednost statistički 
opravdana, ako je veća od trostruke srednje pogreške. U našem 
je slučaju trostruka srednja pogreška manja od srednjih vrijed¬ 
nosti 3 odnosno 1. Dakle srednje su vrijednosti varijaciono- 
statistički opravdane. 

ć) U gornjem smo primjeru komparirali teoretske rela¬ 
tivne brojke s eksperimentalno pronađenim relativnim brojka¬ 
ma. Ima i drugi način, a taj je, da apsolutne brojke dobivene kod 
pokusa uspoređujemo s teoretski očekivanim apsolutnim broj¬ 
kama. 


Teoretske apsolutne brojke izračunali smo već prije po 

formuli -H-. Srednju pogrešku ovih brojeva naći ćemo tako, da 

postavimo dvije alternative: jedna je apsolutna brojka (q) vari¬ 
janata za jedno svojstvo, a druga je zbroj svih ostalih varija¬ 
nata za drugo svojstvo (n—q), ako je n zbroj svih varijanata. 


m = 
brojaka. 


2-lH-Hl je formula za srednju pogrešku apsolutnih 
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Za teoretsku apsolutnu brojku M = 780 je m = 

_ 1 / r 780< 1020 - 780) = 1s ^ 

1020 

I ovdje je trostruka srednja pogreška manja od srednjih 
vrijednosti M = 780 odnosno 240. 

II. VIŠE ALTERNATIVA 

U gornjem smo primjeru varijaciono-statistički istražili 
dvije alternative. Valja nam proučiti i one slučajeve, gdje ima 
više alternativa, a to je vrlo često baš u poljoprivrednoj genetici. 

Kako je već spomenuto, križao sam jari rijetkoklasi še- 
steroredni ječam golog zrna s dvorednim ozimim ječmom obu* 
venog zrna. 

Dosada smo istražili, kako se nasljeđuje broj redova na 
klasu. Proučit ćemo, kako se nasljeđuje i golo, odnosno obu¬ 
veno zrno. 

Za broj redova u našem su primjeru 2 alternative, a za svoj¬ 
stva zrna (golo odnosno obuveno) opet 2 alternative. 

U takvom slučaju treba svaku pojedinu alternativu prouča¬ 
vati u odnosu prema svim ostalim alternativama. 

Nas interesira, u kojem se omjeru pojavljuju četiri spome¬ 
nuta fenotipa u F 2 generaciji. 

a) Najprije želimo usporediti teoretske relativne brojke s 
relativnim brojkama dobivenim kod pokusa. Izračunat ćemo 
to tako, da apsolutne brojke dobivene eksperimentalnim putem 
K 

pomnožimo sa —. 

U našem se primjeru (str. 130) pokusom dobivena relacija 
fenotipopa najviše približava teoretskoj relaciji 9 : 3 : 3 : 1 s 
kombinacionim brojem K = 16. 

Teoretskim relativnim brojevima 9 : 3 : 3 : 1 odgovaraju 
pokusom dobivene relativne brojke i to: 


Za 9. «XI« , 9,145 dvoredaca obuvenog zrna 

Za 3...—' [ = 3,0901 dvoredaca golog zrna 




I 


Za 3 
Za 1 


184 X 16 _ 2 8g6 gesteroredaca obuvenog zrna 
1020 

56 X 16_ — 0.8784 šesteroredaca golog zrna. 
1020 


Srednje su pogreške za teoretsku relaciju po formuli 


m = 


m = 


m = 


m = 


/ N (K 

N) 

i " 


/ 9X7 

= 0 , 

jf 1020 

1/3 X 13 


1/ 1020 


1/1 X 15 


|f 1020 



ove: 


Budući da su srednje vrijednosti veće od trostruke srednje 
pogreške, to je naša teoretska pretpostavka, da se u F 2 gene¬ 
raciji stvaraju pojedini fenotipovi u omjeru 9 .3.3 . 1 varija- 

_D_ 

ciono-statistički opravdana. To proizlazi i iz relacije mD > 

(u tab. XLVIII.) koje su manje od 3, a to znači, da razlika iz¬ 
među teoretskog omjera 9 : 3 : 3 : 1 i pokusom dobivenog omjera 
nije statistički opravdana. 


Tab. XLVIII. 


i Svojstva 

Broj bi¬ 
ljaka u 
pokusu 

A 

omjer za 
K = 16 

B 

očekivano 
za K = 16 
M ± m 

(A—B) = D 

_D_ 

m 0 

dvoredac 
obuveno zrno 

583 

9,145 

9 ± 0,2487 

0,145 

0,583 

dvoredac 
golo zrno 

197 

3,090 

3 ± 0,1955 

0,090 

0,460 

šesteroredac 
obuveno zrno 

184 

2,886 

• 

3 ± 0,1955 

— 0,113 

- 0,578 

šesteroredac 
golo zrno 

56 

0,878 

1 ± 0,1213 

— 0,122 

-1,005 
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b) Za uspoređenje apsolutnih teoretskih i pokusom dobi¬ 
venih brojaka služimo se slijedećom metodom: 

Pokusom dobivene apsolutne brojke 583 : 197 : 184 : 56 
(str. 130) odgovaraju ovim teoretskim apsolutnim brojkama 

izračunatim na temelju 

K 


Za 

583 .. . 

■ ■ »x>7- 

573,75 

Za 

197 .. . 

. . . 3XT = 

191,25 

Za 

184 .. . 

o V 1020 _ 

. . . 3X ig 

191,25 

Za 

56 . . . 

1 x 1020 - 
• • • X 16 

63,75 


Srednje su pogreške za teoretske apsolutne brojke po for¬ 
muli m = -9) ove: 

i n 


. 1/~573,75 X 446,25 1CC 

za q = 573,75 je m = / - 1 q - ^q -—— = 15,8 


za q = 191,25 je m 


za q = 63,75 je m 


= 191,25 ) 

. f 63,75 X 
‘ 102 


X 828,75 


= 12,465 


75 X 956,25 


= 7,7304 


Apsolutne srednje vrijednosti ovdje su veće od trostruke 
srednje pogreške, pa je naša teoretska pretpostavka, da se u F 2 
generaciji stvaraju pojedini fenotipovi u omjeru 573,75 : 191,25 : 
: 191,25 :63,75 varijaciono-statistički opravdana. To pro¬ 
izlazi i iz relacije (u tab. XLIX.). Ove su relacije manje od 3, 

a to znači, da razlike između teoretskih apsolutnih brojaka 
573,75:191,25:191,25:63,75 i pokusom dobivenih brojaka 
583 : 197 : 184 :56 nisu statistički opravdane. 
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Tab. XLIX. _ ' 

_D_ 

m D 

0,583 
0,461 
-0,581 
-1,005 


III. VJEROJATNOST — PROBAB1LITET — P. 

X 2 _ METODA 

Do koje su mjere razlike između eksperimentalno dobive¬ 
nih i teoretskih brojaka u granicama vjerojatnosti P, lzraču- 
namo pomoću x 2 metode. 

x 2 je suma (2) kvocijenata kvadrata razlika eksperimen¬ 
talnih (e) i teoretskih (t) brojaka, podijeljeno s teoretskim 
brojkama (t). 

( č ti 2 

Dakle: x 2 = 2 -—^— 

Što se eksperimentalno dobiveni brojevi bolje slažu s teo¬ 
retskim, to je i x 2 manji. Izvjestan rezultat dobiven u pokusu 
ili opažanjem, može se više ili manje podudarati s teoretskom 
pretpostavkom. Uzroci velikih razlika medu spomenutim bi oj 
kama mogu biti različiti. 

Kod monohibridnog cijepanja očekujemo u F 2 generaciji 
fenotipove u teoretskom omjeru 3:1. Ako je na pr. kod pokusa 
ustanovljen apsolutni omjer fenotipova 84 : 30 t. j. 2,8 : 1', onda 
smatramo, da odgovara našoj pretpostavci t. j. teoretskom 
omjeru 3:1. Da smo kod pokusa pronašli omjer 3900 : 1000 
t. j. 3,9 : 1, onda taj omjer ne odgovara cijepanju 3:1, a na¬ 
stao je zbog različitih poremećenja. Jedno od takvih je u tome, 
što je životna sposobnost individua recesivnog razreda u nekih 
organizama slabija od individua dominantnog razreda. Pa i 
onda, kad ne bi bilo poremećenja, možemo teoretski očekivati 
u više bilijuna cijepanja omjer 3,9 : 1 samo jedamput. 


A 

B 


Broj biljaka 

očekivano 

A — B = 

kod pokusa 

M ± m 


583 

573,75 ± 15,843 

9,25 

197 

191,25 ± 12,465 

5,75 

184 

191,25 ± 12,465 

-7,25 

56 

63,75 ± 7,7304 

-7,75 






140 

Da prosudimo, slaže li se naša pretpostavka (teoretske 
brojke) s eksperimentalnim podacima ili opažanjima, i do koje 
se mjere s njima podudara, treba utvrditi granicu, kod koje 
razlika između jednih i drugih brojaka nije više slučajno 
»moguća«, nego je posljedica nekih poremećenja. 

Spomenuli smo, da pomoću x 2 određujemo, do koje se 
mjere podudaraju teoretske i empirijske brojke. Može se izra¬ 
čunati vjerojatnost — probabilitet P, kada će x- 
prekoračiti izvjesnu teoretsku vrijednost — pretpostavku. 

Svakom x 2 odgovara izvjesna vrijednost za P. Ako se x 2 
diže od 0 do neizmjerno, P se smanjuje od 1 do 0. Obično sma¬ 
tramo da je razlika između teoretske i empirijske srednje vri¬ 
jednosti onda signifikantna, ako je veća od trostruke srednje 
pogreške. 

U biometrici najčešće uzimamo za graničnu vrijednost 
p = o.Ol (= 1%). To znači, ako je razlika između teoretskih i 
empirijskih brojaka takva, đa je u 100 slučajeva možemo oče¬ 
kivati jednom, onda kažemo, da se eksperimentalni podaci po- 
duraju s teoretskim. 

Kod genetskih istraživanja služimo se graničnom vjero¬ 
jatnošću P = 0,01 obično onda, kad želimo odrediti, da li je 
pokusni materijal heterogen ili nije. • 

Kad uspoređujemo empirijske brojke s teoretskim u poje¬ 
dinim generacijama nakon križanja, onda se služimo granič¬ 
nom vjerojatnošću P = 0,05. Dakle, ako je razlika između teo¬ 
retskih i empirijskih brojaka takva, da se u 100 slučajeva po¬ 
javljuje barem 5 puta, onda smatramo, da su eksperimentalni 
podaci u skladu s teoretskim. 

Izračunavanje x 2 i vjerojatnosti — P. 

Računski postupak za x 2 prikazat ćemo na slijedećem pri¬ 
mjeru: 

U F 2 gen. od križanja obuvenog ječma crnih pljevica s go¬ 
lim ječmom žutih pljevica izrasli su fenotipovi u- ovom omjeru. 

298 biljaka obuvenog i crnog zrna, 

85 biljaka obuvenog i žutog zrna, 

92 biljke golog i crnog zrna, 

27 biljaka golog i žutog zrna. 


Treba istražiti, da li ovaj omjer odgovara cijepanju 9:3:3:1. 
Da izračunamo x 2 valja nam ustanoviti razliku (e-t) između 
eksperimentalno (e) dobivenih i teoretskih (t) brojaka, i to za 
svaki fenotip posebno. Najprije ćemo svaku od ovih razlika 


Tab. L. 


Svojstva zrna 

eks. - e 

teoretski -t 

e-t 

(e-t) 2 

(e-t)*: t 

obuveni, crni 

298 

282,375(9) 

+15,625 

244,1406 

0,8645 

obuveni, žuti 

85 

94,125(3) 

- 9,125 

83,2656 

0,8846 

goli, crni 

92 

94,125(3) 

— 2,125 

4,5156 

0,0470 ‘ 

goli, žuti 

27 

31,375(1) 

- 4,375 

19,1406 

0,6100 

Ukupno 

502 

502 



2,406 = 


kvadrirati, a onda dividirati s odgovarajućom teoretskom broj¬ 
kom (t). Suma svih ovih koeficijenata je tražena vrijednost x 2 . 
U našem je slučaju x 2 = 2,406. 

Nadalje treba izračunati, koliko je u našem slučaju ne¬ 
ovisnih brojčanih skupina (degrees of freedom). 

Uspoređivali smo četiri empirijske skupine brojaka s četiri 
teoretske. Od četiri empirijske brojčane skupine tri su neovisne, 
a četvrta rezultira iz ostalih-triju. Dakle 4—1=3 ili općenito: 
broj neovisnih skupina (degrees of freedom) = n — 1. U Fi- 
sherovoj tabeli LI. str. 142 naći ćemo u rubrici za 3 neovisne 
skupine, da je od x 2 =3,665 vjerojatnost P = 0,3, a za x 2 =2,366 
vjerojatnost P = 0,5. Prema tome je za naš x 2 =2,406 vjerojat¬ 
nost P=0,49. To znači, da možemo u 100 naših F 2 generacija 
očekivati barem 49 s isto takvom razlikom, koju smo dobili u 
našem primjeru između empirijskih i teoretskih podataka. Pre¬ 
ma tome su naši empirijski podaci u skladu s teoretskim. 

Da je P bio 0,01, onda bismo izvjesnu razliku u 100 F 2 gen. 
morali očekivati najviše jednom. Budući da se u tom slučaju 
empirijski podaci ne bi slagali s teoretskim, značilo bi, da je 
neslaganje uvjetovano različitim smetnjama. 

Da je vjerojatnost P = 0,999, trebalo bi takav rezultat pri¬ 
miti s velikom rezervom, on bi značio, da će se u hiljadu slu¬ 
čajeva izvjesna razlika pojaviti 999 puta. U takvom slučaju treba 


Ukupno: 502 



Vjerojatnost 


142 
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istražiti, nema li izvjesnih računskih ili metodskih pogrešaka. 
Ako se konačno uvjerimo, da ih nema, onda je takvo poduda- 
ranje između empirijskih i teoretskih podataka moguće samo 
onda, ako pokus obuhvaća vrlo velik broj ■ individua. 

Kad su kod dihibrida geni za dva para svojstava neovisni 
t. j. lokalizirani svaki u drugom nehomolognom hromosomu, 
stvaraju se u F 2 generaciji četiri različna fenotipa u omjeru 
9 :3 :3 : 1. Uzroci, što se katkada omjer fenotipova u F 2 gene¬ 
raciji mnogo razlikuje od teoretskog omjera 9 :3 :3 : 1 i što je 
vjerojatnost P < 0,01 mogu biti slijedeći: 

1. pokusni su podaci heterogeni, 

2. poremećenja: a) u cijepanju jednog para svojstava uvje¬ 
tovanih genima Ab, ili 

b) u cijepanju drugog para svojstava, koje uvjetuju geni 

Bb, 

3. geni za dva različita para svojstava lokalizirani su u 
istom hromosomu. 

ad 1. Istraživanje pokusnih podataka na 
heterogenost. 

Rezultate pokusa, koje želimo varijaciono-statistički obra¬ 
diti, sačinjavaju obično podaci od više pojedinačnih podjednakih 
pokusa na pr. od više F x gen, F 2 gen, od povratnih križanja 
(Aa Bb X aa bb) i t. d. Iz podataka dobivenih pokusom možemo 
stvoriti pravilne zaključke samo onda, ako varijaciono-stati¬ 
stički obradimo homogeni materijal. 

Da li izvjesni podaci dolaze od homogenog pokusnog ma¬ 
terijala ili ne, možemo najbolje istražiti pomoću x 2 metode. 

Ako je pokusni materijal homogen, podaci pojedinačnih 
pokusa ne smiju se razlikovati od teoretskih podataka više nego 
što dopušta granična vjerojatnost t. j. P = 0,01. 

Kod homogenog materijala možemo podatke svih pojedi¬ 
nih pokusa proučiti zajedno. 

Ako je pokusni materijal posve heterocigotan, ne možemo 
ga uopće ni upotrijebiti. 

U slučaju, da se od većeg broja pojedinačnih pokusa samo 
nekoliko njih ne slaže (a to možemo istražiti pomoću x 2 me¬ 
tode), treba ih odstraniti iz daljnjeg biometričkog proučavanja, 
kad ustanovimo, zbog čega su razlike nastale. Na pr. križanjem 
četiriju genotipa pšenice-golice s druga četiri genotipa brkulje, 
razvile su se u F 2 generaciji biljke bez osja i biljke s osjem u 
slijedećem omjeru: (Tab. LII.) 






Mjesto očekivanog omjera 3 : 1 (588,75 :196,25) omjer je 
fenotipova 3,42 : 1 (608 : 177). Dakle se on ne podudara s teo¬ 
retskim brojkama. 

Pomoću x 2 — metode možemo istražiti, da li je materijal 
homogen ili nije. 

Imamo podatke od četiri F 2 generacije (m) i za jedan par 
(2) svojstava (n). Prema tome je broj neovisnih skupina 
(<n—1). (n—1) = (4—1) . (2—1) = 3. Budući da je suma svih 
x 2 = 13,536 to je vjerojatnost P < 0,01 (Tab. LI.), dakle je ma¬ 
terijal heterogen. U tom se slučaju križanje broj 2 ne podudara 
s ostalima, jer je za x 2 = 10,34, vjerojatnost P = 0,01 —- 0,001 
(Tab. LI.). 

Dakle treba da drugo potomstvo izlučimo, a za ostala je 
tri potomstva onda omjer fenotipova: 524 : 167 t. j. 3,14 : 1. 

Istražimo li sad heterogenost, naći ćemo, da je broj neovi¬ 
snih skupina (m—1). (n—1) = (3—1). (2—1) = 2, dok je suma 
svih x 2 = 3,196 a vjerojatnost P > 0,2. Dakle pokusni je ma¬ 
terijal homogen i empirijski se podaci slažu s teoretskim. 

ad 2) Poremećen je u neovisnom cijepanju 
svojstava d i h i b r i d a u F 2 gen. 

Kod neovisnog cijepanja svojstava dihibrida očekujemo u 
F 2 generaciji cijepanje genotipova u omjeru 9 :3 :3 : 1. Križa¬ 
njem jednog jarog ječma s kapucom na klasu (furcatum) i go- 
log zrna, s nekim ozimim ječmom s osjem na klasu i obuvenim 
zrnom, dobili smo u F 2 generaciji 4 fenotipa. 

U tab. Lili. sakupljeni su brojčani podaci, eksperimentalni 
i teoretski, za omjer cijepanja svojstava: a) 9 :3 :3 : 1, te b) za 
k a p u c u na klasu s obuvenim i golim zrnom omjera 3 : 1 i 
osjem na klasu s obuvenim i golim zrnom omjera 3:1. 

Nadalje su u istoj tabeli i x 2 . 

Jesu li pokusi u skladu s teoretskom pretpostavkom? 

U omjeru 9 :3 :3 : 1 je x 2 = 14,954 (tab. Lili). Budući da 
ima četiri razreda različitih fenotipova, to su samo (4—1) = 3 
razreda neovisna. U Fisherovoj tabeli ćemo naći, da je za 3 ne¬ 
ovisne skupine i x 2 = 14,954 vjerojatnost P između 0,01 i 0,001. 

Uzmemo li u križanja graničnu vjerojatnost P — 0,05, 
to u našem primjeru eksperimentalni podaci nisu u skladu s teo¬ 
retskim. 

Biljke s kapucom manje su otporne protiv zime negoli 
biljke s osjem na klasu. U skupini biljaka s kapucom usporedimo 
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eksperimentalne brojke za golo i obuveno zrno s teoretskim 
omjerom 3 :1. Istoj ćemo provesti i sa biljkama s osjem. 

Broj neovisnih skupina je sada samo 2. Budući da x 2 = ; 

= 2,451, to je vjerojatnost P = između 0.5 i 0.3. Dakle se eks¬ 
perimentalni podaci s^da mnogo bolje slažu s teoretskim nego 
prije. 


Tab. LIH. 

a 

t 


Svojstvo 

eks. 

e 

teor. t 
9:3:3:1 

(e -t) 2 
t 

teor. t 
(3:1): (3:1) 

(e -t) 2 
t 

Furbatum 

obuveni 

286 

326,25 

4,966 

298,5 

0,523 

goli 

112 

108,75 

0,097 

99,5 

1,570 

Aristatum 

obuveni 

133 

108,75 

5,407 

136,5 

0,089 

goli 

49 

36,25 

4,484 

45,5 

0,269 

' 

Ukupno 


580 

580,00 

14,954= x a 

580,0 

2,451 = x» 


ad 3. Geni za 2 različita para svojstava loka¬ 
lizirani su u istom hromosomu. 

a) F 2 generacija: 

Križanjem kukuruza normalnog endosperma (Sh) i aleu- 
rona u boji (C) s kukuruzom deformiranog endosperma (sli) i 
bezbojnog aleurona izrasla su u F 2 generaciji četiri fenotipa u 
slijedećem omjeru: 

AB : Ab : aB : ab = 747 : 20 : 25 j 208 = 1000 = n 1 :n 2 :n 3 :n 4 = n 

Na prvi pogled vidimo, da se ovaj omjer fenotipova ne 
slaže s omjerom 9 :3 :3 : 1, koji dobijemo kod neovisnog cije¬ 
panja t. j. od neovisnih gena. 

Da li se empirijski podaci slažu s teoretskim o neovisnom 
cijepanju, istražit ćemo pomoću x 2 — metode. 

Ako ustanovimo, da se oni ne slažu, iako su inače homo¬ 
geni, onda mogu biti još tri uzroka nepodudaranja i to: 

a) poremećenje u cijepanju za strukturu endosperma Sh 
sh (A), 

f) poremećenje u cijepanju za boju aleurona Cc (B), 


y ) su vezani, dakle lokalizirani u istom hromosomu 

Bez obzira, da li su spomenuta dva gena za Sh i C lokali¬ 
zirana u istom hromosomu ili nisu, omjer je fenotipova, ako 
nema smetnja u jednom ili drugom paru svojstava ili u’ oba 
ovaj: 

(Sh C -f- Sh c) : (sh C + sh c) = (n 4 + n 2 ) : (n 3 + nj == 3 : 1 
(Sh C + sh C) : (Sh c + sh c) = (n 4 + n 3 ) : (n 2 + n 4 ) = 3 : 1 
U našem je slučaju za svojstvo Sh sh omjer fenotipova 
(Sh C + Sh c) : (sh C + sh c) = (747 + 20) : (25 + 208) = 3,29 :1 
a za svojstvo Cc: 

(Sh C + sh C) : (Sh c + sh c) = (747 + 25) : (20 + 208) = 3,38 : 1 
I za jedan i za drugi par svojstava ima razlika u cijepanju 
između empirijskih podataka i teoretskog omjera 3:1. Razlike 
doduše nisu velike, ali ne možemo olako prijeći preko njih, a 

^. a . s v e , ne uv J er * rno . J esu li smetnje u cijepanju varijaciono-sta- 
tistički opravdane ili nisu. Zato ćemo izračunati x 2 . 

Od spomenuta tri slučaja, koji mogu utjecati na nepodu¬ 
daranje empirijskih podataka s teoretskim, mogu u pokusu dje¬ 
lovati sva tri ili dva ili pak samo jedan. 

Svaki od ovih slučajeva može da uzrokuje jedan dio od 
cjelokupne razlike. Prema tome se x 2 ili bolje x 2 „ (n = ukupno) 
sastoji od ovih članova: 

Ako se brojčani podaci odnose na F 2 generaciju t. j. na 
križanje AaBb X AaBb kao što je to u našem primjeru, onda 
izračunamo x 2 po ovim formulama: 

x 2 a = (nj + n 2 - 3n s — 3n 4 ) 2 : 3n = [747 + 

+ 20 - (3 X 25) -(3 X 208)] 2 : 3 X 1000 = 1,512, 
iz toga proizlazi P > 0,05 

x 2 b = (Hi — 3n 2 n s — 3n 4 ) 2 : 3n = [747 - (3 X 20) + 

-(-25 (3 X 208)]*: 3 X 1000 = 2,5812 dok je P > 0,05 

x z v = (n, — 3n 2 3n 3 -(- 9n 4 ) 2 : 9n = [747 — (3 X 20) — 

(3 X 25) + (9 X 208)] 2 : 9 X 1000 = 685,584; x 2 A , x 2 B i x\ 
pripada po jedna neovisna skupina. U Fisherovoj tabeli LI. po¬ 
tražimo vjerojatnost P od gornjih brojeva. 

Na poremećaj u cijepanju jednog ili drugog para svojstava 
ili na vezane gene zaključujemo već onda, ako je P < 0,05. 

Budući da je u našem slučaju za strukturu endosperma i za 
boju aleurona vjerojatnost P > 0,05, to znači, da ova svojstva 
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nisu bila u cijepanju tako poremećena, da empirijski podaci ne 
bi odgovarali teoretskim. 

Ako ustanovimo, da nema smetnje u cijepanju za jedan i 
drugi par svojstava, a pokusni se podaci razlikuju od teoretskog 
omjera 9 :3 :3 :1, onda nam još treba izračunati i x 2 v , pa da 
ustanovimo, jesu li geni vezani ili nisu. U slučaju, da je P < 0,05, 
a ti našem je primjeru zaista tako, onda to znači, da su geni 
vezani. 

b) Povratno križanje: 

Imamo li podatke o povratnom križanju, za izračunavanje 
x 2 upotrebit ćemo druge formule. 

Primjer: 

. Po povratnom križanju kukuruza normalnog endosperma 
i aleurona u boji (Sh C . sh c) s kukuruzom deformiranog endo¬ 
sperma i bezbojnog aleurona, razvila su se četiri fenotipa u 
ovom omjeru: 


ShC 

Shc 

shC 

shc 

Ukupno 

418 

16 

14 

402 

850 

Hi 

n 2 

n 3 

"4 

n 


U tom je slučaju: 

x 2 a = (nj + n 2 — n 3 — n 4 ) 2 : n = (418 + 

+ 16 — 14 - 402) 2 : 850 = 0,3812 
i P > 0,05 

x 2 b — (n x — n 2 + n 3 — n 4 ) 2 : n = (418 — 16 + 

+ 14 — 402) 2 : 850 = 0,2306 
iz toga je P > 0,05 

x 2 v = (nj — n 2 — n 3 + n 4 ) 2 : n = (418 — 16 — 

— 14 + 402): 850 = 734,235 
dakle je P < 0,01. 

Dakle i iz povratnog križanja vidimo, da nema biometrički 
opravdanih poremećenja u cijepanju za strukturu endosperma 
(Sh sh) i' boju aleurona (Cc), jer je za Pa i za Pb > 0,05, 

Iz x 2 v = 734,235 je vjerojatnost P <0,05, to znači, da se 
geni (Sh sh) za strukturu endosperma i geni za boju (Cc) aleu¬ 
rona nalaze u istom hromosomu. 


C. - 2. BIOMETRIKA KORELACIJE KVALITATIVNIH 
ILI ALTERNATIVNIH SVOJSTAVA 

1. Korelacioni koeficije n t — r 

Kao kod kvantitativnih svojstava, tako i kod kvalitativnih 
međusobnu zavisnost istražujemo pomoću korelacionog koefi¬ 
cijenta i regresije. Na slijedećem primjeru izračunat ćemo ko¬ 
relacioni koeficijent za kvalitativna svojstva. 

Križali smo rajčicu visoka rasta i sferičnih plodova s raj¬ 
čicom niska rasta kruškolikih plodova zato, da u F 2 generaciji 
i daljnjim generacijama pronađemo po koju biljku niska rasta 
i sferičnih plodova. Osim toga htjeli smo istražiti, da li se vi¬ 
sina biljke i oblik ploda nasljeđuje ovisno ili neovisno jedno od 
drugoga. U F 2 generaciji razvile su se biljke ovih svojstava i u 
ovom broju: 

69 visoka rasta i sferičnih plodova 
8 visoka rasta i kruškolikih plodova 
7 niska rasta i sferičnih plodova 
19 niska rasta i kruškolikih plodova 

Je li visina biljke i oblik plodova u međusobnoj ovisnosti 
ili nije, naći ćemo tako, da izračunamo korelacioni koeficijent r. 
Na isti način kao kod kvantitativnih svojstava, treba i ovdje va¬ 
rijante poredati u korelacionu tabelu. Za suponovano svojstvo 
(x) uzet ćemo na pr. visinu biljke, a za relativno (y) svojstvo 
oblik ploda. 

Korelaciona tabela je ovakva: 


Tab. LIV. 




Oblik ploda y 




sferičan 

kruškolik 

Zbroj 



0y 

iy 


Rast 

visoki rast 0x 

69 = pl 

8 = pll 

Pox 

77 

X 

niski rast lx 

7 = pili 

19 = pIV 

PlX 

26 


Zbroj 

76 = p0y 

27 = P ly 

103 
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Korelacioni koeficijent izračunamo po Bravaisovoj formuli. 

Ako korelaciona tabela ima samo 4 rubrike, kako je to u 
našem slučaju, onda se za alternativna svojstva gornja formula 
promijeni ovako: 

PiPiv —PiiPiii 
V Pox' Pix ' Poy ’ Piy 

U našem je primjeru: 

(69 X 19) ~ (8 X 7) _ 1311-56 nrin 
; 1/77X26X76X27 1/4108104 

f Dakle između visine biljke i oblika ploda je jaka korelacija. 


2. Srednja pogreška korelacionog 
koeficijenta — m r 


Mogu li se iz izračunatog korelacionog koeficijenta izvo¬ 
diti pravilni zaključci, pronaći ćemo pomoću srednje pogreške 
takvog koeficijenta. 

Srednja pogreška korelacionog koeficijenta izračunava se 
i za kvalitativna svojstva po formuli, koju smo upotrebili za 


kvantitativna svojstva: m — 



U našem je primjeru m r = 


1 — 0,619* 
1/103 


= ° l61683 P = 0,0607. 
10,149 


Trostruka ova pogreška manja je od našeg korelacionog 
koeficijenta r = 0.619 tako, da je ta pogreška svakako realna 
t. j. signifikantna brojka za izvođenje pravilnih zaključaka. 


3. Regresija 

Regresija se i za alternativna svojstva izračunava po for 
mulama za kvantitativna svojstva, dakle Ry/x = r “> 

odnosno R x / y == — 

°y 
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4. Srednja pogreška regresije 


Ovu vrijednost izračunavamo po formulama, koje smo 
upotrebili za kvantitativna svojstva i to: 

m R = m r — za R x / y i po formuli: 

°y 

m R = m r za R y/X 
o x 

5. Asocijacioni koeficijent: 


Jesu li dva svojstva korelativno vezana ili nisu, možemo 
približno ustanoviti pomoću asocijacionog koeficijenta, koji je 
izradio engleski statističar Y u 1 e. 

^ , . . , pl plV — pII pIII 

Ovaj koeficijent glasi: As. k. = pl pIV + p n p m 

Za prijašnji primjer (str. 149) za visinu rasta i oblik plo¬ 
dova kod rajčice taj koeficijent je: 


As. k. = 


69X 19 — 8X7 1255 

69X 19 + 8 X 7 1367 


= 0,918 


Asocijacioni koeficijent kreće se u pozitivnom i negativ¬ 
nom smjeru od 0—1. Za mjerenje njegove jačine može nam i 
ovdje poslužiti R o e m e r—O r p h a 1 o v a tabela (Tab. 
XXXVIII., str. 103). 


C — 3. KORELACIJA IZMEĐU KVANTITATIVNIH 
I KVALITATIVNIH SVOJSTAVA 

Kad na istim organizmima istražujemo kvantitativna i kva¬ 
litativna svojstva, katkada je potrebno znati, jesu li ta svojstva 
međusobno ovisna, pa ako jesu, kako je jaka ta ovisnost. Ra¬ 
čunski postupak, kad se radi o istraživanju korelacije između 
kvalitativnih i kvantitativnih svojstava, prikazat ćemo na ovom 
primjeru: 

Križali smo jednu čistu liniju graška, koji ima veliko žuto 
sjeme, s drugom linijom koja ima malo i zeleno sjeme. Na bilj¬ 
kama F 2 generacije istražili smo boju i težinu sjemena, da pro¬ 
učimo, ima li korelacije između tih svojstava. Uz pomoć kore- 
lacione tabele izračunali smo korelacioni koeficijent. 
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Boju sjemena sam uzeo u korelacionoj tabeli za suponovano 
(x) svojstvo, a težinu sjemena za relativno (y) svojstvo. 

Prema tome korelaciona tabela je ovakva: (Tab. LV.) 


Tab LV. 



A y 


Korelacioni koeficijent izračunamo i ovdje po formuli: 

r=rr £ P a * a y — n b x b y 
n p x X Py 

Stoga se najprije mora izabrati polazna točka A x i A y 
(aproksimativne srednje vrijednosti) za svojstvo x i y. Uzmimo 
na pr. za Ax=Ox žutu boju i za Ay = 25,5 cg. 

b x izračunamo po formuli za kvalitativna svojstva b x = ^ = 

125 

~ 448 = 0,279, by pak po formuli za kvantitativna svojstva 

b y =^B5 = -^ = 0,116. Kako je A x = O x to su produkti 

p • a x • a y —za-žotu boju i težinu zrna ravni ništici, a isto tako i 
produkti u vertikalnom redu od A y . 

Prema tome preostaju samo produkti pa x a y za zelenu 
boju, a izračunaju se ovako: 

Pozitivni produkti su oni desno od A y , a negativni oni 
lijevo od A y , dakle: 


Pozitivni 

produkti: 


^ P ■ 3x . a y 


17. 

1 

. 1 

17 

13. 

1 

. 2 

26 

19. 

1 

. 3 

27 

5. 

1 . 

. 5 

M 

1 . 

1 . 

, 5 

5 


Ukupno: 95' 


Negativni 


produkti: 
P • 3x • a y 


19. 1 . 1 

19 

20. 1 . 2 

40 

13. 1.3 

39 

1.1.4 

4 

1.1.5 

5 


Ukupno: — 107 


2pa x a y = — 107 +• 95 = — 12. 

px izračunamo po formuli za standardnu devijaciju kvali¬ 
tativnih svojstava: 


,-j _1 VPo Pi 1/323 X 125 

n 448 

: 


± 0,448 = p x 


Formula za p y glasi: p y = _ b 2 _ 

1 n v 

f l844 

“448 (0.116) 2 = 2,026 

= ?P fl x a y —nb x b y _ ~J2~ 448 X 0,279 X 0,116 
n P* Py 448~X 0,448 _ X'2,026 = 


A . . J n.ucucijem r = u,Učsa sliiedi, da 

u b ° je 1 tezine s J' em ena graška nema nikakve korelacije 
(Vidi tab. XXXVIII., str. 103). J 





D. NEKE BIOMETRIČKE OZNAKE I FORMULE 

a) OZNAKE 

A . . . polazna točka - aproksimativna srednja vrijednost 
(str. 22) 

a • • • areal (str. 22) 

As. k.. asocijacioni koeficijent (str. 151.) 

b ' ' 1 uđa, i enos t srednje vrijednosti M od A izražena u area¬ 
lima. A + b a = M (str. 22.) 

D . . . razlika - diferencija. Udaljenost jedne varijante V od M 
D — V — M (str. 64.) 

E . . . eksces (str. 53.) 

E. P. . vjerojatna pogreška (str. 39.) 

F L , F a itd. filijalne generacije kod bastarda (str. 127.) 

K . . . kombinacioni broj za F 2 generaicju (str. 131.) 

M . . . srednja vrijednost (str. 21.) 

m • • • srednja pogreška (str. 36.) 

m D • • srednja pogreška diferencije (str. 72.) 

m s • • srednja pogreška sume (str. 72.) 

ffl " ' ’ srednja pogreška srednje vrijednosti M (str. 36.) 

m v . • srednja pogreška varijacionog koeficijenta v (str. 38.) 

”1° • • srednja pogreška standardne devijacije (str. 38.) 

iM r . . srednja pogreška korelacionog koeficijenta r (str. 107) 

niR . . srednaj pogreška regresije (str. 109.) 

Med . Medijana (str. 20.) 

Mo . . Modus (str. 55.) 

n . . . broj varijanata (individua, parcela, slučajeva i t. d.) 

(str. 12.) 

P . . . probabilitet (str. 139.) 
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p . . . frekvencija (broj individua, parcelica, slučajeva i t. d.) 

u pojedinim razredima (str. 21.) 

Q . . . kvartil (str. 28.) 

q, . . . apsolutna vrijednost zadnje varijante u prvoj četvrtini 
.varijacionog reda (str. 28.) 

q 2 . . . apsolutna vrijednost zadnje varijante u drugoj četvrtini 
varijacionog reda (str. 21.) 

q 3 . . . apsolutna vrijednost zadnje varijante u trećoj četvrtini 
varijacionog reda (str. 28.) 

q . . . apsolutni broj individua jedne alternative (str. 135.) 

R . , . regresija (str. 108.) 
r . . . korelacioni koeficijent (str. 100.) 

S ... . . asimetrija krivulje (str. 50.) 

2 . . . znak zbrajanja (str. 21.) 

a. . . . standardna devijacija (str. 29.) 

p . . . . (str. 100.) 

V . . . varijanta (str. 21.) 

Vs • • • varijaciona širina (str. 28.) 

v ... varijacioni koeficijent (str. 34.) 

x s . . . x 2 metoda za izračunavanje probaliteta P (str. 139.) 


b) FORMULE 


A) Kvantitativna svojstva: 


M e d i a n a: str. 20. 

Med = q 2 

Srednja vrijednost: str. 21. 

M = A + B = A + ab = A + a^ 

n 

Varijaciona širina: str. 28. 

V S = V„-V, 

Kvartil: str. 28. 

n — ~ 9i ■ 

2 

Standardna devijacija: str. 33. 


-±.V^ 

n 


^P£!_ h 2 ili 


b 3 ili točnije 


o= ±a b 2 

n — 1 


/ z pa* a 52 

a = ± a —- b2 —12 Sheppardova formula 

Kvartilni koeficijent: str. 34. 

M 

V a r i j a c i o ni koeficijent: str. 34. 


Modus: str. 55. 

Mo = 3 Med — 2 M 

Asimetrija: str. 50. 


s = (ll^_3b^L a +2 b3): 
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Eksces: str. 53. * 

E = [(^rP~ ~ 4b ^7T" + — 3b4 ) : p 4 ] — 3 

Korelacioni koeficijent: str. 100. 

a) Bravaisov: str. 100. 

jg»pa x a y — n b x b y 

n p x Py 

za mali broj varijanata: str. 106. 

r SV x V y -nM x M y 

II Oj Oy 

b) parcijalni korelacioni koeficijent: 
a) 3 različita korelirana svojstva: str. 109 i 110. 
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Vjerojatna pogreška — E: 

a) srednje vrijednosti: str. 39. 

EPm = 0,6745 ihm 

b) standardne devijacije: str. 39. 

EP 0 = 0,6745 m 

• c) varijacionog koeficijenta: str. 39. 

EP V = 0,6745 m v 

d) Peterova vjerojatna pogreška: str. 63. 
a) za jednu varijantu sg: 

EP sg = 0,8453 


S(dzD) 

/ n (n — 1) 




i6a 


j>) za srednju vrijednost M: 

EP„ = 0,8453 

nj/n —1 

e) Besselova vjerojatna pogreška: str. 63. 

a) za jednu varijantu sg: 

EPs g = 0,6745 VJLOj 

/>) za srednju vrijednost M: 

EP„ = 0,6745 ~\[ SD2 
r n (n — 1) 

f) korelacionog koeficijenta: str. 

EP r = 0,6745 m r 

g) regresije: str. 

EPr = 0,6745 m R 

Grafički prikaz korelacije: str. 112. i 113. 

1. Linearna podjela varijanata: y — a + bx 

a) parni broj podataka: 

b=^[2S(xy)-(n + l)(S(y)J 

a = -I[2(y)-2(n)b] 

b) neparni broj podataka: 

b = ji[2(xy)-k 2 2(y)] 

a = js(y)-k s b 

2. Parabolična podjela varijanata: y = a + bx + cx 2 

c = i- [« 2 (y) - (3 2 (xy) + ? 2 (x 2 y)] 

K 3 

b = ± [2 2 (xy) - (n + 1) 2 (y)] - (n+1) • c (parni 

broj podataka) 

b=^[2 (xy) — k, 2 (y)J — (n + 1) ■ c (neparni 


broj podataka) 
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a = -j-[2(y)- 2(n)b-2(n 2 ) c] 


B) Kvalitativna svojstva: 
Srednja vrijednost: str. 131. 


Standardna devijacija za srednju vrijednost: str. 132. i 133. 

m 

n 

100VSS Hi 
n 

0 = ± \ Po% ' Pi°/o 

Standardna devijacija kombinacionog broja: str. 134. 

3 ■ K Vp^ 

n 

Srednja pogreška m: 

a) srednje vrijednosti: str. 134. 

m M = 

b) kombinacionog broja: str. 135. 

_ 1/N(K - N) 


c) apsolutnog broja: str. 135. 


/ r q (n - 

n 


Korelacioni koeficijent: str. 150. 

r = P 1 pIV pH pili 

V Po* Pix Poy Piy 

Asocijacioni koeficijent: str. 151. 

As k _ Pl PlV ~ PH PIH 
' ‘ pl pIV + pll pili 

11 Tavear: Biometrika u poljoprivredi 


LITERATURA 




Feldman W. M.: Biomathematic, London, 1923. 

Fisher, R. A.: Statistical methods for research workers. 
London, 1938. 

Fisher R. A. and Y-ates F.: Statistical tables. Edinburgh, 1938. 

Isenbeck K.: Die Priifung des Zuchterfolges. Roemer u. Ru- 
dorf: Handbuch der Pflanzenziichtung. Bd. I., Berlin 1941. 

Jastremski, Hotimskii, Stavrovskii i Bojarskii: Teorija ma- 
tematičeskoj statistiki. Planhozgiz, 1931. 

Janko J.: Zahlady statističke indukce. Praha 1937. (Nakl. 
St. uradu statistickeho). 

Johannsen W.: Elemente der exakten Erblichkeitslehre, 
Jena 1926. 

Just G.: Praktische Obungen zur Vererbungslehre. Frei- 
burg 1923. 

Kapteyn J. C.: Ske\v frequency curves in biology and sta- 
tistics, Groningen 1903. 

Kapteyfl J. C.: Definition of the Correlation Coefficient 
(Monthly Notices of R. Art Society, 1912.) 

Kohn S.: Zaklady teorie statističke metody. Praha 1929. 
(Nakladem St. uradu statistickeho). 

Koller S.: Graphische Tafeln zur Beurteilung statistischer 
Zahlen. Dresden u. Leipzig 1943. 

Kronacher C. u. Pataw C.: Biometrik, Berlin 1930. 

Lang A.: Die experimentelle Vererbungslehre in der Zoo- 
logie seit 1900. Jena 1914. 

Mikič F.: Smjer gibanja (trend); tablice za računanje kod 
ekvidistantnih apscisa (rukopis). 

Patau K.: Eine neue X 2 Tafel. Zeitschrift ftir induktive Ab- 
stammungs-u. Vererbungslehre. Bd. 80. Berlin 1942. 

Pearson R.: Tables for Statisticians and Biometricians, Cam- 
bridge 1914. 

Romanovskii: Elementi teorii korreljacii. Taškent 1928. 

Qućtelet: »Anthropomćtric«, Pariš 1871. 

Sapeg'in A.: Variacionnaja statistika, 1929. 

* 




164 


Sinnot E. and W. Dunn L. C.; Principle of genetics. New 
York and London 1939. 

»Student«: The probable error of a mean. Biometrika. 

VI. 1908. 

Tavčar A.: Varijaciona statistika u eksperimentalnoj poljo¬ 
privredi. Izdanje Ministarstva poljoprivrede. Beograd 1929 
Tavčar A.: Temelji nauke o nasljedstvu. Izdanje Ministar¬ 
stva poljoprivrede, Beograd 1930. 

Tavčar A.: Caklavost i njen odnos napram diugim svoj¬ 
stvima pšeničnog zrna. Poljodjelska znanstvena smotra br. 4, 

Zagreb 1941. . 

Tavčar A.: Dvostruko-recesivno nasljeđivanje opadanja osja 

kod ječma. Poljodjelska znanstvena smotra br. 6 - Za S reb 

Tedin O.: Biologische Statistik. Roemer u. Rudorf. Hand- 

buch der Pflanzenzuchtung. Bd. I. Berlin 1941. 

Udny Yule G.: An Introduction to the Theory of Statistics, 

L ° n< Udny Yule G.: On the Theory of Correlation. Journ. Royal 

Stat ZMer W?ForS'und Tabellen zur Errechnung des mit- 

^^^Biometriba? a'Journal for the S.udy of Biological ProMems. 


SADRŽAJ 


Strana 

5 

7 


Predgovor 

UVOD 


. . 9 

A. VARIJAB1LITET 9 

a) Općenito . • •.' ' 14 

b) Svojstva varijanata . ■ • jg 

c) Valjanost biometrickih istraživanja • 


B. B10METR1CKA ISTRAŽIVANJA KVANTITATIVNIH SVOJSTAVA 

B-1. Biometrika varijacije 

a) Varijacioni red • • .. ■ • 

b) Grafićki prikaz podjele varijanata . . ■ • 

Biometričke vrijednosti: 

I. Apsolutna vrijednost svojstva 

1. Medijana — Med 

2. Srednja vrijednost — M . • • ' 

a) Srednja vrijednost razrednih varijanata . 

rharlierova metoda kontrole • • • .. • * 

b) Srednja vrijednost cijelih (diskretnih) varijanata 

II. Računanje varijabiliteta 

1. Varijaciona širina — .. 

2. Kvartil — Q . 

3 Standardna devijacija — 0 ■ ■ ■ ' ' 

a) Standardna devijacija razrednih varijanata . 

b) Standardna devijacija cijelih (disloetnih) vari 

janata . . 

4. Varijacioni koeficijent — v .. 

5. Kvartilni koeficijent — .. 

III. Srednja i vjerojatna pogreška m i E. .. 

1 Srednia pogreška — m . • ■ • • ■ 

a) Srednia pogreška srednje vrijednosti . • • 

a) Gianice za biometrički opravdanu sred j 

b) Siednja pogreška standardne devij^z7^" 

c) Srednja pogreška vanjacionog koeficijenta 

2. Vjerojatna pogreška — E. .. 

IV. Poredenje varijacionih redova sa teoretskom raspodjelom 
(distribucijom) 

1. Poredenje brojaka 

2. Grafički prikaz 

3. Hazenov dijagram . ■ 


17 

18 


19 

20 
21 
23 

25 

26 


28 

28 

29 

31 

32 
34 
34 

34 

36 

36 

36 

38 

38 

39 


40 

41 
43 
47 








Strana 


V. Abnormalne varijacione krivulje .... 48 

1. Asimetrične krivulje genotipa.49 

2. Ekscesivne krivulje.52 

3. Varijacione krivulje sa dva ili više vrhova ... 55 

VI. Modus ..55 

VII. Neke specijalne varijaciono-statističke metode u poljo¬ 

privredi 55 

1. Kako se može unaprijed izračunati broj varijanata 

potrebnih za pokus.56 

2. Računsko izjednačenje varijacije priroda izazvane 

heterogenošću tla.58 

a) Standardna metoda.59 

b) Kako se kod sortnih pokusa s okopavinama 

računski korigira prirod modificiran uslijed 
praznih mjesta.62 

3. Peterova i Besselova metoda za izračunavanje vje¬ 

rojatne pogreške.. . 62 

4. Modificirana Peterova metoda za srednju pogrešku 65 

5. Studentova metoda . 67 

6 . Mjerenje točnosti neke razlike pomoću Kollerove 

tabele i dijagrama.71 

7. Srednja i vjerojatna pogreška diferencije i sume . 72 

8 . Varijanca.77 

a) Općenito .77 

b) Primjena Fisherove metode na nekvađratnoj 

podjeli pokusnih parcelica.82 

c) Kvadratna podjela pokusnih parcelica ... 90 

B-2. Biometrika korelacije 

Općenito . 95 

Linearna korelacija.99 

1. Bravaisov korelacioni koeficijent za mnogo varijanata 99 

2. Izračunavanje korelacionog koeficijenta za mali 

broj varijanata . 106 

3. Srednja i. vjerojatna pogreška korelacionog koefi¬ 

cijenta .107 

4. Regresija.108 

5. Srednja i vjerojatna pogreška regresije ■ . . 109 

6 . Korelacija između više svojstava — parcijalna ko¬ 

relacija ..109 

Nelinearna korelacija 111 

7. Grafički prikaz korelacije.112 

a) Linearna korelacija . 113 

b) Parabolična korelacija.120 

C. BIOMETRIČKA ISTRAŽIVANJA KVALITATIVNIH ILI ALTER¬ 
NATIVNIH SVOJSTAVA 

C-l. Biometrika varijacije 

Općenito . 127 

I. Dvije alternative.130 

1. Srednja vrijednost.130 

2. Standardna devijacija .132 

3. Srednja pogreška srednje vrijednosti .... 133 


I 



Ii. Više alternativa 


Strana 
. 136 


III. Vjerojatnost — »probabilitet« x- — metoda 

1. Istraživanje pokusnih podataka na heterogenost 

2. Poremećenje u neovisnom cijepanju svojstvava dihi- 

brida u Fa gen. 

3. Geni za dva različita para svojstava lokalizirani u 

istom hromosomu. 

a) Fa generacija.. . . 

b) Povratno križanje . -. 

C-2. Biometrika korelacije kvalitativnih ili alternativnih 
svojstava.... 

1. Korelacioni koeficijent — r. 

2. Srednja pogreška korelacionog koeficijenta — m r 

3. Regresija. 

4. Srednja pogreška regresije. 

5. Asocijacioni koeficijent. 


139 

143 

145 

146 
146 
148 


149 

149 

150 

150 

151 
151 


03. Korelacija između kvantitativnih i kvalitativnih svojstava 151 


D. NEKE BIOMETRIČKE OZNAKE 1 FORMULE 

Literatura . 

Tabelice za l/n . 


155 

163 

169 






































TABLICE ZA j/n * 


F i s h e r, R. A, anu Y a t e s, F.: Statistica! tables. Edinburgh, 1938 


Uputa za upotrebu tablica 


Zadatak: Treba izvaditi drugi korijen: 

]/39a ]/3Q, j/395; j/393 j/3953 j/39,58. 

Rješenje: U tablici lijevo od vertikalne linije potražimo 
najprije brojku 39. U horizontalnom smjeru od ove brojke su 
dva reda brojeva, prvi za neparni, a drugi za parni broj mjesta 
lijevo odnosno desno (ako je brojka manja od 1) od decimalne 
točke broja iz kojega vadimo drugi korijen. 

U glavi tablice su u horizontalnom smjeru najprije brojke 
od 0 do 9, a zatim od 1 do 5. 

Budući da j/390 ima neparni broj mjesta, to potražimo u 
rubrici pod 0 u vertikalnom smjeru kod 39 za j/3y0 = 19,749, 
a za parni broj mjesta j/39 = 6,245. 
j/395 = 19,875 (vertikalna rubrika 5), 
j/39,5 = 6,2849 (vertikalna rubrika 5). 


j/395 . 19,875 

za 0,8 o ,5 . 13 

_ desni dio 

J/ 395,8 tablice 

(rubrike 5 i 3) 0,3.8 

J/3953 = 19,896 

J/393”.6.2849 

za 0,08 0,05 _... 40 

_ desni dio 

j/39,58 tablice 

_ .... A no O A 


(rubrike 5 i 3) (0-03.24 

]/39/3 = 6,2913 
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